Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	729,608	A→G	Y18C (TAT→TGT)	ybfA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	729,608	0	A	G	100.0%	109.8 / NA	57	Y18C (TAT→TGT)	ybfA	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (0/57);  total (0/57)

AAGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCG  >  W3110S.gb/729552‑729612                                                         |     aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:68778/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:10612/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:52279/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:56752/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:57021/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:57356/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:58728/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:59140/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:64685/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:66083/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:6648/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:66687/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:67290/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:67383/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:51665/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:69495/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:70363/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:73134/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:73676/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:75598/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:77659/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:78276/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:79145/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:81492/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:82311/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:84478/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:8766/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:9882/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:23506/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:11384/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:11840/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:13178/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:14271/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:151/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:17027/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:17161/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:17815/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:18685/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:19405/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:20624/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:20976/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:2290/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:48472/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:24494/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:26243/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:27813/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:28571/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:28620/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:34572/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:34859/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:35985/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:40275/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:41001/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:46758/61‑1 (MQ=255) aaGTATGGAACTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:47611/61‑1 (MQ=255)          aCTCTACAAAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:68960/52‑1 (MQ=255)                          ccTGGATGGCTTATCTTTTTACGCCGTACTTGTGCg  <  1:67429/36‑1 (MQ=37)                                                         |     AAGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCG  >  W3110S.gb/729552‑729612

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