Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	693,942	A→G	intergenic (‑103/+11)	ybeZ ← / ← miaB	hypothetical protein with nucleoside triphosphate hydrolase domain/isopentenyl‑adenosine A37 tRNA methylthiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	693,942	0	A	G	100.0%	89.0 / NA	49	intergenic (‑103/+11)	ybeZ/miaB	hypothetical protein with nucleoside triphosphate hydrolase domain/isopentenyl‑adenosine A37 tRNA methylthiolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (49/0);  total (49/0)

AAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATA  >  W3110S.gb/693918‑693967                         |                          aaaaaCAAGGCCCACCGGAACTGCGGGCCTGAGAATTACGGCTGataata  >  1:79958/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGc          >  1:41795/1‑42 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:60009/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:45306/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:46831/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:48391/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:48425/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:5140/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:53018/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:53226/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:55600/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:56369/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:56612/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:59873/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:45188/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:62817/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:63396/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:64556/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:64828/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:66587/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:67332/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:82321/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:86331/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:88325/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:9557/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:26345/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:13748/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:1376/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:13974/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:14145/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:14775/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:15268/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:15820/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:17006/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:17429/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:21886/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:11904/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:28176/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:29735/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:30863/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:30873/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:32350/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:32593/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:35545/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:36104/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:36684/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataata  >  1:43518/1‑50 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataat   >  1:26996/1‑49 (MQ=255) aaaaaCAAGGCCCACCGGAACGGCGGGCCTGAGAATTACGGCTGataa    >  1:66598/1‑48 (MQ=255)                         |                          AAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATA  >  W3110S.gb/693918‑693967

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: ‑