Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,609,129	(GAT)3→4	coding (234/624 nt)	bglJ →	DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,609,120	1	.	G	100.0%	14.7 / NA	6	A75? (GCG→GCN)	bglJ	DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (0/6);  total (0/6)
*	W3110S.gb	4,609,120	2	.	A	100.0%	14.0 / NA	6	A75A (GCG→GCA)	bglJ	DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (0/6);  total (0/6)
*	W3110S.gb	4,609,120	3	.	T	100.0%	13.7 / NA	6	A75A (GCG→GCT)	bglJ	DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (0/6);  total (0/6)

CGGCGTTTAGTTATTGCG‑‑‑GATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGC  >  W3110S.gb/4609103‑4609160                   |||                                         cGGCGTTTAGTTATTGCGGATGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGc  <  1:1247990/61‑1 (MQ=255) cGGCGTTTAGTTATTGCGGATGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGc  <  1:1371021/61‑1 (MQ=255) cGGCGTTTAGTTATTGCGGATGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGc  <  1:2026832/61‑1 (MQ=255) cGGCGTTTAGTTATTGCGGATGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGc  <  1:2667616/61‑1 (MQ=255) cGGCGTTTAGTTATTGCGGATGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGc  <  1:289824/61‑1 (MQ=255) cGGCGTTTAGTTATTGCGGATGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGc  <  1:300781/61‑1 (MQ=255)                   |||                                         CGGCGTTTAGTTATTGCG‑‑‑GATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGC  >  W3110S.gb/4609103‑4609160

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