Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,582,247	A→T	R86R (CGT→CGA)	ydeM ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,582,247	0	A	T	100.0%	120.0 / NA	67	R86R (CGT→CGA)	ydeM	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (0/67);  total (0/67)

ATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAA  >  W3110S.gb/1582234‑1582294              |                                                tttcATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:67688/58‑1 (MQ=255) aTGCATTAAAAATTCTTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:41884/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCGTGCATAGCGTTGTTGATTGTGAATAACTTTACGGaa  <  1:14199/61‑1 (MQ=38) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTTTTTATAGTGAATAACTTTACGGaa  <  1:43960/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:60820/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:70032/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:68347/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:66296/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:64662/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:64046/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:61182/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:70611/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:59333/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:5730/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:54910/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:53550/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:53202/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:53141/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:49250/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:71472/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:73951/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:75859/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:78446/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:81448/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:82063/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:82117/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:82519/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:82679/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:83717/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:85383/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:8620/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:9048/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:9732/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:9916/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:36367/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:11450/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:11465/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:11930/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:13225/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:13692/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:14940/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:15792/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:17512/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:26432/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:27630/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:28120/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:28164/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:31103/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:34530/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:36195/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:42957/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:10924/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:46621/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:44407/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:44338/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:43666/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:43474/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:52664/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:42836/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:42314/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:41121/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:39853/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:38020/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:36683/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTTAGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:8909/61‑1 (MQ=255) aTGCATTAAAAATTCGTTTATGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:63842/61‑1 (MQ=255)  tGCATTAAAAATTCGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGaa  <  1:45253/60‑1 (MQ=255)              |                                                ATGCATTAAAAATACGTTTTTGGCCTGCATAGCGTTGTTGATAGTGAATAACTTTACGGAA  >  W3110S.gb/1582234‑1582294

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