Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,498,089	0	T	C	61.9%	11.8 / 19.6	21	V204V (GTA→GTG)	frc	formyl‑CoA transferase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (8/0);  new base C (0/13);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 4.91e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCC  >  W3110S.gb/2498032‑2498098                                                          |          aTAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCATTACgg       >  1:1548541/1‑62 (MQ=255) aTAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACgg       >  1:822508/1‑62 (MQ=255) aTAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACgg       >  1:3647916/1‑62 (MQ=255) aTAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACgg       >  1:2987891/1‑62 (MQ=255) aTAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACgg       >  1:2940899/1‑62 (MQ=255) aTAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACgg       >  1:1910909/1‑62 (MQ=255) aTAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACgg       >  1:1892164/1‑62 (MQ=255) aTAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACgg       >  1:174519/1‑62 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAATGGTTCAACACGGCATcc  <  1:550746/62‑1 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAAAGGTTCAACACGGCATcc  <  1:2087379/62‑1 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAAAGGTTCAACACGGCATcc  <  1:2149072/62‑1 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAAAGGTTCAACACGGCATcc  <  1:2394711/62‑1 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAAAGGTTCAACACGGCATcc  <  1:2973861/62‑1 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAAAGGTTCAACACGGCATcc  <  1:3502039/62‑1 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAAAGGTTCAACACGGCATcc  <  1:3541786/62‑1 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAAAGGTTCAACACGGCATcc  <  1:3926108/62‑1 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAAAGGTTCAACACGGCATcc  <  1:4474845/62‑1 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAAAGGTTCAACACGGCATcc  <  1:697372/62‑1 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAAAGGTTCAACACGGCATcc  <  1:735794/62‑1 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAAAGGTTCAACACGGCATcc  <  1:737847/62‑1 (MQ=255)      ccAACTTATCGAGACGCTGCTGGTCGCGTAATTTCACGCGGCAAAGGTTCAACACGGCATcc  <  1:1443454/62‑1 (MQ=255)                                                          |          ATAACCCAATTTATCGAGACGCTGCTGGTCACGTAATTTCACGCGGCAAAGGTTCAATACGGCATCC  >  W3110S.gb/2498032‑2498098

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