Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,346,205	T→C	T879A (ACG→GCG)	yfaL ←	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,346,205	0	T	C	75.6%	74.7 / 22.9	45	T879A (ACG→GCG)	yfaL	adhesin
Reads supporting (aligned to +/- strand):  ref base T (11/0);  new base C (0/34);  total (11/34)
Fisher's exact test for biased strand distribution p-value = 9.85e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAA  >  W3110S.gb/2346146‑2346220                                                            |                ccccATATTGATGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:895320/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:377150/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:213414/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:2157126/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:229810/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:2499652/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:2642046/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:2651673/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:2802220/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:2842226/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:2906127/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:1121176/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:428656/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:457318/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:553102/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:71587/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:755950/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:898580/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:921292/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:1938550/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:119081/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:1205748/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:1308442/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:1327784/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:152801/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:154006/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:1569719/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:1609239/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:1610251/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:1621362/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:1655572/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:177232/61‑1 (MQ=255) ccccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:1854717/61‑1 (MQ=255)   ccATATTGACGTAGCCGCTTCCTGTCAGACTGAACTGCGCATTGTTTTGAAACTGCGCgg               <  1:2874910/60‑1 (MQ=255)               aGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGaa  >  1:2782764/1‑61 (MQ=255)               aGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGaa  >  1:268856/1‑61 (MQ=255)               aGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGaa  >  1:2678731/1‑61 (MQ=255)               aGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGaa  >  1:1609143/1‑61 (MQ=255)               aGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGaa  >  1:540563/1‑61 (MQ=255)               aGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGaa  >  1:2458180/1‑61 (MQ=255)               aGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGaa  >  1:1870678/1‑61 (MQ=255)               aGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGaa  >  1:1885936/1‑61 (MQ=255)               aGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGaa  >  1:1950770/1‑61 (MQ=255)               aGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGaa  >  1:1058699/1‑61 (MQ=255)               aGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGAGTGGGATCAGCTGCGaa  >  1:2177802/1‑61 (MQ=255)                                                            |                TCCCATATTGACGTAGCCGCTGCCTGCCAGACTGAACTGCGCATTGTTTTGAAACTGCGTGGGATCAGCTGCGAA  >  W3110S.gb/2346146‑2346220

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