Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,755,871	0	T	A	42.0%	65.4 / 32.9	50	G58G (GGT→GGA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (28/1);  new base A (1/20);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 9.06e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGG  >  NC_000913/1755721‑1756013                                                                                                                                                       |                                                                                                                                               cACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGt                                                                                                                                                >  2:887487/1‑151 (MQ=255)   ccATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTaa                                                                                                                                              <  1:520923/151‑1 (MQ=255)       cGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAcc                                                                                                                                          >  1:1364794/1‑151 (MQ=255)          aCCGAAAACAGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccgcc                                                                                                                                       >  2:800011/1‑151 (MQ=255)             gAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCt                                                                                                                                    >  1:860693/1‑151 (MQ=255)              aaaaCCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCt                                                                                                                                    >  2:356886/1‑150 (MQ=255)                    gAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGc                                                                                                                             >  1:481731/1‑151 (MQ=255)                           aaGAGTTGTTAGAAAAAATGCTGGACGCGGGCATGAACTTTATGCGTCTGACCTTCTCTCAGGGTGACTATGAAAACCGCGGTCAGCGAATTCAAAATCTGCGCAACGTGATAAGAAAACCTGGAAAAACCGCCGCTATCCTGCTTGa                                                                                                                         <  2:699803/148‑1 (MQ=255)                           aaGAGATTTTAGAAAAAAGGCTGGACGCTGGAATGACCTTTATGCGTCTGACCTTCTCTCATGGTGACTATGAAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATAc                                                                                                                      <  2:1114083/151‑1 (MQ=255)                                    tAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTcc                                                                                                             >  1:1384319/1‑151 (MQ=255)                                     agaaaaaaTGCGGGACGCTGGCATAAACGTTAGGCGTCTGACCTTCTCTCGGGGTGACTATGAAGAACGCGGTCAGCGCATGAAAATTCTGCGCAAGGTGAGGAAAAAAATGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTc                                                                                                              <  2:856676/145‑1 (MQ=255)                                           aaTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATc                                                                                                      >  1:925887/1‑151 (MQ=255)                                             tGCGGGACGAGGGCATGAAGTTTATGCGTCTGAACTTCTCGCAGGGTGACTATGCAAAACCCGGTCAGCGAATTAAGAACCGGCGAAAGGTGATAAACAAAATTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATcc                                                                                                     <  1:445211/150‑1 (MQ=255)                                              gcgggACGCGGGCTGGAACGTTAGGCGTCTGATCTTCTCTCATGTTAACTAGGCAGAACCCGGTCAGCGCATTCAGAATCGGCGCAAGGTAATGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATcc                                                                                                     <  2:387305/146‑1 (MQ=255)                                                  gACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCAcc                                                                                                >  2:108572/1‑150 (MQ=255)                                                           aTGAACGTTATGCGTCTGACCTTCTCTCAGGGTGACTATGCAAAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAAAAAACCGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACt                                                                                        <  2:852741/149‑1 (MQ=255)                                                              aaCGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGa                                                                                     >  2:1402791/1‑149 (MQ=255)                                                                           tGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGt                                                                      >  2:313815/1‑151 (MQ=255)                                                                                ttCTCTCGTGGTGACTAAGCAAAACAGGGTCAGCGCATAAAGAATCGGCGCAAGGTAAGGAAAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacga                                                                         <  1:715637/143‑1 (MQ=255)                                                                                     tcATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGa                                                              >  1:88085/1‑149 (MQ=255)                                                                                       aTGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                          >  1:759741/1‑151 (MQ=255)                                                                                         ggTGACTATGCAGAACACGGTCAGCGCATTCAAAATCTGCGCAACGTGATGAGAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTg                                                        <  2:1134259/151‑1 (MQ=255)                                                                                           tGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGt                                                      >  1:350551/1‑151 (MQ=255)                                                                                             aCTAGGCAAAACCCGGCCAGCGCATTCAAAATCTGCGCAGCGTGAAGAAAAAACCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTCCTCTGAAAGCTGGt                                                      <  2:531802/149‑1 (MQ=255)                                                                                                aggaaaaaCCCGGCCAGCGAATTAAAATTCGGCACAACGTGAGGAGAAAACCGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCa                                                    <  1:1416666/142‑1 (MQ=255)                                                                                                    agaccgcGGTCGGCGCTTTCAAAATCGGCGAAACGTGAGGAGAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAg                                                   <  2:369836/139‑1 (MQ=255)                                                                                                    aGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAAATGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCtt                                              >  1:603508/1‑150 (MQ=255)                                                                                                         acGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACttt                                        >  2:1255465/1‑151 (MQ=255)                                                                                                           ggTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTc                                       >  2:383960/1‑150 (MQ=255)                                                                                                           ggTCAGCGCATTAAGATTCGGCGCAGCGTGATGAAAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTca                                      <  2:1079023/151‑1 (MQ=255)                                                                                                           ggTAAGCGCATTAAGAATCGGCCAAGGGTGAGGAGAAAACGTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTca                                      <  2:368619/151‑1 (MQ=255)                                                                                                            gccAGCGCATAAAAAATCGGCGAAACGTGATGAAAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcac                                     <  2:451701/149‑1 (MQ=255)                                                                                                            gTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcac                                     >  2:1095271/1‑151 (MQ=255)                                                                                                                 cgcATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTg                                >  2:725234/1‑151 (MQ=255)                                                                                                                  gcATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTg                                >  2:1345237/1‑150 (MQ=255)                                                                                                                     ttCAAAATCTGCGAAAGGTGAGGAGAAAAATTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                              <  2:1358668/149‑1 (MQ=255)                                                                                                                      tCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATaa                            >  1:807794/1‑150 (MQ=255)                                                                                                                          attcggcgaAACGTGAGGAGAAAAAGTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAAt                          <  2:1312950/143‑1 (MQ=255)                                                                                                                          aaTCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTg                       >  2:721757/1‑151 (MQ=255)                                                                                                                          aaTCTGCGCAACGTGAGGAGAAAACCGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCt                        <  1:34834/150‑1 (MQ=255)                                                                                                                             cTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGtt                     >  2:1381776/1‑150 (MQ=255)                                                                                                                              tGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTAt                   >  2:632774/1‑151 (MQ=255)                                                                                                                                cgcAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCg                 >  1:234511/1‑151 (MQ=255)                                                                                                                                 gcAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCgg                >  2:134700/1‑151 (MQ=255)                                                                                                                                   aaggtgaacagaaaaaCGGGAAACACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGtt                     <  1:1366541/132‑1 (MQ=255)                                                                                                                                   aaCGTGATGAACAAAACTGGAAAAAACGCCGCTATCCTGCCAGAACCCATAGGTCCGGCAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTCTCTCTGAAACCTGGTCAGACCTTTACTTTCACCATTGATAAATCTGCTATTgg                >  1:1302770/1‑149 (MQ=255)                                                                                                                                      gtgaggaaaaaaaCTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCaa             <  2:1377991/140‑1 (MQ=255)                                                                                                                                      gTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAAc            >  2:395868/1‑150 (MQ=255)                                                                                                                                            aGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGaa      >  2:1262176/1‑150 (MQ=255)                                                                                                                                               aaaaCTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATgg  >  1:457055/1‑151 (MQ=255)                                                                                                                                               aaaaCGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATg   <  2:542196/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                               CACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGG  >  NC_000913/1755721‑1756013

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