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

CCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGA  >  NC_000913/125251‑125523                                                                                                                                         |                                                                                                                                         gcTGGCGAAAGCTTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGGGATGGTCAGGATTGTGAACGCTGGAATATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtga                                                                                                                                    <  2:71208/142‑1 (MQ=255)   tGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGAATATCGCTTAGGtgtgg                                                                                                                                  >  2:334733/1‑140 (MQ=255)    ggCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgatga                                                                                                                                 <  1:399469/143‑1 (MQ=255)     gCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGaa                                                                                                                                <  1:530358/143‑1 (MQ=255)     gCGAAAGATTACGGCGTAGGTTCTGACGTTTAGAGCGTGACCTCCTTCACCGAGCTGGCGCGCGATGGTCAGGATTGGGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGaa                                                                                                                                <  1:286067/143‑1 (MQ=255)      cGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGACTGGGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGGGCGTTCCGTATATCGCCTAGGagatgaca                                                                                                                               >  2:450366/1‑141 (MQ=255)               aCGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGGACGCTGGAACATGCTGCACCCGCTGGAACCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGCCGCTCCg                                                                                                                      >  2:197944/1‑143 (MQ=255)                        gtgcTGCCGTTTATAGCGTGACCGCCTTCTCCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGACACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCa                                                                                                             <  2:34085/140‑1 (MQ=255)                         ttCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCAt                                                                                                            <  1:96803/143‑1 (MQ=255)                                     taGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGGCAGTATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACtata                                                                                                <  1:447242/143‑1 (MQ=255)                                                ccctCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTc                                                                                     <  1:467601/140‑1 (MQ=255)                                                       cGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGc                                                                              >  1:139884/1‑143 (MQ=255)                                                           cgggCGCGTGAGGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGt                                                                          <  2:68849/141‑1 (MQ=255)                                                              gcgcGTGATGGTCAGGAGTGTGACCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCCGACATCGCTTAGGTGATGAACGGCGCTCCGGCCGTGGCATCCACCGACCATATGAAACCGTCCGCTCAGcgcgcccg                                                                       >  2:441817/1‑136 (MQ=255)                                                                       ggTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGCTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTCACGt                                                              >  2:280056/1‑143 (MQ=255)                                                                                                 gctgcCCCCGCGGGAACCTCCGCGCGGTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                    <  2:124161/143‑1 (MQ=255)                                                                                                 gctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                    >  1:383654/1‑143 (MQ=255)                                                                                                        ccGCTGGAACCTCCGCGCGGTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGGACCGGCTGACGACTACCGCGTACTGGGTACTGa                             <  1:526527/143‑1 (MQ=255)                                                                                                        ccGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCGGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGACCGTACTTACGTACCGGCGGACGACTACCGCGTACTGGGCACTGa                             >  2:6383/1‑143 (MQ=255)                                                                                                           cTGGAACGTCTGCGCGTTCCGGATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATgg                          <  2:310651/143‑1 (MQ=255)                                                                                                                 ccTCCGCTCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTACATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCgg                    <  1:144943/142‑1 (MQ=255)                                                                                                                            tCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGAGTACTGGGGACTAATGGCTTCGCACGTTCCGCCa         >  2:17365/1‑143 (MQ=255)                                                                                                                                   atCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCCGGTCCGTACTTACGTACCGGCTGACGCCTACCGCGTACTGGGTCCTCATGGGCTCGGTCGTTTCGACAGACGTga  >  2:95081/1‑143 (MQ=255)                                                                                                                                         |                                                                                                                                         CCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGA  >  NC_000913/125251‑125523

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