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.9%	79.6 / 32.8	56	G58G (GGT→GGA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (28/4);  new base A (0/24);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 9.40e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGT  >  NC_000913/1755724‑1756020                                                                                                                                                    |                                                                                                                                                      cATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTaa                                                                                                                                                     >  2:390762/1‑150 (MQ=255)     ggACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccg                                                                                                                                                >  2:871740/1‑151 (MQ=255)      gACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccgc                                                                                                                                               >  1:802010/1‑151 (MQ=255)           aaaaCCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCt                                                                                                                                           >  2:856616/1‑150 (MQ=255)               ccGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCt                                                                                                                                      >  1:636532/1‑151 (MQ=255)                cGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTg                                                                                                                                     >  2:284687/1‑151 (MQ=255)                cGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTg                                                                                                                                     >  1:541811/1‑151 (MQ=255)                        aaGAGATGTTAGCAAAAATGCGGGCCGCTGGAATAAACGTTATGGGTCTGAACTTCCCTCAGGGTGACTATGCAAAACAGGGTCAGCGCATTCAGAATCTGCGCACCGTGATGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATa                                                                                                                              <  2:474941/150‑1 (MQ=255)                           agaTGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGGGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCa                                                                                                                           >  2:477734/1‑150 (MQ=255)                           agaTGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCa                                                                                                                           >  2:896494/1‑150 (MQ=255)                              tGTTATCTAAAAGGCGGGACGTGGGCATAAACGTTAGGCTTCTAATCTTCTCTAATGGTGACTATGAAAAACCCGGTAAGCTAATTCAAATTCTGCGCAACGTGATGAGAAAATGTTAAAAAACCGCCGCTATCCTGCtt                                                                                                                                  <  1:604861/140‑1 (MQ=255)                              tGTTAGCAAAAATGCTGGCCGCGGGCATGAGCGTTATGCGTCTGAACTTCTCTCAGGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGAAAAGCTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAg                                                                                                                        <  2:118763/150‑1 (MQ=255)                              tGTTAGAAAAAATGCTGGCGGCTGGCAGAAGCTTGATGCGTCTAAACTCCTCTCAGGGTGACTAGGAAAAACGCGTTCGGCTCATTCAAAGTCTGCGCAAGTTAAGGACAAAAGCTGGAAAATCCGCCGCTATCCTGCTTGATACCAAAg                                                                                                                        <  2:413939/150‑1 (MQ=255)                                ttAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACCCGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAAAAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAgg                                                                                                                       <  2:220066/149‑1 (MQ=255)                                 tagaaaaaaTGCTGGACGTTGGCATGAATTTTTTGCTTTTGATTTTCTCTCATGGTGACTAGGAAAAACACGGTCAGCGCTTTCAAATTCTGAGAAACGTGATGAGAAAATCTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTc                                                                                                                     <  1:661537/145‑1 (MQ=255)                                       aaaTGCTGGCGGCGGGCATGAACGTTATGGGTCTGACCTCCTCTCATGGTGCCTATGAAGAACACGGCCGGCGCATTCAGATTCGGCGCAAGGTGAGGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAAt                                                                                                              <  2:888706/151‑1 (MQ=255)                                          tGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAGAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCg                                                                                                           >  1:172054/1‑151 (MQ=255)                                             tGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCAc                                                                                                        >  1:845570/1‑151 (MQ=255)                                                aCGCGGGAATGAACGTTATGCGTCTGAACTTCTCTCAGGGTGACTATCCAAAACACGGTCAGCGCATTCAGAATCTGCCAAACGTGATGAGCAAATCTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGc                                                                                                          <  2:554437/146‑1 (MQ=255)                                                 cGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATg                                                                                                    >  1:418909/1‑151 (MQ=255)                                                     ggCATAAACGTTATGCGTCTGAACTTCTCTCGTGGTGACTAGGCAAAACACGTTCAGCGCATTCAGAATCTGCGCAGCGTGAAGAGAAAAACTGATAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGaa                                                                                                  <  1:158863/149‑1 (MQ=255)                                                            aCGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAgg                                                                                         >  1:534880/1‑151 (MQ=255)                                                             ctttATGCGTCTGAACTTCTCTCAGGGTGACTATGAAGAACACGGTCAGCGCATTCAGAATCTGCGCAAGGTGATGAAAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAgg                                                                                         <  1:505375/148‑1 (MQ=255)                                                               ttATGCTTATGACCTTCCCTAAGGGTGAATAGGAAAAACAGGGCCAGCGATTTAAAAATCGGCGAAACGTGAGAAAAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCg                                                                                       <  1:516709/150‑1 (MQ=255)                                                               ttAGGCTTCTGACCTTCTCTCAGGGTGCCTAGGCAGAACCGGGTCAGCGCATTAAGAATCTGCGCACCGTAATGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCg                                                                                       <  2:721055/150‑1 (MQ=255)                                                               ttAGGCGTCTGACCTTCTCTAAGGGGGACAAGGAAAACCACGGTCAGCGCATCCAGATTCTGCCCAACGTGATGAGAAAACGGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGc                                                                                        <  2:15829/149‑1 (MQ=255)                                                                    cGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacg                                                                                 >  1:225850/1‑151 (MQ=255)                                                                       cTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgacg                                                                              >  2:60002/1‑151 (MQ=255)                                                                        tGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACGATGAAACTGGAAGGCGGTAACGACGt                                                                             >  1:622707/1‑151 (MQ=255)                                                                        tGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGt                                                                             >  2:67044/1‑151 (MQ=255)                                                                         gAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAGCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgacg                                                                              <  1:216840/149‑1 (MQ=255)                                                                            cTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtct                                                                         >  2:336220/1‑151 (MQ=255)                                                                               ctctcATGGTGATTATGCAAACCACGGTCAGCCCATACAGAATCTGCGCAACGTGAGGAAAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtctc                                                                        <  1:779579/149‑1 (MQ=255)                                                                                tctcATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTg                                                                      >  2:874393/1‑150 (MQ=255)                                                                                     tGGTGACTATGAAGACCCCGTTCAGCGCATTCAGAATCTGCGAAACGTGAGGAGAAAACCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaaa                                                                   <  2:281732/148‑1 (MQ=255)                                                                                     gggTGACTAGAAAAAACACGGTCAGCCCATTCAAAATCTGCACAAGGTGAGGAAAAAAATGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAg                                                                  <  1:451635/148‑1 (MQ=255)                                                                                      ggTGACTATGCAAAACACGGTCAGCACATTCAAAATCTGCGAAACGTGATGAAAAAAACGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaaa                                                                   <  1:745464/147‑1 (MQ=255)                                                                                      ggTGACTATGCAAAACACGGCCAGCCCATTAAAAATCGGCGCAACGTGATGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCt                                                                <  2:185776/150‑1 (MQ=255)                                                                                             aTGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGAc                                                        >  2:212810/1‑151 (MQ=255)                                                                                                 aaaaCACGGTCAGCGCATCAAAAATCTGCGCAACGTGATGAGAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGAcc                                                       <  2:190249/146‑1 (MQ=255)                                                                                                 aGACCAGGGCCAGCGCATACAAAATCGGCGAAACGTGAAGAAAAAACCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCtt                                                     <  1:365524/150‑1 (MQ=255)                                                                                                   aaCACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTAc                                                  >  2:808201/1‑151 (MQ=255)                                                                                                   aaCACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCtt                                                     <  2:901181/148‑1 (MQ=255)                                                                                                      acGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACttt                                               >  2:560982/1‑151 (MQ=255)                                                                                                       cGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTc                                              >  1:375519/1‑151 (MQ=255)                                                                                                        ggTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTc                                              >  1:678399/1‑150 (MQ=255)                                                                                                         gTCAGCGAATTCAGAATCTGCGCAACGTGAGGAGAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcac                                            <  2:207355/151‑1 (MQ=255)                                                                                                               ccATTAAAAATCTCCACAACGTAAGAACCAAATTTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcacca                                          <  2:774039/146‑1 (MQ=255)                                                                                                                 aTTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                                     >  1:761554/1‑150 (MQ=255)                                                                                                                     aaaaTCTGCGCACCGTGATGAGAAAAACTGTAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAAt                                 <  2:263783/148‑1 (MQ=255)                                                                                                                       aacctgcaaaaCGTGAGGAGCAAATCTGGAAAAGCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATc                                <  2:97896/140‑1 (MQ=255)                                                                                                                       aaTCTGCGCAATGTGAAGAGAAAAACGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTg                              <  2:912656/151‑1 (MQ=255)                                                                                                                             cgcAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCg                        >  2:452518/1‑151 (MQ=255)                                                                                                                                aaCGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGc                      >  2:862690/1‑150 (MQ=255)                                                                                                                                aaCGTGAGGAGAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTa                           <  1:233134/145‑1 (MQ=255)                                                                                                                                            aaaaCTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATgg         >  2:220563/1‑151 (MQ=255)                                                                                                                                                    tAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGt  >  2:42598/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGT  >  NC_000913/1755724‑1756020

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