Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	125,335	A→G	E773E (GAA→GAG)	aceE →	pyruvate dehydrogenase E1 component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	125,335	0	A	G	100.0%	49.4 / NA	15	E773E (GAA→GAG)	aceE	pyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (6/9);  total (6/9)

CTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACC  >  NC_000913/125209‑125466                                                                                                                               |                                                                                                                                    cTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAACAtgc                                                                                                                        >  1:40127/1‑140 (MQ=255)      gTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCAc                                                                                                                   <  2:371841/140‑1 (MQ=255)                   tCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACTc                                                                                                      >  2:348436/1‑140 (MQ=255)                      cGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGCTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACTCcgc                                                                                                   <  1:293585/140‑1 (MQ=255)                        tCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACTCcgcgc                                                                                                 <  1:309577/140‑1 (MQ=255)                                       agaTCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCt                                                                                  <  1:348436/140‑1 (MQ=255)                                          tCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACGCCGCGCGTTCCGTATATCGCGCAg                                                                               <  2:1987/140‑1 (MQ=255)                                              ggCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGtga                                                                           >  2:40364/1‑140 (MQ=255)                                                    aGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAacg                                                                     >  2:148556/1‑140 (MQ=255)                                                               tAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCa                                                          <  1:40364/140‑1 (MQ=255)                                                                             tataGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGa                                            >  2:269851/1‑140 (MQ=255)                                                                                         tcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACt                                <  1:269851/140‑1 (MQ=255)                                                                                                        gcgcGTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGt                 <  1:148556/140‑1 (MQ=255)                                                                                                           cgTGATGGTCAGGATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCg              <  2:191206/140‑1 (MQ=255)                                                                                                                       gATTGTGAGCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAcc  >  2:342248/1‑140 (MQ=255)                                                                                                                               |                                                                                                                                    CTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACC  >  NC_000913/125209‑125466

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑