Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,664,820	1	.	G	20.4%	114.0 / 38.6	44	coding (602/1263 nt)	rffD	UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (35/0);  new base G (9/0);  total (44/0)
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
Rejected as consensus: Frequency below/above cutoff threshold.

TCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTC  >  W3110S.gb/3664763‑3664827                                                           |        tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACt   >  1:861239/1‑64 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:584668/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1637019/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1664718/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1862189/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:272158/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:290409/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:303244/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:355482/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:505444/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:53868/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1116084/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:78429/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:808931/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:846772/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:852243/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:861904/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:908762/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:985303/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1394113/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1037455/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1137639/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1147968/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1160355/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1168943/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1182816/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1232048/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1247563/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1307944/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1381503/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1606504/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1434054/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1483649/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1502472/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTc  >  1:1569249/1‑65 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGAAGAGcacacg   >  1:49438/1‑64 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGAAGAGcacacg   >  1:404130/1‑64 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGAAGAGcacacg   >  1:672000/1‑64 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGAAGAGcacacg   >  1:708953/1‑64 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGAAGAGcacacg   >  1:1440598/1‑64 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGAAGAGcacacg   >  1:1482853/1‑64 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGAAGAGcacacg   >  1:286773/1‑64 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGAAGAGcacacg   >  1:1497594/1‑64 (MQ=255) tCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGAAGAGcacacg   >  1:181911/1‑64 (MQ=255)                                                           |        TCGCGGAAGCTGTTTTCGGTGAGCTTACACATTTCCGCCGTCCGCGAGTTAGTGACGA‑CACACTC  >  W3110S.gb/3664763‑3664827

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