Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,999,954	0	G	T	70.7%	49.3 / 23.4	41	A318A (GCG→GCT)	xdhA	xanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/12);  new base T (29/0);  total (29/12)
Fisher's exact test for biased strand distribution p-value = 1.27e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTA  >  W3110S.gb/2999897‑3000001                                                          |                                                ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggg                                              >  1:1113001/1‑61 (MQ=255) ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggg                                              >  1:1068468/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCCCGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:779483/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCCCTCCATTGCTTCTGCTGGGGGGAATAAAg                                    >  1:306496/1‑58 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:880264/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:8456/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:795911/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:944970/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:68371/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:667695/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:660044/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:548618/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:546869/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:495293/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:484883/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:41191/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:405941/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:344938/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:990373/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:284141/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:1117538/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:1153464/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:1164205/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:1222738/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:1276843/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:1308156/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:1431997/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:1457452/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:210184/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:251313/1‑61 (MQ=255)              aaCACCGGCGCTTATGCATCTCACGGGCACTCCATTGCTTCTGCTGGGGGGAATAAAGTCg                                 >  1:279252/1‑61 (MQ=255)                                             ccATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa  <  1:534157/61‑1 (MQ=255)                                             ccATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa  <  1:202924/61‑1 (MQ=255)                                             ccATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa  <  1:627577/61‑1 (MQ=255)                                             ccATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa  <  1:1340483/61‑1 (MQ=255)                                             ccATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa  <  1:681009/61‑1 (MQ=255)                                             ccATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa  <  1:468665/61‑1 (MQ=255)                                             ccATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa  <  1:1237108/61‑1 (MQ=255)                                             ccATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa  <  1:787833/61‑1 (MQ=255)                                             ccATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa  <  1:1235724/61‑1 (MQ=255)                                             ccATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa  <  1:362179/61‑1 (MQ=255)                                             ccATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa  <  1:30814/61‑1 (MQ=255)                                              cATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa  <  1:777829/60‑1 (MQ=255)                                                          |                                                GGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTA  >  W3110S.gb/2999897‑3000001

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