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.1%	39.8 / 21.6	38	G58G (GGT→GGA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (18/4);  new base A (0/16);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 2.25e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAAC‑TT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTG  >  NC_000913/1755729‑1756016                                                                                                                                                 |                                                                                                                                                   gACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAAC‑TT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑ccgc                                                                                                                                           >  2:204507/1‑151 (MQ=255)               tctaaagattttttAGCAAAAATGCTGGACGCGGGCATGAACTTTATGCGTCTAACC‑TC‑CTCTCATGGTGACTATGCAAAACACGTTCACCCCATCCAAATTCTGCCCAGCGTGAGGAGCAAACCTGAAAAAA‑CCGCCGCTATCCTGCtt                                                                                                                              <  1:1046615/138‑1 (MQ=255)                   aaGAGATGTTAGCAAAAAGGGTGGACGCTGGCATGAACGTTATGCGTCT‑AAC‑TTCCTCTCAGGGAGACTGTGCAGAACACGGCCGGCGTATCAAAAACCGGCAAACCGTAATGAGAAAAACTGGAAAAA‑CCGCCGCTATCCTGCTTGAt                                                                                                                           <  2:251991/149‑1 (MQ=255)                        aTGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAAC‑TT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCaa                                                                                                                      >  2:722441/1‑149 (MQ=255)                            tAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAAC‑TT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTcc                                                                                                                >  2:185438/1‑151 (MQ=255)                                tAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAAC‑TT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGa                                                                                                             >  2:457930/1‑150 (MQ=255)                                   aaTGCTGGACGCTGGCATGAACGTTATGCGTCTGAAC‑TT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATc                                                                                                         >  2:257917/1‑151 (MQ=255)                                           aCGCTGGCATGAACGTTATGCGTCTGAAC‑TT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCAt                                                                                                 >  2:145670/1‑151 (MQ=255)                                            cgctgcaaggcagtttATGCGTC‑GAACTTT‑CTCTCGTGGTACTTACGCAAAACCGGGACCGCCAATTAATAATCGGCGACAAGTGATGAAAAAAACTGGTAAAT‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCAt                                                                                                 <  1:427359/136‑1 (MQ=255)                                            cgcgGGCATGAACTTTATGCGTCTGAAC‑TT‑CTCTCATGGTGACTAGGCAAAACGCGGTCGGCGCTTTAAGAATCTGCGAAACGTGATGAGAAAACCTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATg                                                                                                <  2:435388/151‑1 (MQ=255)                                                   aTGAACGTTATGCGTCTGAAC‑TT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTgg                                                                                         >  1:911877/1‑151 (MQ=255)                                                      aaCGTTATGCGTCTGAAC‑TT‑CTCTCATGGTGACTATGAAGAACCCGGTCAGCGCATTAAAAATCGGCGCAACGTGATGAGAAAAACGGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGaa                                                                                       <  2:469992/150‑1 (MQ=255)                                                               cGTCTGAAC‑TT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacg                                                                             >  1:657281/1‑151 (MQ=255)                                                                gTCTGAAC‑TT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacg                                                                             >  1:869774/1‑150 (MQ=255)                                                                  cTGACC‑TT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAGGTGATGAGAAAAACTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacg                                                                             <  2:342050/148‑1 (MQ=255)                                                                         tt‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtct                                                                     >  2:9497/1‑150 (MQ=255)                                                                             tctcATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTg                                                                  >  2:614945/1‑150 (MQ=255)                                                                                   ggtaattATGCAGAACGCGGCCAGCGCATTAAAAATCTGCGAAAGGTGATGAAAAAAACTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaaa                                                               <  1:854147/141‑1 (MQ=255)                                                                                     tGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGt                                                         >  1:611476/1‑151 (MQ=255)                                                                                      gACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTc                                                        >  2:350594/1‑151 (MQ=255)                                                                                           ggCAAAACACGGCCAGCCCATTAAGAATCTGCGCACCGTAATGAGAAAAACGGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAg                                                      <  1:27723/147‑1 (MQ=255)                                                                                              aaagcccGGTCTGCCCATCAAAAATGTGCGCAGCGTTATATACAACGCTGTAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGt                                                         <  2:1093461/136‑1 (MQ=255)                                                                                              aaaaCACGGTCAGCACATTCAAAATCTGCACAACGTAATGAGCAAAACTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCt                                                  <  1:62414/147‑1 (MQ=255)                                                                                              aGAACACGGCCAGCGCTTCAAGAATCTGCCAAAGGTGAGAACAAAAACGGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCt                                                  <  1:552156/149‑1 (MQ=255)                                                                                                           cgcATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTg                                   >  1:17629/1‑151 (MQ=255)                                                                                                             cATTCAAATTCGGCGCAACGTGATGAGCAAAACTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGa                                  <  2:634791/150‑1 (MQ=255)                                                                                                                  aGAATCTGCGCAAGGGAATGAGAAAACCTGATAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcacca                                      <  1:1093094/141‑1 (MQ=255)                                                                                                                    aaTCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTg                          >  1:110664/1‑151 (MQ=255)                                                                                                                             accgggatgaGAAAACCTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCt                           <  2:652074/136‑1 (MQ=255)                                                                                                                             aacgagaggaaaaaaaCTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTAt                      <  1:574770/134‑1 (MQ=255)                                                                                                                                   atgaCCAAAACTGGAAAAT‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCg           <  1:1040533/151‑1 (MQ=255)                                                                                                                                    tgaGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCg           >  1:916503/1‑150 (MQ=255)                                                                                                                                    gaagaaaaaTTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACa              <  1:475722/142‑1 (MQ=255)                                                                                                                                     gaGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGa          >  2:921829/1‑150 (MQ=255)                                                                                                                                        aaaaaTTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACa              <  2:942761/142‑1 (MQ=255)                                                                                                                                         aaaaCTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATgg     >  2:196201/1‑151 (MQ=255)                                                                                                                                         aaaaCTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACa              <  1:35130/142‑1 (MQ=255)                                                                                                                                              cggTAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTg  <  1:550778/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAAC‑TT‑CTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTG  >  NC_000913/1755729‑1756016

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