Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,755,871	0	T	A	54.5%	0.6 / 27.5	23	G58G (GGT→GGA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (8/2);  new base A (0/12);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 1.41e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-03

CATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGG  >  NC_000913/1755724‑1756019                                                                                                                                                    |                                                                                                                                                     cATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTaaa                                                                                                                                                   <  1:256240/151‑1 (MQ=255) cATAGGACCGAAATCCGAATCATAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGATGACTATGCAGATCACGGTCACCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTaaa                                                                                                                                                   <  1:3697/151‑1 (MQ=255)                   aTCTGAAGATTTATTATCAAAAATGTTGGCCGCGGGCAGGACCGTTATGCGTCTAACCTTCTCTAGGGGTGACTAAAAAAAACACGGCCACCCCATTCAGATTCGCCCCATCGTAAGGAGAAAACCTGGAAAAACCGCCGCTATCCTGCtt                                                                                                                                 <  1:326195/151‑1 (MQ=255)                             ttGTTAGCTAAAAGGCTGGACGCGGACATAACCGTTATGCTTCTAAACTTCTCTCATGGTGACTATGAAAAACCGGGCCAGCCCATAAAGAATCTGCGAAAGGTGAGAAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAg                                                                                                                       <  1:101398/150‑1 (MQ=255)                              tGTTAGAAAAAAGGGGGGACGCGGGAAGGAACGTTAGGCGTCTAACCTTCTCTGATGGGAAAAAGGAAAAACCCGGCCACCCCATTAAGAATCGGCCAAACGTGAGAAAAAAAACTGTAAAAACCGCCGCTATCCTGCTTGATACCAAAgg                                                                                                                      <  2:365906/151‑1 (MQ=255)                                    aaaaagggggggACGCGGGAAGAAACGTTAGGCGTCTAACCTTCTCTCATGTTGAAAAGACAAAACCCGGCCAGCCCATCCAGAATCGGAGAAAGGTGATGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGa                                                                                                                <  2:369766/141‑1 (MQ=255)                                           gCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCg                                                                                                          >  1:302268/1‑150 (MQ=255)                                                            aCGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAgg                                                                                        >  1:366193/1‑151 (MQ=255)                                                            aCGTTATGCGTCTGAACTTCTCTCATGGGGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAgg                                                                                        >  2:325567/1‑151 (MQ=255)                                                                    cGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacg                                                                                >  1:75640/1‑151 (MQ=255)                                                                     tctaaacctcctctcAGGGTAATTAGGCAAAAAAGGGCAAGCGCTTTCAAAATCGGCCAAACGTGAAAAAAAAACTTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacga                                                                               <  2:362217/141‑1 (MQ=255)                                                                               ctctcATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtct                                                                        >  1:117592/1‑148 (MQ=255)                                                                                   cATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAg                                                                 >  1:276934/1‑151 (MQ=255)                                                                                    atggtcaatataaaaaaCACGCCCAGCCCAGTAAGAATCGGCCAACCGTGAGAAGAAAAGCGGAAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                                <  2:200527/136‑1 (MQ=255)                                                                                          actaggaaaaaCACGGCCAGCCCATACAGAATCGGCGCACCTTAAAGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTc                                                           <  2:230977/141‑1 (MQ=255)                                                                                                   accccGGTCAGCGCATTAAAAACCTGCGAAACGTGAAGAGAAAAACGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTAc                                                 <  1:201097/147‑1 (MQ=255)                                                                                                    acccGGTCGCCGCATTCAAATCCGGCAAAAGGTGAGGAGAAAACCTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACt                                                <  1:325483/148‑1 (MQ=255)                                                                                                    acacGGTCAGCGCATTCAAAATCGGCGCAACGGGATAAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACt                                                <  2:250839/151‑1 (MQ=255)                                                                                                     cctggacacctcATTAAGAATCGGCGCATTGTGAAGAAAAAAACGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAAGCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACtt                                               <  1:305912/140‑1 (MQ=255)                                                                                                                  ttCAGAATCTGCGCAACGTGATGAGCAAAACTGGCAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCAGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATaa                                  >  2:397238/1‑151 (MQ=255)                                                                                                                       attCTGCGAAACGGAATGAGAAACACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATc                               <  1:335997/147‑1 (MQ=255)                                                                                                                                        gaGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGaa            >  2:18809/1‑151 (MQ=255)                                                                                                                                                  ggTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCgg  >  2:367206/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGG  >  NC_000913/1755724‑1756019

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