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% 97.9 / ‑6.7 37Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/20);  minor base G (0/1);  total (16/21)
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.

AGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTC  >  NC_000913/125247‑125507
                                                                                                                                            |                                                                                                                        
agaTCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAg                                                                                                                        <  1:103258/143‑1 (MQ=255)
  aTCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGCACAACGCTcttgt                                                                                                                      <  2:261961/143‑5 (MQ=255)
     cgggCGAAATTTTTCGTCGTAGGTTCTGATGTTTATAGCGTGACCTCCTTCACCGGGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCGGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgat                                                                                                                   <  1:346924/141‑1 (MQ=255)
     cTGGCGAAAGATTACGGCGTAGGTTCTGGCGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgat                                                                                                                   <  2:10975/143‑1 (MQ=255)
         cGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTAAGCGCGATACGTATATAGCTTAGGAGATGAac                                                                                                               >  1:201647/1‑143 (MQ=255)
           aaaGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGCGATGAacga                                                                                                             >  2:287459/1‑143 (MQ=255)
           aaaGACCACGCCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacga                                                                                                             <  1:292416/143‑1 (MQ=255)
                ttgcGGCGTAGGTTCTGACGTTTATAGCGTGAACTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTc                                                                                                        <  1:15336/140‑1 (MQ=255)
                ttACGGCGGAGGTTCTGACGTTTATAGCATGAGCTCATTCACAGAGGTGGCGCGGGAGGGTGAGTAGTGGGAACTCTGGCACAGGCTGCACGCGGTGGAAACTCCGCGCGTACCGTATATCGCATAGGTGGTGAACGACGCta                                                                                                        >  2:316330/1‑142 (MQ=255)
                   cGGCGTGGGTGCGGACGTGTATAGCGTGACCTCCTTCACCGAGCTGGCGCGCGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACGCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCgg                                                                                                     <  2:122271/143‑1 (MQ=255)
                   cGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCgg                                                                                                     <  2:505236/143‑1 (MQ=255)
                          gttgcTGACGTGTATAGCGTGGCCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCATATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGc                                                                                              <  2:194588/139‑1 (MQ=255)
                               cggcGTTTATGGTGCGGCCTCCCACACCGAGCTGGCGCGTGATGGTCAGGAGTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTCTACCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTa                                                                                         <  2:130609/140‑1 (MQ=255)
                                gACGTTTATAGCGTGACCGCCTGCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAc                                                                                        <  2:581003/143‑1 (MQ=255)
                                      tatgGCGTGACGTCTCTCGGCGGGCTGGGCCGGGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACGCGCTGGAAACTCCGCGCGTTCCGAATAACAAATCGAGAATAAACGACGCTCCGGCAGTGGCATCTACCGACta                                                                                  <  2:287957/143‑1 (MQ=255)
                                      tataGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACta                                                                                  <  1:218368/143‑1 (MQ=255)
                                         aGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTatat                                                                               <  1:513237/143‑1 (MQ=255)
                                          cccgcaCCTCCTTCCCCGACCTGGCGCGTCACGGTCAGGATTGTGCACGTTGGAACCTGCTGCACCCGCCGTAAACGCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATg                                                                              <  1:262940/138‑1 (MQ=255)
                                             tGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGAGTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGAGATAAACGACGCTCCGGCAGTGGCCTCTACCGACTATATGaaa                                                                           >  2:190698/1‑143 (MQ=255)
                                                cctcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTGCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGCAACTg                                                                        >  1:270888/1‑143 (MQ=255)
                                                  ctctTCACCTAGAGGACGCGTGATGGAAAGAATTTTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGAATAGCGCTGAGGTGATGAACCACGCTCCGGCAGTGGCATCTACCGACGATATGAAACtttt                                                                      <  1:206350/141‑1 (MQ=255)
                                                    ctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTCCACCCGCTGGAATCTCCGCGCGTTCCGTATATCGCTTAGGTGATCAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg                                                                    <  1:541518/143‑1 (MQ=255)
                                                          cGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAGACTGGGCGCTGAGc                                                              >  2:65215/1‑143 (MQ=255)
                                                                ggCGCGTGATGGTCAGGATTGTGAACGCTCGAACATGCTTCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTcc                                                        <  2:21312/143‑1 (MQ=255)
                                                                              aGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCg                                          >  2:293900/1‑143 (MQ=255)
                                                                              aGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATACGAAACCGTTCGCTGAGCAGGTCCGTACTCAGGTACCg                                          >  2:439508/1‑143 (MQ=255)
                                                                                gATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGCCGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGc                                        >  1:454655/1‑143 (MQ=255)
                                                                                            tGGAACTTCCTGCCCCCGTGGGCACCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGTTCGGTCAGTGACATCTACCGACTATAGTAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcg                            <  2:350802/143‑1 (MQ=255)
                                                                                                 cATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCGGGTCCGTACTTACGTACCGGCTGACGACTAGcgtctac                       >  2:341861/1‑138 (MQ=255)
                                                                                                  gtgcggcACCGGCTGGAACCTGCGGCGGGTCCGTATATCGCTTAGGTGAGGGACGACGCTCCGGCAGTGGCATCTACCGACTAGATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACt                      <  2:484157/142‑1 (MQ=255)
                                                                                                       gcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTa                 >  1:81837/1‑143 (MQ=255)
                                                                                                        ctcccGCCGGATACTCCGCGCGGTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAATCTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTAc                <  2:469040/141‑1 (MQ=255)
                                                                                                        cACCCGCTGGAAACTCCGCGCGCTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCGACCGACGATATGAAACTGTCCGCGGAGCAGGTCCGTCATTACGAACAGGAGGCCGACAACCGCGTACGGGGCAc                >  2:592946/1‑143 (MQ=255)
                                                                                                          cccGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGACATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACGGGGTACTg              >  2:373542/1‑143 (MQ=255)
                                                                                                              cTGGAAACTCCGCGCGTGCCGTATATCGCTTAGGTGATGAGCGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATgg          <  1:188162/143‑1 (MQ=255)
                                                                                                                   aaCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCg     >  2:33325/1‑143 (MQ=255)
                                                                                                                   aaCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTAACGACTAGATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTGCTGGGCAATGAAGACTTCg     >  2:547014/1‑143 (MQ=255)
                                                                                                                      tCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGCCGACTACCGCGTACTGGGGAGTGATGGCTTCggcc  >  1:292582/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                        
AGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTC  >  NC_000913/125247‑125507

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