Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,733,280	0	T	A	53.3%	10.1 / 19.9	15	N360K (AAT→AAA)	xylG	D‑xylose ABC transporter dual domain ATPase
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base A (8/0);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 1.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03

CGGACAATGGGAAGGAAAAATTTATATTGATGGCAAACAGGTAGATATTCGTAACTGTCAGCAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTT‑ACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCAAAAATGTATTCTGGAATCAATCCAGCAACTCAAAGTTAAAACGTCGTCCCCCGACCTTGCTA  >  NC_000913/3733136‑3733407                                                                                                                                                 |                                                                                                                                 cGGACAATGGGAAGGAAAAATTTATATTGATGGCAAACAGGTAGATATTCGTAACTGTCAGCAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGGGAAAAAAtttc                                                                                                                             >  1:155039/1‑148 (MQ=255)          ggaaggaaAAATTTATATTGATGGCAAACAGGTAGATATTCGTAACTGTCAGCAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAAAATTAACCCcgcc                                                                                                                      >  2:543600/1‑148 (MQ=255)          ggaaggaaAAATTTATATTGATGGCAAACAGGTAGATATTCGTAACTGTCAGCAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGGAAAAATATTACCCCcgccgc                                                                                                                    >  2:322547/1‑150 (MQ=255)            aaggaaAAATTTATATTGATGGCAAACAGGTAGATATTCGTAACTGTCAGCAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTTGTAAAAAAATTACCCTccccgc                                                                                                                    >  2:84158/1‑148 (MQ=255)                 aaaaTTTATATTGATGGCAAACAGGTAGATATTCGTAACTGTCAGCAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATa                                                                                                            <  2:276157/151‑1 (MQ=255)                 aaaaTTTATATTGATGGCAAACAGGTAGATATTCGTAACTGTCAGCAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAAAATTAACCCCCCCCCACTc                                                                                                                >  1:490876/1‑147 (MQ=255)                   aaTTTATATTGATGGCAAACAGGTAGATATTCGTAACTGTCAGCAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATaaa                                                                                                          <  2:518800/151‑1 (MQ=255)                         tatTGATGGCAAACAGGTAGATATTCGTAACTGTCAGCAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGGAAAAAAATTACCCTCGCCGCCCTCAATAAATTT‑Acc                                                                                                   >  2:38569/1‑151 (MQ=255)                                   aaaCAGGTAGATATTCGTAACTGTCAGCAGGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTT‑ACCGGTGGCAtt                                                                                          <  1:520598/150‑1 (MQ=255)                                    aaCAGGTAGATATTCGTAACTGTCAGCAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTTGTTAAAAAAATAACCACGGCGCCCCCAATAAATTTGACAGGTGTCAt                                                                                           >  1:336461/1‑149 (MQ=255)                                                           cAGCAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAAAATTACCCCCGCCCCACCCAATAAATTT‑AACCGTGGGCTTAACCCGCTTTATGACGCGGc                                                                      >  1:373013/1‑146 (MQ=255)                                                                    aTCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTT‑ACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCAAAAATGTa                                                        <  2:392236/151‑1 (MQ=255)                                                                                          tCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGGAAAAAAATTAACCCCCCCGCCCTCAATAAATTT‑TCCGGTGGGATTAGCCAGCTTGATGACGCGGCCGAGCAAAAATGTATTCTGGAATCAATCCCGCaa                                    >  1:197658/1‑149 (MQ=255)                                                                                                aaGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTT‑ACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCAAAAATGTATTCTGGAATCAATCCAGCAACTCAAAGt                            <  1:401009/151‑1 (MQ=255)                                                                                                                          ccGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTT‑ACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCAAAAATGTATTCTGGAATCAATCCAGCAACTCAAAGTTAAAACGTCGTCCCCCGACCTTGCTa  <  1:171555/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                 CGGACAATGGGAAGGAAAAATTTATATTGATGGCAAACAGGTAGATATTCGTAACTGTCAGCAAGCCATCGCCCAGGGGATTGCGATGGTCCCCGAAGACAGAAAGCGCGACGGCATCGTTCCGGTAATGGCGGTTGGTAAAAATATTACCCTCGCCGCACTCAATAAATTT‑ACCGGTGGCATTAGCCAGCTTGATGACGCGGCAGAGCAAAAATGTATTCTGGAATCAATCCAGCAACTCAAAGTTAAAACGTCGTCCCCCGACCTTGCTA  >  NC_000913/3733136‑3733407

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