Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	631,633	G→A	W580* (TGG→TGA)	cstA →	pyruvate transporter CstA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	631,633	0	G	A	100.0%	29.0 / NA	12	W580* (TGG→TGA)	cstA	pyruvate transporter CstA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/5);  total (7/5)

CTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTC  >  NC_000913/631451‑631879                                                                                                                                                                                       |                                                                                                                                                                                                                                                       ctgGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACAATCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGAGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGcagcacggaagagcacac                                                                                                                                                                                    >  1:116613/1‑234 (MQ=255)  tgGCANNCGNNNTGNNCGTACTGGCGTGGGGCTACTTCCNCCATCAGGGCGTGGTCGATCCGCTGGGCGNCNTNAACANNCTGNGGCCGCTGTTNGGTATTGCCAAACAGCTGCTGGCCGGGATGGGGCTGATGCTCTGGGCCNGNGTGTTGTTCAAGATGAAACGTCAACGTTCCGNCTGAGTGGGGCCGGAACCAACGGCNCGgcnnctngn                                                                                                                                                                                                                        >  2:130902/1‑211 (MQ=255)    gCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAAc                                                                                                                                                                                                                                       <  1:37238/197‑1 (MQ=255)                        ggCGTGGGGCTACTTCCTCCATCAGGGCGTGTTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGGTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGGACCAACGGTCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTc                                                                                                                                                              <  1:130902/250‑1 (MQ=255)                               ggCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCCACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGcc                                                                                                                                                                                  >  1:58845/1‑223 (MQ=255)                                        ctccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGAATGCAATTGCCAACCAGATGCTGGCAGAGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGACGCTGGTACCAACGGCATGGCTGCTGCTTTGTACCCTGACCCCAGGCTGGCAGAAAGCGTTTAGCCCGGATGAGCAAGTAGGCTTCCTGGGCATAGGtaata                                                                                                                                              >  1:28308/1‑250 (MQ=255)                                                           aTCCGCTGGGCGGCATTAACACTCTCNGGCCGCNNTTNCGAATTGCCACCCAGATGCNGGNNNGNNNGNCGCCGNCANNCTATCCCGTNGTGTTGTTCAAGATGAAACCTCAACGTTACGCCAGAATGGCGCCGATACCCAGNGCCCCGCTGNNGATTTTTACCCTGACCGCAGGNTGGNNANNANCggntan                                                                                                                                                                                    >  2:37078/1‑188 (MQ=255)                                                                                                                     ggCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAATCGTCAACGTTACGCCGGAGTGGCGCTGGTACCAACGGCCGGGCTGCTGATTTGTACCCGGACCGCAGGCTGTCAGAAAGCGTGTAGCCCGGATGGGAAAGTCGGCTGCCTGGCCAGTGCGAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTc                                                                                                             <  1:88274/206‑1 (MQ=255)                                                                                                                                                    gtgtTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGAGGCGAAAGTCGGCTTCCTGTCCATTGCTAATAAGTGCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAGCAACCGTCTGGATGCCGGGTTAAc                                          <  1:51320/242‑1 (MQ=255)                                                                                                                                                              aGATGNAANGNCANNGNTACGNCTGAGTGGCGCTGGTACCAACGGCCTg                                                                                                                                                                                                                                 >  2:78828/1‑49 (MQ=255)                                                                                                                                                                                   ccTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTGCAACAACCGTCTGGAGGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTtc  <  1:70514/251‑1 (MQ=255)                                                                                                                                                                                   ccTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGGATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGttt  >  1:11917/1‑250 (MQ=255)                                                                                                                                                                                       |                                                                                                                                                                                                                                                       CTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTC  >  NC_000913/631451‑631879

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