Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,642,892	(T)5→6	coding (48/1155 nt)	yfgB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,642,887	1	.	T	100.0%	73.9 / NA	28	I18N (ATC→AAC)	yfgB	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (0/28);  total (0/28)

GGTTGAGATCCAGCAGGTTGA‑TTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAA  >  W3110S.gb/2642867‑2642927                      |                                         ggTTGAGATCCAGCTGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:1133900/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTTGTGACGTTTTCAGGTGTGACTaa  <  1:1348943/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:14847/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:997672/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:995706/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:994010/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:978195/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:857530/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:83909/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:796624/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:553157/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:454070/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:269117/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:266698/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:187808/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:176998/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:1023005/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:1353607/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:1343981/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:1325523/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:1242470/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:118137/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:1162191/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:1157916/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:1121993/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:1117922/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:1115303/62‑1 (MQ=255) ggTTGAGATCCAGCAGGTTGATTTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTaa  <  1:1094074/62‑1 (MQ=255)                      |                                         GGTTGAGATCCAGCAGGTTGA‑TTTTTCCATCTTTCGTGGTGACGTTTTCAGGTGTGACTAA  >  W3110S.gb/2642867‑2642927

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