Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 125,387 C→T 100% Q791* (CAG→TAG)  aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,3870CT100.0% 137.7 / ‑6.7 48Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/22);  minor base A (1/0);  total (26/22)
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.

ATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGA  >  NC_000913/125249‑125525
                                                                                                                                          |                                                                                                                                          
aTCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGTGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATAGCTTAGGt                                                                                                                                        >  1:836144/1‑143 (MQ=255)
aTCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGt                                                                                                                                        <  2:477526/143‑1 (MQ=255)
  ccTGTCGAAACATTTTGTCGACGGTTCTAACCTTCATAGCGTTGTTTCCGGTACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtga                                                                                                                                      <  1:121726/143‑1 (MQ=255)
                 cGGCGTAGGTTCGGCCGTTGATAGCGTGCCCTCCTGCACCGAGCTGGCGCGTGATGGTCAGGATTGTGCATGCTGGAACATGCCGCACCCGCTGGAAACGCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCgg                                                                                                                       <  2:872341/143‑1 (MQ=255)
                      tAGGTTCTGACGTTTATCGCGTGACCACCTTCACCGAGCTGGCGCTTGATGGTCAGGATTGTCAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTg                                                                                                                  <  2:145829/143‑1 (MQ=255)
                      tAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAggg                                                                                                                  >  2:171660/1‑141 (MQ=255)
                       aGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTgg                                                                                                                 >  1:354809/1‑143 (MQ=255)
                               gggTTTATGGCGTGACCTCCTTAACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAcc                                                                                                         <  2:474377/141‑1 (MQ=255)
                                 gTTTATAGCGTGACCTCCTGCACCGAGGTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGa                                                                                                       <  1:505258/143‑1 (MQ=255)
                                   ttATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCGTAGGTGATTAACGACGCTCCGGCAGTGGCGTCGACCGACt                                                                                                     >  2:686541/1‑143 (MQ=255)
                                      taGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGAGTGTGAACGCTGGAACATCGTGCACCCGGTGGAAACTCCGCGCGTTGCGTGTACCACTTCGGGGATGAACGACGCTCCGGCAGTGGCATCTCCCGACcata                                                                                                  >  2:696496/1‑143 (MQ=255)
                                      taGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGAATGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCCTCTACCGACtata                                                                                                  >  1:8978/1‑143 (MQ=255)
                                          gTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGaa                                                                                              >  1:772863/1‑143 (MQ=255)
                                                  ctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                                                                      <  1:581015/143‑1 (MQ=255)
                                                      aCCGAGCTGGCGCTTAATGGTTATGTTTGTGATCTCTGTAACATGCGGTACTCGCTGGAATCTCCGCGCGTTCCGCATAGCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGa                                                                                  <  1:123253/143‑1 (MQ=255)
                                                        cgagtcggCGCGTGATGGTCATGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGc                                                                                <  1:518670/137‑1 (MQ=255)
                                                         gAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCa                                                                               >  2:243744/1‑143 (MQ=255)
                                                          aGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAg                                                                              <  1:711999/143‑1 (MQ=255)
                                                                cgcgTGATGGTCAGGATTGTGACCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGt                                                                        <  1:571041/143‑1 (MQ=255)
                                                                 gcgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCGGCGCGTTCCGTATATCGCTAAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTGTATGAAACTGCCCGCCGAGAGGGTCCGTa                                                                       >  2:813773/1‑143 (MQ=255)
                                                                  cgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTAc                                                                      >  2:808492/1‑143 (MQ=255)
                                                                      aTGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTAc                                                                  <  1:471993/143‑1 (MQ=255)
                                                                        ggTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACAGACTATATGAAACTGTTCGCTGAGCAAGTCCGTAATTACgg                                                                >  2:224259/1‑142 (MQ=255)
                                                                         gTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGCACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTAGTTACGTa                                                               >  2:553024/1‑143 (MQ=255)
                                                                         gTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTa                                                               <  2:72897/143‑1 (MQ=255)
                                                                          tCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAc                                                              >  1:488197/1‑143 (MQ=255)
                                                                             ggATTGTGAACGCTGGAACATGCTGCACCCGCTGGAACCTCCGCGCGTTCCGTATATCGCTTCGGATGTGAACGACACTCCGGTAGAGTTTTGTGCCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCgg                                                           <  1:206482/143‑1 (MQ=255)
                                                                                 tgtgggcGCTGGAACATGCGGTACGCGGTGGAAGCGTCGCGCGTTCCGTATATCGCTTAGGTGATGAACGAGGCTCCGGCAGTGGCATCTACCGAGTATATGAAGCTGTTCGCTGAGCAGGTCCGTAGTTACGTACCGGCTga                                                       <  2:752846/137‑1 (MQ=255)
                                                                                     aaCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTgacgac                                                   >  2:130067/1‑143 (MQ=255)
                                                                                         cTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCTGGTCCGTCCTTCCGTTCCGGCTGCCGACTAcc                                               >  1:307327/1‑143 (MQ=255)
                                                                                                   ctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTggg                                     <  1:386723/143‑1 (MQ=255)
                                                                                                    tgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGACAGAAGCATCTACCGACTATATGCAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTCCCGCGTACTCGGt                                    >  1:126756/1‑143 (MQ=255)
                                                                                                      cACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTAc                                  <  2:488197/143‑1 (MQ=255)
                                                                                                         ccGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTCCCGCGTACTGCGGACTGa                               >  2:752305/1‑143 (MQ=255)
                                                                                                            cTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATgg                            <  2:125656/143‑1 (MQ=255)
                                                                                                               gAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGcgt                         >  1:264851/1‑141 (MQ=255)
                                                                                                                   cTCCGCGCGTTCCGTATATCGCTTAGTTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGt                     <  2:710555/143‑1 (MQ=255)
                                                                                                                    tCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTc                    >  1:515712/1‑143 (MQ=255)
                                                                                                                    tCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACGTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGGTGGCATCGGTc                    >  1:107327/1‑143 (MQ=255)
                                                                                                                          cgTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATATACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGTCTGACGACTACCGCGTACTGGGGACTGATGGCTTCGGTCGTTcct              >  1:854899/1‑142 (MQ=255)
                                                                                                                              ccGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAg          <  1:130067/143‑1 (MQ=255)
                                                                                                                              ccGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGGCGCTCCGGCAg          >  1:489859/1‑143 (MQ=255)
                                                                                                                              ccGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACCGATGGCCTCGGCCGGTCCGAccg          >  1:491835/1‑141 (MQ=255)
                                                                                                                               cGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGc         >  1:131479/1‑143 (MQ=255)
                                                                                                                                 tataTCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCg       >  2:64021/1‑143 (MQ=255)
                                                                                                                                  atatCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGt      <  1:606733/143‑1 (MQ=255)
                                                                                                                                     tCGCTTAGGTGATGAACCACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTCCCGTCTTACTTTTTTTTTTTACTGTGTACTGATGGCTTCGGTCGTTCCGACAGCCGTgag   <  1:168956/143‑1 (MQ=255)
                                                                                                                                      cGCTTAGGTGATGAACGACGCTCCGGCAATGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTgaga  >  1:418451/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                          
ATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGA  >  NC_000913/125249‑125525

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