Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,471,772	A→G	F424S (TTC→TCC)	btuB ←	vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,471,772	0	A	G	100.0%	102.9 / NA	54	F424S (TTC→TCC)	btuB	vitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (0/54);  total (0/54)

ACTGTTTGCTTTTCTCCGGGTCCAGATTCGGATTTCCGTAGAAGCCATACAGTTGCCCCAG  >  W3110S.gb/3471731‑3471791                                          |                    attGTTTGCTTTTCTCAGGATTCAGATTCGTATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:39143/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGGCCAAg  <  1:88318/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:62245/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:57339/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:6274/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:64988/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:65107/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:66574/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:67886/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:68315/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:69250/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:70374/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:70782/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:72031/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:72610/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:74266/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:76626/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:7767/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:78530/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:78625/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:79420/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:80079/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:80972/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:81695/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:81834/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:82628/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:86454/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:88284/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:32648/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:12017/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:15697/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:16534/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:22153/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:25402/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:26662/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:26957/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:27457/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:31277/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:31824/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:10606/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:36569/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:38113/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:38197/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:39290/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:41559/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:47129/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:48384/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:53289/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:54010/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:54020/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:55888/59‑1 (MQ=255) attGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTCCCCAAg  <  1:18309/59‑1 (MQ=255)  ttGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:46374/59‑1 (MQ=255)  ttGTTTGCTTTTCTCAGGATTCAGATTCGGATTACCGTAGGAGCCGTACAGTTGCCCAAg  <  1:13857/59‑1 (MQ=255)                                          |                    ACTGTTTGCTTTTCTCCGGGTCCAGATTCGGATTTCCGTAGAAGCCATACAGTTGCCCCAG  >  W3110S.gb/3471731‑3471791

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