Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,631,788	G→A	Q352Q (CAG→CAA)	slt →	soluble lytic murein transglycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,631,788	0	G	A	100.0%	37.1 / NA	14	Q352Q (CAG→CAA)	slt	soluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/6);  total (8/6)

TGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACAGCGTGGTTTCTACCCGATGGTTGCAGCACAACGCATCGGCGAAGAGTATGAGCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTG  >  NC_000913/4631591‑4631950                                                                                                                                                                                                      |                                                                                                                                                                   tGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATTGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAAGCGGATTTATTGCTGGAAcgn                                                                                                                                               <  1:49920/219‑2 (MQ=255)               cAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCGCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATAGCGTTACGGGCAAGCGGATTTATTGCTGGAACGCGGACGTGAAGCTTAAGCAAAAGAGAGTTTGCATCAAc                                                                                                          <  1:361353/242‑1 (MQ=255)                       tCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAAGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACAGCGTGGt                                                                                         >  1:26793/1‑251 (MQ=255)                        cACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTCTTGAACGACGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGCGAAAGATGAATGGCGTTACTGGCAAGCGGATTTATTGCTGGAACGAGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACAGCGTGGtt                                                                                        >  1:302102/1‑251 (MQ=255)                                       gAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAAGCGGATTTATTGCTGGAAcgcgn                                                                                                                                             <  1:51985/183‑2 (MQ=255)                                                                       nnnnnnnnnngCNNNNNGAANGCCGTGTACGAAGGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAAGCGGATTTATTGCTGGAACGCGGACg                                                                                                                                          <  2:267542/144‑1 (MQ=255)                                                                                       tGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGTCTTGCGCGTCTGCCTATGGAAGCGAAAGAGAAATATGAATGGCGTTACTGGCAAGCGTATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACAGCGTg                                                                                           <  1:301762/185‑1 (MQ=255)                                                                                          aCGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAAGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACAGCGTGGTTTCTACCCGATGGTTGCAGCACAACGCATCGGCGAAGAGTATGAGCTGAAGATTGATAAAGcgccgc                      >  1:158439/1‑251 (MQ=255)                                                                                                          ncgcTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGGTGACGCGTCTGCCTATGGCAGCGAACGAGAAAGATGAATGGCGTTACTGGAAAGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACAGCGTGGTTTCTACCCGATGGTTGCAGCACAACGCATCGGCGAAGAGTATGAGCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAgcgc      >  1:83602/2‑251 (MQ=255)                                                                                                          gcgcTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAAGCGGCTTTATTGCTGGAACGCGGACGTGCAGCTGAAGCAAAAGAGATTTTGCATCCACTCATGCAACAGCGTGGTTTCTACCCGATGGTTGCAGCACAACGCATCGGCGAAGAGTATGAGCTGAAGATTGATAAAGCGCCGCAg                    >  1:213028/1‑237 (MQ=255)                                                                                                               tGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGCAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAAGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACAGCGTGGTTTCTACCCGATGGTTGCAGCACAACGCATCGGCGAAGAGTATGAGCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTg  >  1:413386/1‑250 (MQ=255)                                                                                                                                   gCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAAGCGGCTTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACAGCGTGGTTTCTACCCGATGGTTGCAGCACAACGCATCGGCGAagag                                                  >  1:367037/1‑182 (MQ=255)                                                                                                                                    ccTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGAGGAATGGCGTTACTGGCAAGCTGATTTATTGCTGGAACGCGTACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATg                                                                                                     <  1:246338/130‑1 (MQ=255)                                                                                                                                                                              gaaagaTGAATGGCGTTACTGGCAAGCGGATTTATTGCTGGAACGCGGACGTGACGCTGAAGCAATAGAGATTTTGCATCAACTCCTGCAACAGCGTGGTTTCTACCCGATGGTTGCAGCACAACGCATCGGCGAAGAGTATGAGc                                           >  1:178532/1‑146 (MQ=255)                                                                                                                                                                                                      |                                                                                                                                                                   TGGCGTTTGATGGGCAACGATGTCACCGACGAGCAGGCGAAATGGCGCGATGACGCCATTATGCGCTCGCAATCTACTTCGCTTATTGAACGCCGTGTACGAATGGCGCTTGGCACCGGCGATCGTCGCGGCCTGAATACCTGGCTGGCGCGTCTGCCTATGGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACAGCGTGGTTTCTACCCGATGGTTGCAGCACAACGCATCGGCGAAGAGTATGAGCTGAAGATTGATAAAGCGCCGCAGAATGTTGACAGCGCCCTG  >  NC_000913/4631591‑4631950

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