Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,992,354	0	G	T	26.4%	132.1 / 35.6	53	G21C (GGT→TGT)	ygeI	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/39);  new base T (0/14);  total (0/53)
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.

AATTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAACA  >  W3110S.gb/2992318‑2992382                                     |                             aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:18417/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:304458/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:742164/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:1790843/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:195251/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:1506575/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:881185/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:1366001/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:413508/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:131836/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:1316182/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:1311418/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATAAaca  <  1:1260252/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAGaca  <  1:1744173/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:688182/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:493513/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:646615/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:601214/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:523516/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:731724/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:490071/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:474500/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:69025/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:726262/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:729664/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:987284/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:732629/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:769402/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:812191/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:820968/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:903654/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:1222475/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:443863/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:1349870/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:1471885/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:1473046/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:1540712/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:1590214/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:1652786/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:1735999/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:1744933/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:1776907/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:1792676/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:1823559/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:245052/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:292330/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:310352/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:330039/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:33901/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:464321/65‑1 (MQ=255) aatTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGCGATGAAGTCGTATCATTAGATAAaca  <  1:1231476/65‑1 (MQ=255)       tctcAAAGTTCAAATATAGCGAATGCAACGTGTGATGAAGTCGTATCATTAGATaaaga  <  1:764138/59‑3 (MQ=255)         tcAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAaca  <  1:193514/57‑1 (MQ=255)                                     |                             AATTTATCTCAAAGTTCAAATATAGCGAATGCAACGGGTGATGAAGTCGTATCATTAGATAAACA  >  W3110S.gb/2992318‑2992382

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