Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,970,117	G→A	R152R (CGG→CGA)	mutH →	DNA mismatch repair protein MutH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,970,117	0	G	A	100.0%	16.2 / NA	8	R152R (CGG→CGA)	mutH	DNA mismatch repair protein MutH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/3);  total (5/3)

GTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATCCCGCTGGCGCAGCGTCGCGTTGGATCACCGTTGCTGTGGAGCCCGAATGAAGAGGAAGACCGGCAGCTACGCGAAGACTGGGAAGAATTAATGGATATGATTGTTCTCGGTCAGGTTGAGCGGATTACCGCTCGTCACGGGGAGTATTTACAGATACGACCGAAAGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTG  >  NC_000913/2969921‑2970342                                                                                                                                                                                                     |                                                                                                                                                                                                                                  gtcgtcCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATCCCGCTGGCGCAGCGTCGCGTTGGATCACCGTTGCTGTGGAGCCCGAATGAAGAGGAAGACCGACAGCTACGCGAAGACTGGGAAGAATTAATGGATATGATTGTTCTCGGTCAGGtt                                                                                                                                                                             >  1:166983/1‑251 (MQ=255)                                   cccGTTAACGGGCAATCGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATCCCGCTGGCGCAGCGTCGCGTTGGATCACCGTTGCTGTGGAGCCCGAATGAAGAGGAAGACCGACAGCTACGCGAAGACTGGGAAGAATTAATGGATATGATTGTTCTCGGTCAGGTTGAGCGGATTACCGCTCGTCACGGGGAGTATTTAc                                                                                                                                           >  1:46103/1‑251 (MQ=255)                                                                       cAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATCCCGCTGGCGCAGCGTCGCGTTGGATCACCGTTGCTGTGGAGCCCGAATGAAGAGGAAGACCGACAGCTACGCGAAg                                                                                                                                                                                                                      >  1:185208/1‑140 (MQ=255)                                                                                                     aCTGTGGATACCGGTTGAAGGCGAGAGCAGCATCCCGCTGGCGCTGCGTCGCGTTGGATCACCGTTGCGGTGGAGCCAGAATGAAGAGGAAGAGCGACAGCTACGCGAAGACTGGGAAGAAGTAAGGGATATGATTGTTCTCGGTCAGGTTGAGCGGATTACCGCTCGTAACGGGGAGTATTTACAGATACGTCCGAAATCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGa                                                                         <  1:93698/251‑1 (MQ=255)                                                                                                                                          cTGGCGCAGCGTCGCGTTGGATCACCGTTGCTGTGGAGCCCCAATGAAACGGAAGACCGACTGCTCAGCGAAGACTGGGAAGAATTAATGGATATGATTGTTCTCGGTCAGGTTGAGCGGATTACAGCTCGTCACGGGGAGAATTTACAGATACGACCGAAAGCAGCGAATGCGAAAGAGCTTACCGAAGCCATTGGTGCCCGGGGc                                                                                >  1:99549/1‑207 (MQ=255)                                                                                                                                                        ngTTGGATCACCGTTGCTGTGGAGCCCGAATGAAGAGGCCGCCCGACAGCTACGAGACGACTGGGCAGAATTAATGGATATGATTGTTCTCGGTCAGGTTGAGCGGATTACCGCTCGTCACGGGGAGTATTTACAGATACGACCGAAAGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCac                      >  1:39784/2‑251 (MQ=255)                                                                                                                                                                nnccNNNNNTGTNGAGCCCGAATGAAGAGGAAGACCGACAGCTACGCGAAGACTGGGAAGAATTAATGGATATGATTGTTCTCGGTCAGGTTGAGCGGATTACCGCTCGTCACGGGGAGTATTTACAGATACGACCGAAAGCAGCg                                                                                                                       <  2:134235/144‑1 (MQ=255)                                                                                                                                                                            ggAGCCCGAATGAAGAGGAAGACCGACAGCTACGCGAAGACTGGGAAGAATTAATGGATATGATTGTTCTCGGTCAGGTTGAGCGGATTACCGCTCGTCACGGGGAGTATTTACAGATACGACCGAAAGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTg  <  1:196068/251‑1 (MQ=255)                                                                                                                                                                                                     |                                                                                                                                                                                                                                  GTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATCCCGCTGGCGCAGCGTCGCGTTGGATCACCGTTGCTGTGGAGCCCGAATGAAGAGGAAGACCGGCAGCTACGCGAAGACTGGGAAGAATTAATGGATATGATTGTTCTCGGTCAGGTTGAGCGGATTACCGCTCGTCACGGGGAGTATTTACAGATACGACCGAAAGCAGCGAATGCGAAAGCGCTTACCGAAGCCATTGGTGCCCGGGGCGAACGGATTCTGACGCTGCCACGCGGCTTTTATTTGAAGAAGAATTTCACCAGTGCACTACTGGCCCGTCATTTTCTG  >  NC_000913/2969921‑2970342

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