Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_012967	649,391	T→A	100%	I471F (ATC→TTC)	mrdA ←	peptidoglycan DD‑transpeptidase MrdA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012967	649,391	0	T	A	100.0%	37.4 / NA	14	I471F (ATC→TTC)	mrdA	peptidoglycan DD‑transpeptidase MrdA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (7/7);  total (7/7)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGG > NC_012967/649362‑649420 | aaCCTTCACGATACCGTCATTAATCAGGAACATCAg < 1:235356/36‑1 (MQ=37) aaCCTTCACGATACCGTCACTAATCAGGAacatcaa > 1:451453/1‑35 (MQ=25) cTTCACGATACCGTCATTAATCAGGAACATCAGTGc < 2:543819/36‑1 (MQ=37) ttCACGATACCGTCATTAATCAGGAACATCAGTGcc < 1:570901/36‑1 (MQ=37) ttCACGATACCGTCATTAATCAGGAACATCAGTGcc > 1:900870/1‑36 (MQ=37) cACGTTACCGTCATTAATCAGGAACATCAGTGCCtt < 1:919479/36‑1 (MQ=25) cACGATACCGTCATTAATCAGGAACATCAGTGCCtt < 1:138512/36‑1 (MQ=37) cACGATACCGTCATTAATCAGGAACATCAGCGCCtt > 2:824486/1‑36 (MQ=21) cGATATCGTCATTAATCAGGAACATCAGTGCCTTAc < 2:727904/36‑1 (MQ=25) cGATACCGTCATTAGTCAGGAACATCAGTGCCTTAc > 1:465473/1‑36 (MQ=37) cGATACCGTCATTAATCCGGAACATCAGTGCCttcc > 1:647247/1‑34 (MQ=25) gATACCGTCATTAATCAGGAACATCAGTGCCTTACt < 1:269849/36‑1 (MQ=37) ccGTCATTAATCAGGAACATCAGTGCCTTACTCATc > 2:806664/1‑36 (MQ=38) aaTCAGGAACATCAGTGCCTTACTCATCTCGATTgg > 2:737523/1‑36 (MQ=37) tCAGGAACATCAGTGCCTTACTCATCTGGATTgggg > 1:279810/1‑36 (MQ=38) | AACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGG > NC_012967/649362‑649420

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