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% 92.3 / NA 36Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/23);  total (13/23)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAG  >  NC_000913/125252‑125524
                                                                                                                                       |                                                                                                                                         
cTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACATCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCCCTTAgtttat                                                                                                                                    >  2:497283/1‑138 (MQ=255)
   gccaaaGCTTACGGCGTAGGTCCGGACGTTTGGAGGTCGACCTCCTTCAACGACCTGGCGCGTGATGGGCAGGATTGCGACCGCTGGAACAGGCTGCCCCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGaa                                                                                                                                 <  2:338233/140‑1 (MQ=255)
    cGAAAGATTACGGCGTAGGTTCTGACGTTTCGCGTCGCACTTCTTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAac                                                                                                                                <  1:366047/143‑1 (MQ=255)
     gAAAGGTTACGGAGTAGGTTCTGTCGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacg                                                                                                                               <  1:826846/143‑1 (MQ=255)
          aTTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGCACGACGCt                                                                                                                          >  1:332928/1‑143 (MQ=255)
                       ttCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGAGGCTCAGGCAGTGGcgt                                                                                                             >  1:13878/1‑141 (MQ=255)
                        ccggatgTTTATAGCGTGATCTACTTGAACGAGTTGGCGTGTGATGGTTAGGCTTGTGATTGATGGATCGTGTATAACCCGCTGGAAACTCCGCGCGTTCAGTCAATCATTTATGTGATGAACGACGCTCCGGCAGTGGCATc                                                                                                            <  2:430712/137‑1 (MQ=255)
                            ctggttATATCGTGACCTCCGTCACCGGGCTGGCGCTTGATGGTCAGTATTGTGATCGCTGGAACATGCTGCTCCCGCTGGAAACGCCGCGCGGTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAcc                                                                                                        <  2:215441/139‑1 (MQ=255)
                            aCGTTGATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCCGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAcc                                                                                                        <  1:156343/143‑1 (MQ=255)
                              gTTTATCGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACACGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGa                                                                                                      <  1:488607/143‑1 (MQ=255)
                              gTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGa                                                                                                      <  1:676459/143‑1 (MQ=255)
                                 tataGCGTGACCTCCTTCACCGGGTTGGGTCGCGGTGGTCCGGATTGTGCACGCTGGGACATGCTGCACCCGCTGGAAACGCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACta                                                                                                   <  2:869423/143‑1 (MQ=255)
                                   taGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtata                                                                                                 <  2:695275/143‑1 (MQ=255)
                                        tGACCTACTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGaaa                                                                                            <  1:423912/143‑1 (MQ=255)
                                                  cACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCACCGGCAGTCGCAGCTACCGACTATATGAAACGGTTCGCTg                                                                                  >  2:384883/1‑143 (MQ=255)
                                                   cccGCGCGGGCGCGTGATGGCCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGa                                                                                 <  2:724988/142‑1 (MQ=255)
                                                    ccGAGCTGGCGCGTGAGGGTCAGGATTGTGAACCCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAg                                                                                <  1:383850/143‑1 (MQ=255)
                                                        gCTGGCGCGTGATGGGCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAgg                                                                            <  1:586757/143‑1 (MQ=255)
                                                            gcgcGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCg                                                                        >  1:657604/1‑143 (MQ=255)
                                                              gcgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCCTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGCCCATATCCAACTGTTCGCTGAGCCGGTCCgga                                                                      >  2:157980/1‑141 (MQ=255)
                                                              gcgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTACGTATATCACTTAGGTGATGAACGACGCTCCGGACCCAACAGCTACCGACTATATGAAACTGTTCGAAGAGCAggacagga                                                                      >  2:620590/1‑137 (MQ=255)
                                                                 tGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGCTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGGACtt                                                                   >  2:587460/1‑143 (MQ=255)
                                                                      gTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCATAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGCCTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGtc                                                              >  2:94392/1‑142 (MQ=255)
                                                                                tgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTgacg                                                    <  1:467310/143‑1 (MQ=255)
                                                                                       tGGGACATGCGGCAGCAGCTGGAATCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcg                                             <  2:561198/143‑1 (MQ=255)
                                                                                        ggAACAGGCTGCACCCGCTGGAAACGCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcgc                                            <  1:663532/143‑1 (MQ=255)
                                                                                        gcaaCATGCCGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcgc                                            <  1:462084/141‑1 (MQ=255)
                                                                                         gAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTCACGTACCGGCTAACGACTaaccag                                           >  2:128605/1‑138 (MQ=255)
                                                                                           aCATGCTGCACCCGGTGGAAGCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCAGCTAGCGACTATAGGAAACTCTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTa                                         <  2:530672/143‑1 (MQ=255)
                                                                                                      ccGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTACGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGa                              >  1:296571/1‑143 (MQ=255)
                                                                                                      ccGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGGGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGAACTGGGGACGGa                              >  1:773914/1‑143 (MQ=255)
                                                                                                                               atgtCGCTTAGGTGACGAACGACGTCCCGGCAGTGGCATCTACCGACTATATGACACTGTTCGCTGAGCAGGTCCGTACTTGCGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGt     <  2:27787/140‑1 (MQ=255)
                                                                                                                               atatGGCTTAGGTGACGAACGACGCTCCGGCAGTGGCGTCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGt     <  2:246518/143‑1 (MQ=255)
                                                                                                                               atatCCCTTAGGTGATGACCGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGt     <  1:300586/143‑1 (MQ=255)
                                                                                                                                 atCGCTTAGGTGATGAACGACGCCCCGGCAGTGGCATCTACCGACTATATCAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTga   <  1:850110/143‑1 (MQ=255)
                                                                                                                                  tCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAGTCGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACCGCCATgag  >  2:744452/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                         
CTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAG  >  NC_000913/125252‑125524

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: