Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,543,104	1	.	T	63.2%	30.5 / 22.9	19	coding (317/375 nt)	yfeK	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/7);  new base T (0/12);  total (0/19)
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.

GGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTTACATGCGTTAATTGCGCAGACGGATAAAA  >  W3110S.gb/2543052‑2543139                                                      |                                    gggAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTTACATg                          <  1:589605/64‑1 (MQ=255) gggAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTTACATg                          <  1:1461423/64‑1 (MQ=255) gggAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTTACATg                          <  1:1692325/64‑1 (MQ=255) gggAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTTACATg                          <  1:1899953/64‑1 (MQ=255) gggAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTTACATg                          <  1:2436405/64‑1 (MQ=255) gggAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTTACATg                          <  1:2638309/64‑1 (MQ=255)  ggAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTTACATg                          <  1:2662957/63‑1 (MQ=255)                         cccGGTAAGAGCGATGAGAACGCACAGCCTTTTTTTACATGCGTTAATTGCGCAGACGGATaaaa  <  1:855932/65‑1 (MQ=255)                         cccGGTAAGAGCGATGAGAACGCACAGCCTTTTTTTACATGCGTTAATTGCGCAGACGGATaaaa  <  1:826030/65‑1 (MQ=255)                         cccGGTAAGAGCGATGAGAACGCACAGCCTTTTTTTACATGCGTTAATTGCGCAGACGGATaaaa  <  1:584690/65‑1 (MQ=255)                         cccGGTAAGAGCGATGAGAACGCACAGCCTTTTTTTACATGCGTTAATTGCGCAGACGGATaaaa  <  1:346508/65‑1 (MQ=255)                         cccGGTAAGAGCGATGAGAACGCACAGCCTTTTTTTACATGCGTTAATTGCGCAGACGGATaaaa  <  1:2761275/65‑1 (MQ=255)                         cccGGTAAGAGCGATGAGAACGCACAGCCTTTTTTTACATGCGTTAATTGCGCAGACGGATaaaa  <  1:2754008/65‑1 (MQ=255)                         cccGGTAAGAGCGATGAGAACGCACAGCCTTTTTTTACATGCGTTAATTGCGCAGACGGATaaaa  <  1:2298585/65‑1 (MQ=255)                         cccGGTAAGAGCGATGAGAACGCACAGCCTTTTTTTACATGCGTTAATTGCGCAGACGGATaaaa  <  1:1717612/65‑1 (MQ=255)                         cccGGTAAGAGCGATGAGAACGCACAGCCTTTTTTTACATGCGTTAATTGCGCAGACGGATaaaa  <  1:1648170/65‑1 (MQ=255)                          ccGGTAAGAGCGATGAGAACGCACAGCCTTTTTTTACATGCGTTAATTGCGCAGACGGATaaaa  <  1:1441641/64‑1 (MQ=255)                          ccGGTAAGAGCGATGAGAACGCACAGCCTTTTTTTACATGCGTTAATTGCGCAGACGGATaaaa  <  1:1307064/64‑1 (MQ=255)                                agagCGATGAGAACGCACAGCCTTTTTTTACATGCGTTAATTGCGCAGACGTATaaaa  <  1:1225848/58‑1 (MQ=255)                                                      |                                    GGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCC‑TTTTTTACATGCGTTAATTGCGCAGACGGATAAAA  >  W3110S.gb/2543052‑2543139

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