Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,269,085	G→A	D394N (GAC→AAC)	dnaB →	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,269,085	0	G	A	77.8%	29.1 / 5.6	18	D394N (GAC→AAC)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base A (0/14);  total (0/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTAT  >  W3110S.gb/4269045‑4269128                                         |                                            ggCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTcc                         <  1:1232897/61‑1 (MQ=255) ggCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTcc                         <  1:1798107/61‑1 (MQ=255) ggCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTcc                         <  1:218905/61‑1 (MQ=255) ggCGCTGGCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTcc                         <  1:2718774/61‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:967102/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:1156263/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:940666/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:29036/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:2702613/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:2657334/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:2481426/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:2072536/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:2059971/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:2057201/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:1857772/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:1537191/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:1452226/59‑1 (MQ=255)                       tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt  <  1:1339546/59‑1 (MQ=255)                                         |                                            GGCGCTGTCCCAGTTGAACCGTTCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTAT  >  W3110S.gb/4269045‑4269128

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