Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,173,325	1	.	T	72.0%	58.4 / 21.9	25	intergenic (+67/‑102)	chiA/bfd	periplasmic endochitinase/bacterioferritin‑associated ferredoxin
Reads supporting (aligned to +/- strand):  ref base . (7/0);  new base T (18/0);  total (25/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 = 8.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATACCTTCCAGTATCGCCGCATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTA‑TTTTTATTCATATTTCTATTCTTT  >  W3110S.gb/4173270‑4173349                                                         |                         aTACCTTCCAGTATCGCCGCATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTA‑TTTTTAtt                  >  1:1191001/1‑64 (MQ=255) aTACCTTCCAGTATCGCCGCATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTA‑TTTTTAtt                  >  1:91667/1‑64 (MQ=255) aTACCTTCCAGTATCGCCGCATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTA‑TTTTTAtt                  >  1:1626602/1‑64 (MQ=255) aTACCTTCCAGTATCGCCGCATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTA‑TTTTTAtt                  >  1:374408/1‑64 (MQ=255) aTACCTTCCAGTATCGCCGCATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTA‑TTTTTAtt                  >  1:1973580/1‑64 (MQ=255) aTACCTTCCAGTATCGCCGCATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTA‑TTTTTAtt                  >  1:219401/1‑64 (MQ=255) aTACCTTCCAGTATCGCCGCATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTA‑TTTTTAtt                  >  1:2349876/1‑64 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:2038427/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:997938/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:921373/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:773274/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:492727/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:351646/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:286590/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:2378437/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:2215288/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:1008041/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:2007886/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:1955404/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:1885493/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:1884546/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:1860870/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:1497784/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:1255532/1‑65 (MQ=255)                 ccgcATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTATTTTTTATTCATATTTCTATTCttt  >  1:1062707/1‑65 (MQ=255)                                                         |                         ATACCTTCCAGTATCGCCGCATTTTTTGTCAATGGGCTACTTTTGGTAGCCCATTA‑TTTTTATTCATATTTCTATTCTTT  >  W3110S.gb/4173270‑4173349

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