Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,210,315	Δ1 bp	coding (1287/2280 nt)	yehM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,210,315	0	A	.	93.3%	51.1 / ‑1.4	15	coding (1287/2280 nt)	yehM	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (0/14);  total (1/14)
Fisher's exact test for biased strand distribution p-value = 6.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CGGCTGGCAATACTGGAGATCCCTGGGATTGTACGCCAGCAGGGAAGTACACTGACACTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCT  >  W3110S.gb/2210259‑2210359                                                         |                                             cGGCTGGCAATACTGGAGATCCCTGGGATTGTACGCCAGCAGGGAAGTACACTGACACTTGc                                         >  1:1251122/1‑62 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:1093878/62‑1 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:1129100/62‑1 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:1363982/62‑1 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:1423041/62‑1 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:223911/62‑1 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:24065/62‑1 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:262826/62‑1 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:368779/62‑1 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:385715/62‑1 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:430580/62‑1 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:822773/62‑1 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:850155/62‑1 (MQ=255)                                       gcagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:863253/62‑1 (MQ=255)                                        cagGGAAGTACACTGAC‑CTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCt  <  1:538635/61‑1 (MQ=255)                                                         |                                             CGGCTGGCAATACTGGAGATCCCTGGGATTGTACGCCAGCAGGGAAGTACACTGACACTTGCAGGCAACGGTGAAGAACGCTGGAAATTAACCCGCCCGCT  >  W3110S.gb/2210259‑2210359

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑