Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,947 G→T 100% R418L (CGT→CTT)  aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,9470GT100.0% 155.7 / ‑6.7 50R418L (CGT→CTT) aceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (26/23);  minor base A (1/0);  total (27/23)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAACGTGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCA  >  NC_000913/126802‑127084
                                                                                                                                                 |                                                                                                                                         
aaaCGTGCAGAAGCAGCTCCGGCAGCGGCTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTaa                                                                                                                                        >  2:756669/1‑149 (MQ=255)
  aCGTGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGGCTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTCGTGCGAACCAGAGGNATAACt                                                                                                                                      >  2:23445/1‑149 (MQ=255)
    gTGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTgg                                                                                                                                    >  2:213726/1‑149 (MQ=255)
       cAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTa                                                                                                                                 <  1:132197/149‑1 (MQ=255)
         gAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAAt                                                                                                                               <  2:672377/149‑1 (MQ=255)
             cagcTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATc                                                                                                                           <  2:529709/149‑1 (MQ=255)
              agcTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGCAATCGAAGAAGTGGAACTGGGCCGCATCCAGAGAATCTCTGGTGCGAACCTGAGCCTTAACGGGGGAATGATcc                                                                                                                          >  2:289309/1‑149 (MQ=255)
                 tCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAAGTAGGTAATGATCCCGc                                                                                                                       >  1:370426/1‑149 (MQ=255)
                 tCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGc                                                                                                                       >  1:564983/1‑149 (MQ=255)
                   cGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCAt                                                                                                                     <  1:657395/149‑1 (MQ=255)
                   cGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCAt                                                                                                                     <  1:649597/149‑1 (MQ=255)
                         cgccTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACt                                                                                                               <  2:522235/146‑1 (MQ=255)
                         cGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACt                                                                                                               <  1:516710/149‑1 (MQ=255)
                          gACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTc                                                                                                              >  1:659970/1‑149 (MQ=255)
                          gACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTc                                                                                                              >  1:278333/1‑149 (MQ=255)
                               gCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACtcc                                                                                                         >  2:498957/1‑147 (MQ=255)
                               gCGGTGGTATCCCTGGCATGCTGCCGTGGACGAAGGTGGACTTCAGCAAGTTTCGTGAAATAGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGGAATGATCCCGCATGTTACTCACTTc                                                                                                         >  1:478895/1‑149 (MQ=255)
                                 ggtggtATCCCTGGCATTCAACCAAAGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCACGTTACTCACTCCGa                                                                                                       >  2:319141/1‑149 (MQ=255)
                                 ggtggtATCCCTGGCATGCTGCCGGGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGa                                                                                                       <  2:457753/149‑1 (MQ=255)
                                     gtATCCCTGGCATGCTGCTGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACaaa                                                                                                   <  1:747504/149‑1 (MQ=255)
                                          ccTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGa                                                                                              >  2:750399/1‑149 (MQ=255)
                                              gCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATc                                                                                          >  2:376694/1‑149 (MQ=255)
                                                  gctgcCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCg                                                                                      <  1:152104/149‑1 (MQ=255)
                                                   ctgcCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGa                                                                                     <  1:213726/149‑1 (MQ=255)
                                                       cGTGGCCGGAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTCCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGACATCACCGAGTTg                                                                                 <  1:563283/149‑1 (MQ=255)
                                                       cGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCCACAAAACCGATAGCACCACGTTg                                                                                 >  2:551396/1‑149 (MQ=255)
                                                        gTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTgg                                                                                <  2:489374/149‑1 (MQ=255)
                                                              gAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATCTCACCGAGTTCGAGGCGt                                                                          >  1:288005/1‑149 (MQ=255)
                                                              gAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATAGCACCGAGTTGGAAGCgg                                                                          >  2:12959/1‑148 (MQ=255)
                                                               aaGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGtt                                                                         <  1:449894/149‑1 (MQ=255)
                                                                  gTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCg                                                                      >  1:525578/1‑149 (MQ=255)
                                                                        ttCAGCAAGTTTGGCGAAATCGAAGAAGTGAAACTAGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAAca                                                                >  2:631987/1‑149 (MQ=255)
                                                                          cAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAAcagc                                                              >  1:501591/1‑149 (MQ=255)
                                                                            gCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTGCCGTAAAcagcag                                                            >  1:411906/1‑149 (MQ=255)
                                                                              tactttGGTGAAACCGAAGAAGTGGAACGGGGCCGCACCCAGAAAATCTCTGGGGCGAACCTGAGCCTTAACTGGGTAATGTGCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGaa                                                          <  2:23148/146‑1 (MQ=255)
                                                                                      tGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGaagaag                                                  <  1:713369/149‑1 (MQ=255)
                                                                                         aaTCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAgcgg                                               <  1:750399/149‑1 (MQ=255)
                                                                                                       aaCTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTAACTTCGACAAACCCGATATCACCGAGTCGGAACCGTTCCGAAAGCACCAGAACAAAGAGACGTCGAAACGAAAGCTg                                 >  2:308920/1‑149 (MQ=255)
                                                                                                         cTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGa                               <  2:564983/149‑1 (MQ=255)
                                                                                                          tGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGAt                              <  2:525578/149‑1 (MQ=255)
                                                                                                            ggCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGAtgt                            <  1:464661/149‑1 (MQ=255)
                                                                                                              ccGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGa                          >  2:470971/1‑149 (MQ=255)
                                                                                                                 cATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACGGCGGCAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGa                       >  2:200369/1‑149 (MQ=255)
                                                                                                                      aGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCGGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCCGAACGTAGAAGCGGCGAAACGTACGCTGCATGTGAAGATCAcc                  >  2:476710/1‑149 (MQ=255)
                                                                                                                       gAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCAccc                 >  1:247843/1‑149 (MQ=255)
                                                                                                                         aaaTCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCg               <  1:716661/149‑1 (MQ=255)
                                                                                                                                tGGTGAGAACCTGAGCCTTAACTGGGTAGTGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCt        >  2:718380/1‑149 (MQ=255)
                                                                                                                                   tGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTtca     <  2:625440/149‑1 (MQ=255)
                                                                                                                                    gCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTtcat    <  1:691511/149‑1 (MQ=255)
                                                                                                                                      gAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAACCAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTtcatca  >  1:413866/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
AAACGTGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCA  >  NC_000913/126802‑127084

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: