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

ATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGC  >  NC_000913/1755725‑1756006                                                                                                                                                   |                                                                                                                                        aTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTaaa                                                                                                                                      >  2:670546/1‑150 (MQ=255)     gACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccgc                                                                                                                                 <  1:333695/151‑1 (MQ=255)       ccGAAAACCGAATCTGAAGAGATGTTGGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAGCTGGTAAAAccgcc                                                                                                                                <  2:67117/150‑1 (MQ=255)          aaacCCGAATCTGAAGAGATTTTAGCTAAAATGCTGGACGCTGGCATGAACTTTATGCGTCTAAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTCCGAAACGTGATAAGCAAAACTGGAAAAACcgccgc                                                                                                                              <  2:87063/149‑1 (MQ=255)                       aagatttGTTAGAAAAAATGTTGGCCGCGGGAATAAACTTTTGGCTTCTAACCTTCTCTCAGGGTGACAAGGCAGAACCCGGTCAGCCCTTTCAGAATCGGCACAACGGGATAAGAAAACGTGAAAAATCCGCCGCTATCCTGCTTGATa                                                                                                                <  1:287522/144‑1 (MQ=255)                       aaGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATAc                                                                                                               >  1:56841/1‑151 (MQ=255)                             tGTTAGCTAAATGGCTGGCCGCCGGCATGAACGTTCTGCGTCTGACCTTCTCCCATGGTGATTATGCAAACCACGATCACCGCTTTCAAAATCTCCGAATCATAATGAACAAAACTGATAAAACCGCCGCTATCCTGCTTGATAcc                                                                                                              <  2:12717/146‑1 (MQ=255)                               taaccaaaaaGGCGGGCCGCGGGAATAAACTTTAGGCGTCAAAACTTCTCAAAGGGTGACTATGCAAAACACGGTCAGCGCATTCAGAATCTGCGCAACGGGATGAAAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCaa                                                                                                            <  1:599388/140‑1 (MQ=255)                                 aGCAAAAATGCTGGACGCGGGCATGAACTTTATGCGTCTGATCTTCTCTCATGGTGACTATGCAAAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAGGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTc                                                                                                       <  1:275470/149‑1 (MQ=255)                                    tAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGa                                                                                                   >  2:624623/1‑150 (MQ=255)                                       aaTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATc                                                                                               >  2:248181/1‑151 (MQ=255)                                       aaTGCTGGACGCTGGAAAGAACGTTCTGCGGCTCACCTTCTCGGATGGTACCTAGACAAAACACGGCCTGCGCATTCATATTCCGCGCAGGGTGAGAAGAAACGCGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATc                                                                                               <  2:741443/151‑1 (MQ=255)                                         tGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCg                                                                                             >  1:38407/1‑151 (MQ=255)                                         tGCGGGCCGCTGGAATGAACGTTATGCGTCTGACCTTCTCTCAGGTTGACTATGAAGAACGGGGTCAGCGCATTAAAAATCGGCCCAACGTGATGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATcc                                                                                              <  1:392868/150‑1 (MQ=255)                                                       aTGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTgg                                                                               >  1:645856/1‑151 (MQ=255)                                                       aTGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTgg                                                                               >  2:483703/1‑151 (MQ=255)                                                          aaCGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGa                                                                              >  1:760604/1‑149 (MQ=255)                                                              ttATGCGTCTAATCTTCTCTCATGGTGACTATGAAAACCACGGTCAGCGCATTAAAATTCGGCGAAACGTGAGGAAAAAAACGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGc                                                                          <  2:592822/149‑1 (MQ=255)                                                                        aaccTTCTCTCGGGGTGACTATGCAAAACCCGGTCACCGCATTAAGAATCGGCGCAACGTAATAAAAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgac                                                                 <  2:28301/145‑1 (MQ=255)                                                                         aaCTTCTCTAATGGAGACTATGCAGAACACCGTCAGCGCATACAAAATCTGCGCAACGTGCTGCGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCTGAAATCCGAACCATGCAACTGGAAGGCGGTCacgac                                                                 >  2:565680/1‑147 (MQ=255)                                                                              ctctcATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTg                                                        >  1:677660/1‑151 (MQ=255)                                                                                      gTGACTATGCAAAACACGGCCAGCGCATTAAGATTCTGCGCAACGTGATGAGAAAAACTGTAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTgg                                                <  1:828210/151‑1 (MQ=255)                                                                                        aactagacaaaaCACGGCCAGCCCATCCAGAATCTGCCAAACGTGATAAGAAAAACTTAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCCGGTc                                              <  2:677476/141‑1 (MQ=255)                                                                                                aGAACCGGGTCAGCCCATACAGAATCTGCGCAAGGAAAAGAGAAAAACGGAAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCtt                                       <  2:1928/150‑1 (MQ=255)                                                                                                  aaCACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTAc                                    >  1:224654/1‑151 (MQ=255)                                                                                                       ggTCCGCACATTAAGAATCTGCGAAACGTGAGAAGAAAAACGGAAAAATCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACttt                                 <  1:385558/149‑1 (MQ=255)                                                                                                                    aGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATc                  >  2:463969/1‑151 (MQ=255)                                                                                                                          tGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTAt            >  1:449668/1‑151 (MQ=255)                                                                                                                             gcAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCgg         >  2:375557/1‑151 (MQ=255)                                                                                                                               acggTGATGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCa       <  1:451996/148‑1 (MQ=255)                                                                                                                                    gaagaaaaaaaCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGc  <  2:324166/144‑1 (MQ=255)                                                                                                                                        aaaaaaaCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAg   <  1:101007/143‑1 (MQ=255)                                                                                                                                          aaaaaCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGc  <  1:386841/144‑1 (MQ=255)                                                                                                                                           aaaaCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAg   <  1:805623/143‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                        ATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGC  >  NC_000913/1755725‑1756006

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