Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,874,542	A→G	I148I (ATT→ATC)	cysD ←	sulfate adenylyltransferase, subunit 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,874,542	0	A	G	100.0%	147.4 / NA	53	I148I (ATT→ATC)	cysD	sulfate adenylyltransferase, subunit 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (0/53);  total (0/53)

AGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGG  >  W3110S.gb/2874536‑2874597       |                                                        agagTAGATTCGCTCTTTAGCGCGGGTTTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:581031/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:683470/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1017839/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:408989/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:453682/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:488694/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:492138/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:497164/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:581926/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:586914/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:62142/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:633252/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:66142/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:31818/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:726657/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:812309/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:830317/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:838325/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:841065/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:879660/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:901029/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:949937/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:961118/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:974162/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:977543/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:993888/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1434037/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1091113/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:110944/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1111280/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1147236/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1171930/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1237667/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:132350/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1343525/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1347771/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1348026/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1408126/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:279043/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1459660/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:182261/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:204342/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:204823/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:222225/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:268100/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:271916/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:276482/62‑1 (MQ=255) agagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGAGCACCACCGAAggcgg  <  1:1211872/62‑1 (MQ=255) agagTAGAATCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:464674/62‑1 (MQ=255)  gagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCTCACCACCGAAggcgg  <  1:646608/61‑1 (MQ=255)  gagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:160831/61‑1 (MQ=255)  gagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:1333397/61‑1 (MQ=255)  gagTAGATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAggcgg  <  1:218063/61‑1 (MQ=255)       |                                                        AGAGTAAATTCGCTCTTTAGCGCGGGATTTCTCTTCGTCACGGCGCGCACCACCGAAGGCGG  >  W3110S.gb/2874536‑2874597

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