Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,406,261	G→A	V95M (GTG→ATG)	smrA →	DNA endonuclease SmrA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,406,261	0	G	A	100.0%	18.3 / NA	9	V95M (GTG→ATG)	smrA	DNA endonuclease SmrA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (4/5);  total (4/5)

CCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGC  >  NC_000913/1406063‑1406493                                                                                                                                                                                                       |                                                                                                                                                                                                                                         ccAACGCGTAACCAACGTACGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCCGCACCGGATTTCTCGACATCATCCAACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGGTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATACGAAACAGGCGAGCCTGAATCTTTTGCGCCCGCCGATGGACGAGTGCCGCAAAATGGTGTTCAGTTTTATTccac                                                                                                                                                                                                   >  1:217300/1‑235 (MQ=255)                                  gACACGCTGCAGCTTGATAATTTCCTCACCACCGTATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGGTGCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCGGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGATGGAAGAGTGCCGCAAAAGGGTGTTCAGGTTTATGCAACAAGCTGTGGCGGATGGttg                                                                                                                                                                                <  1:89822/224‑2 (MQ=255)                                      ngCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTACGCCGCTGGAGTTTCGGCGGGAAGGGCTGCAACATGGGGTGCTGGATATGCTGCGCAGTGGTACATATCCGCAACAGGAGCGCCTGAATCTTTTGCGCCAGCCGATGGAAGAGTGCCGCACAATGGTGTTCAGTTTTATTCAACAAGCCCTGGGGGATGGTTTGCGTACCGTGCTGATTATTCATGGTAAAggg                                                                                                                                                 >  1:33498/2‑251 (MQ=255)                                                          cTCACCACCGGATTTCTCGACATCATCCCACTAAGTCCGAAGCTGGAGTTTCGGCGGGAAGGGGTGCAACATGGGGTGCTGGCTAAGCTGCGCAGTGGTAAATATCCGAAACAGGAGAGAATGAATCATTTGCGCCAGCCGATGGACGAGTGCCGCAAAATGGAGTTCAGTTTTATTCAACAAGCCCTAGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGAGGGAAGATAAATCGCATGc                                                                                                                             >  1:55175/1‑251 (MQ=255)                                                                                                                            cAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGATGGAAGAGTGCCGCAAAATGGTGTTCATTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATAGTGCAc                                                           <  1:182981/251‑1 (MQ=255)                                                                                                                                                                 atCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGATGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATAGTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATAGTGCACCGCGCTACCGCATCATGGCGGCAGCgggg                              <  1:123479/243‑1 (MQ=255)                                                                                                                                                                                     aaTCTTTTGCGCCAGCCGATGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTGAAGGGCGGGATGATAAATCGCATGCCAATATGGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTGCAGGCATATTGCACCGCGCGACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGgc  <  1:327404/251‑1 (MQ=255)                                                                                                                                                                                                ccAGCCGATGGAAGAGTGCCGCAAAATGGTNTTCAGTTTTANTCANNAAGCNCNGNCTGATGGTTTGCGTAACGTGCTGATTATTNATGGTAAAGGGCGGGATGATAAANNGCNNNCCNNTATTg                                                                                                                     <  2:244509/125‑1 (MQ=255)                                                                                                                                                                                                ccAGCCGATGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTg                                                                                                                     >  1:244509/1‑125 (MQ=255)                                                                                                                                                                                                       |                                                                                                                                                                                                                                         CCAACGCGTAACCAACGTGCGCCGCAGCGTATCGACACGCTGCAGCTTGATAATTTCCTCACCACCGGATTTCTCGACATCATCCCACTAAGTCAGCCGCTGGAGTTTCGGCGGGAAGGGTTGCAACATGGGGTGCTGGATAAGCTGCGCAGTGGTAAATATCCGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATAAATCGCATGCCAATATTGTCCGCAGCTATGTGGCGCGCTGGCTGACCGAATTTGATGATGTTCAGGCATATTGCACCGCGCTACCGCATCATGGCGGCAGCGGGGCGTGTTACGTCGCACTACGTAAAACGGC  >  NC_000913/1406063‑1406493

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