Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,227,909	(A)5→6	coding (366/2352 nt)	ygjK →	predicted glycosyl hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,227,904	1	.	A	87.5%	40.7 / 1.5	16	coding (361/2352 nt)	ygjK	predicted glycosyl hydrolase
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (14/0);  total (14/2)
Fisher's exact test for biased strand distribution p-value = 8.33e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

GAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGC‑AAAAACTGACAGCAAAAGATGTGC  >  W3110S.gb/3227855‑3227928                                                   |                         gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:1417128/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:18106/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:1810925/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:1912642/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:201272/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:2031697/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:2048138/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:231076/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:250021/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:2657945/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:3094643/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:327837/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:380249/1‑61 (MQ=255) gAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGCAAAAAACTGAc                >  1:506170/1‑61 (MQ=255)               cGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGC‑AAAAACTGACAGCAAAAGATGTGc  <  1:153483/60‑1 (MQ=255)               cGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGC‑AAAAACTGACAGCAAAAGATGTGc  <  1:2275159/60‑1 (MQ=255)                                                   |                         GAAAGTCGACTTTACGCTGGAGGCATACAGTATTCCCGGTGCGCTGGTGC‑AAAAACTGACAGCAAAAGATGTGC  >  W3110S.gb/3227855‑3227928

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