Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	125,387	C→T	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_000913	125,387	0	C	T	100.0%	82.0 / NA	31	Q791* (CAG→TAG)	aceE	pyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/19);  total (12/19)

GAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTT  >  NC_000913/125246‑125510                                                                                                                                              |                                                                                                                            ttgtttgtGGCGAAAGTTTTCGTCGTAGGTTCT‑ATGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGTATTGTGAACGCTGGAACATGCTGCACCCGCTGTAAACTCCGCGCGTTCCGTATATCGCTTAg                                                                                                                           <  2:116526/136‑1 (MQ=255)   gaTCCTGTCGAAGGTGTGTGGCGTCGGTGCTGGCGTCTATAGCGTGGCCGCCTTCACCGAGCTGGCGCGTCATCGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAgg                                                                                                                          <  2:428604/143‑1 (MQ=255)      ccTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtta                                                                                                                       <  1:252631/143‑3 (MQ=255)                     cGGCGTAGGTGCGGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCgg                                                                                                        <  2:358317/143‑1 (MQ=255)                                  gACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACCTGCTGCACCCGCCGGCAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGACCGACCCTCCCGCACTGGCacccac                                                                                           >  2:190380/1‑138 (MQ=255)                                             cGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACAGGCTGCACCCGCGGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGa                                                                                <  1:522613/143‑1 (MQ=255)                                                   ctcctTCACCGAGCTGGTGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCAGGCAGTGGCATCTACCGACTATATGGAACTGt                                                                          >  1:423654/1‑143 (MQ=255)                                                        tCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCt                                                                     <  1:127829/143‑1 (MQ=255)                                                            cGAGCTGGCGGGTGCTGCCCATGTTGTTGAACTGTGGAACATGGTGTATCCGCTGGAAGCTCCGCGCGTTCCGTATTTCGCTTAGGTGATGAACGACGCTTCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGc                                                                 <  2:333156/143‑1 (MQ=255)                                                            cGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCCGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGc                                                                 >  2:303342/1‑143 (MQ=255)                                                                     gcgGGATGGTACGGAGTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGGGGCATCTACCGACTATAGGAAACTGTTCGCTGAGCAGGTCCGTa                                                        <  2:315531/143‑1 (MQ=255)                                                                                     tgtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCGCCGGCAGGGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTga                                        <  2:201501/143‑1 (MQ=255)                                                                                      gtgAACATTGCAACATGCTGCACCCGCTTTTTACTGCGCGCTTTCCGTATCACTCTTAGATGCTGCACGTTGATATGGCATTGGCATCTACCCACTATAAAATTTAGTTAGGTGGTCATGTCCTTACTTACGTACCGGCTgac                                       <  2:270939/143‑1 (MQ=255)                                                                                      caagctgCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATAGCGCTTAGGTGATGAGCGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTgac                                       <  1:532412/137‑1 (MQ=255)                                                                                         aaCGCTGGAACATGGTGACCCTTCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTgacgac                                    <  2:119691/143‑1 (MQ=255)                                                                                         aaCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTGTATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTTCCGACTATATGAAACTGTTCGCTGCGCAGGTCCGCACTTACGCACCGGCTgacgcc                                    >  2:685699/1‑141 (MQ=255)                                                                                            gCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTAc                                 >  1:699845/1‑143 (MQ=255)                                                                                             cTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGGCTGcc                                >  1:467711/1‑143 (MQ=255)                                                                                             cTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGACAGGGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACCTACGTACCGACcgacgagagcg                                >  2:532412/1‑138 (MQ=255)                                                                                                gAACACGCTGCACCCGCTGGAATCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcgcg                             <  2:477346/143‑1 (MQ=255)                                                                                                      gctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                       <  1:461699/143‑1 (MQ=255)                                                                                                       ctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTggg                      >  1:10934/1‑143 (MQ=255)                                                                                                           aCCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATAGTAAACTGTCCGCTCCCCACCGCCGCACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACt                  >  1:321430/1‑143 (MQ=255)                                                                                                             gcGCTGGAAACGCCGCGCGTTCCGTGTATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGa                <  2:713862/142‑1 (MQ=255)                                                                                                                cTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATgg             >  1:387979/1‑143 (MQ=255)                                                                                                                cTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGAGGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTGCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATgg             <  1:306057/143‑1 (MQ=255)                                                                                                                   gAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGCTGGcgt          >  2:697508/1‑141 (MQ=255)                                                                                                                    aagcTCCGCGTGTTCCGTATATTTCTTAGTTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATTAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGCCTATTGCGTACTGGGTACTGATGGCTTc         <  2:466081/140‑1 (MQ=255)                                                                                                                       cTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGt      <  2:380539/143‑1 (MQ=255)                                                                                                                        tCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGCCTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTCCCGCGGTCTCGGTACTGATCGCTTCGGTc     >  1:38679/1‑143 (MQ=255)                                                                                                                           gcgcgGTCCGTATATCGCTTAGGTGATGGACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGtt  <  1:190380/143‑1 (MQ=255)                                                                                                                                              |                                                                                                                            GAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTT  >  NC_000913/125246‑125510

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: ‑