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%	106.6 / NA	40	Q791* (CAG→TAG)	aceE	pyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/20);  total (20/20)

TCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAG  >  NC_000913/125250‑125524                                                                                                                                          |                                                                                                                                          tCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtg                                                                                                                                      >  2:298763/1‑143 (MQ=255) tCCTGGCGAAAGATTACGGCGTAGGTTCTGACGGTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGGTTAGGtg                                                                                                                                      >  1:506577/1‑143 (MQ=255)  ccTGGCGAAAGATGACGGGGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtga                                                                                                                                     <  2:593689/143‑1 (MQ=255)      gCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGaa                                                                                                                                 <  1:298763/143‑1 (MQ=255)        gAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCGAGCAACTCCCCGCGTTCGCTACAGCGTTTAAGTGATGAaca                                                                                                                               >  2:109901/1‑142 (MQ=255)         aaaGATTCCGGCGCAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacga                                                                                                                              <  1:362452/143‑1 (MQ=255)         aaaGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGCacga                                                                                                                              >  2:143170/1‑143 (MQ=255)                tagtcGTAGGTTATGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGCACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCg                                                                                                                       <  1:260232/139‑1 (MQ=255)                  ggCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTACACCTGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGc                                                                                                                     >  1:527333/1‑143 (MQ=255)                     gTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTAGCGCGTGATGGTCAGGATTGTGCTCGCTGGAACATGCGGCACCCGCGTGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGt                                                                                                                  <  1:185194/143‑1 (MQ=255)                      ggggTTCTGTTGTTTATGGCTTGACCTCCTCCATCCATTCGGCGTGTGATGGTCAGGATTGTGTAGGCTCGTGCACTCCCAACAAGAATGAAGTTCCGTGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTg                                                                                                                 <  2:283399/141‑1 (MQ=255)                      gaGGTTCTGACGTTTATGGCGTGACCTCCTGCACCGGGCTGGCGCGTGGTGGTCAGGATTGGGAACGCTGGTACATGCTGCTCCCGCGGGAACCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTg                                                                                                                 <  2:633758/142‑1 (MQ=255)                       aGGTTCTGACGATTGTAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTgg                                                                                                                <  1:708097/143‑1 (MQ=255)                          ttCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCAt                                                                                                             >  1:490105/1‑143 (MQ=255)                          ttCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGAGCGTTCCGTATATCGCTTAGGTGATGAACGCCGCTCCGGCAGTGGCAt                                                                                                             >  2:121162/1‑143 (MQ=255)                              gACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATGTCGGATAGCTGATGAACGACGCTCCGGCCGTGGCATCTAc                                                                                                         >  2:133176/1‑143 (MQ=255)                                      ttgCGTGATCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtata                                                                                                 <  1:608642/141‑1 (MQ=255)                                             aCCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACt                                                                                          >  1:654933/1‑143 (MQ=255)                                                 cctTCACCGAGCGGGCGCGTGATGGTCAGGAGTGGGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                                                                      <  1:103189/143‑1 (MQ=255)                                                 cctTCACCGACCTGGCGCGTGCTGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                                                                      <  1:46782/143‑1 (MQ=255)                                                  ctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGATGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGCTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGGTCg                                                                                     >  2:511728/1‑143 (MQ=255)                                                       ccGAGCTTGCGCGTGATGGTCATGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTAGATGAAACTGTTCGCTGAg                                                                                <  1:136255/143‑1 (MQ=255)                                                               gcgcGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCg                                                                        >  1:236951/1‑143 (MQ=255)                                                                 gcgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTa                                                                      <  2:200949/143‑1 (MQ=255)                                                                            aGGATTGTGAACGCGGGAACATGCTGCACCCGCTGGAAACGCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCg                                                           <  2:368257/143‑1 (MQ=255)                                                                              gATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAGACTGTTCGCTGAGCAGGTCCGTGCTTACGTACCGGc                                                         >  1:64504/1‑143 (MQ=255)                                                                                          tGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACcg                                             >  2:464204/1‑143 (MQ=255)                                                                                               cATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTAc                                        <  1:762582/143‑1 (MQ=255)                                                                                                aTGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCCGCGCAGGTCCGTACTTACTTACCGGCTGACCACTACCGCGTCCt                                       >  2:430957/1‑143 (MQ=255)                                                                                                  gctgcACCCGCTGGAAACCCCGCGCGTTCCGTATATCGCTTAGGGGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                     <  2:334348/143‑1 (MQ=255)                                                                                                    tgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGt                                   >  1:609405/1‑143 (MQ=255)                                                                                                    tgcACCCGCTGGAAACTCCGCGCGTGCCGGATATCGCATAGGTGATGAACGACGCGCCGGCCGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGCCCGTACGTACGTACCGGCTGACGAATACCGCGAACtacgc                                   >  2:1819/1‑139 (MQ=255)                                                                                                        cccGCTGGAGACGCCGCGCGTTCCGTATGTCGCTTAGGTGACGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTg                               <  1:303826/143‑1 (MQ=255)                                                                                                                 aaCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCg                      >  2:95315/1‑143 (MQ=255)                                                                                                                  aCTCCGCCCGTTCCGTAGATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCgg                     <  1:133176/143‑1 (MQ=255)                                                                                                                           gTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGa            <  1:650707/143‑1 (MQ=255)                                                                                                                               cGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGGTTCGGTCGTTCCGACAGc        >  1:41897/1‑143 (MQ=255)                                                                                                                                gTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACTGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGcc       <  2:236951/143‑1 (MQ=255)                                                                                                                                gTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGcc       >  1:14474/1‑143 (MQ=255)                                                                                                                                     tCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGACCAGGTCCGTACTTACGTACCGGCTGCCGACTACCGCGTCCTGGGTACTGATGGCTTCGGTCGTTCCACCAGCCgcgcg  >  1:223635/1‑139 (MQ=255)                                                                                                                                          |                                                                                                                                          TCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAG  >  NC_000913/125250‑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: ‑