Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,098,457	G→A	T159T (ACG→ACA)	hisF →	imidazole glycerol phosphate synthase, catalytic subunit with HisF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,098,457	0	G	A	84.9%	99.4 / 8.0	53	T159T (ACG→ACA)	hisF	imidazole glycerol phosphate synthase, catalytic subunit with HisF
Reads supporting (aligned to +/- strand):  ref base G (0/8);  new base A (45/0);  total (45/8)
Fisher's exact test for biased strand distribution p-value = 1.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAGGAAGTGCAAAAAC  >  W3110S.gb/2098414‑2098485                                            |                             ataCCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAGGaa            <  1:824234/62‑1 (MQ=255) ataCCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAGGaa            <  1:206993/62‑1 (MQ=255) ataCCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAGGaa            <  1:337891/62‑1 (MQ=255) ataCCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAGGaa            <  1:919339/62‑1 (MQ=255) ataCCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAGGaa            <  1:1322438/62‑1 (MQ=255) ataCCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAGGaa            <  1:1306950/62‑1 (MQ=255) ataCCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAGGaa            <  1:1274049/62‑1 (MQ=255) ataCCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAGGaa            <  1:578709/62‑1 (MQ=255)                                     gggAACCACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:628015/1‑36 (MQ=25)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCaaaaa   >  1:1384381/1‑35 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:820186/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:821014/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:911376/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:837603/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:861634/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:869373/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:815233/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:877453/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:778058/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:889360/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:908030/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:575003/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:913186/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:918087/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:927495/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:937896/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:944834/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:949279/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:959101/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:974276/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:978943/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1426238/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1032481/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1162102/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1184477/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1267659/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1308949/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1313548/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1316274/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1360417/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1386193/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1424578/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:676765/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1437486/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1463449/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:265851/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:341175/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:350085/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:360036/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:442352/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:524004/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:1004410/1‑36 (MQ=37)                                     gggAAACACTCGACTGGGTGCAGGAAGTGCAAAAAc  >  1:643736/1‑36 (MQ=37)                                            |                             ATACCGGCGATGAAAGCCGCACCCGCGTCACTCAATGGGAAACGCTCGACTGGGTACAGGAAGTGCAAAAAC  >  W3110S.gb/2098414‑2098485

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