Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,269,124	T→G	S407A (TCT→GCT)	dnaB →	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,269,124	0	T	G	82.3%	47.5 / 1.5	17	S407A (TCT→GCT)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base G (0/14);  total (3/14)
Fisher's exact test for biased strand distribution p-value = 1.47e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCA  >  W3110S.gb/4269067‑4269151                                                          |                            tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt                         <  1:2059971/59‑1 (MQ=255) tcgctGGAACAACGCCCCAACAAACGTCCCGTGATGTCCGACCTGCGTGAATCCGGCGCTAt                         <  1:967102/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:1156263/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)                        aCGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATc   >  1:1277194/1‑61 (MQ=255)                        aCGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCa  >  1:2516492/1‑62 (MQ=255)                        aCGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCa  >  1:1681276/1‑62 (MQ=255)                                                          |                            TCTCTGGAACAACGTGCCGACAAACGCCCGGTCAACTCCGACCTGCGTGAATCTGGCTCTATCGAGCAGGATGCGGACTTGATCA  >  W3110S.gb/4269067‑4269151

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