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

ATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGT  >  NC_000913/1755725‑1756020                                                                                                                                                   |                                                                                                                                                       aTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTaaaa                                                                                                                                                    <  2:255070/151‑1 (MQ=255) aTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTaaa                                                                                                                                                     >  1:704162/1‑150 (MQ=255)   cGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑cc                                                                                                                                                 <  1:416849/151‑1 (MQ=255)         gAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCt                                                                                                                                           >  1:791282/1‑151 (MQ=255)          aaaaCCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCt                                                                                                                                           >  1:783132/1‑150 (MQ=255)                gAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGc                                                                                                                                    >  1:161224/1‑151 (MQ=255)                gAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGc                                                                                                                                    >  1:806182/1‑151 (MQ=255)                       aaGAGATGTTAGATAAAATGCGGGACGCGGGAATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAAAACACGGCCAGCGAATTCAAAATCGGCGCAACGTGATGAGAAAAACTGGAAAAA‑CCGCCGCTATCCTGCTTGAt                                                                                                                               <  1:672867/149‑1 (MQ=255)                            aTGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCaa                                                                                                                          >  1:775690/1‑149 (MQ=255)                                                      cATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTg                                                                                              >  1:384268/1‑151 (MQ=255)                                                       aTGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTgg                                                                                             >  2:482747/1‑151 (MQ=255)                                                       aTGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTgg                                                                                             >  1:81618/1‑151 (MQ=255)                                                              ttctgcgtttaactttCCCTCAGGGTGATTAGGCAAAACAGGGCCACCCCATTAAAATTCTGCGAACCGTAAGAAGAAAAGGTGGAAAAC‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGc                                                                                        <  2:342848/135‑1 (MQ=255)                                                                   cGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacg                                                                                 >  2:592895/1‑151 (MQ=255)                                                                      cTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgacg                                                                              >  2:846881/1‑151 (MQ=255)                                                                         aaCTTCTCTCAGGGTGACTAGAAAAAACACGGACGGCGCTTTAAAAACCGGCACAAGGGAAAAAGAAAAACTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCg                                                                                       <  2:432728/139‑1 (MQ=255)                                                                            ttCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtctc                                                                        >  2:859506/1‑151 (MQ=255)                                                                             tctctcATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTctct                                                                       >  1:616554/1‑151 (MQ=255)                                                                                     ggTGACTATGCAAAACACGGTCAGCGCATTCAAAATCTGCGAAACGTGATGAGAAAAACTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                                 <  1:739290/149‑1 (MQ=255)                                                                                     ggTGACTATGAAAAACACGGTCAGCGCTTTCAAAATCGGCGAAACGTGAGAAAAAAAATTGTAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCt                                                                <  1:752890/150‑1 (MQ=255)                                                                                     ggTGAATATGCAGAACCCGGTCAGCGTTTAAATACTCCGCGCATGGAAAGGAGAAAACTTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTg                                                               <  2:656176/151‑1 (MQ=255)                                                                                      gTGACTATGCAAAACGCGGTCAGCGCATTCAGATTCTGCGAAAGGTGATGAGCAAAACTGTAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                                 <  1:341462/148‑1 (MQ=255)                                                                                       tGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGt                                                             >  2:676725/1‑151 (MQ=255)                                                                                       tGACAATGAAGAACGCGGGCAGCCCATTAAGAATCTCCCCAGAGAGAAGAGAAAAAGGAAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaa                                                                    <  2:493416/144‑1 (MQ=255)                                                                                        gACTATGCAGAACAGGGTCAGCACATTAAGAATCTGCGCAACGTGATGAAAAAAAGTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGt                                                             <  1:601876/150‑1 (MQ=255)                                                                                         aCTATGCAGAACACGTCCAGCTAATTCAGAATCTGCGCAGGGTGATGAGCAAACCGGGAAAAC‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGt                                                             <  2:541609/149‑1 (MQ=255)                                                                                                agaactcaatCACCGCATTCAGATTCTGCCCAACTTGATGAGCAAAACTGGAAATA‑CCGCCGCTATCCTGCTTGATTCCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCtt                                                     <  2:348274/141‑1 (MQ=255)                                                                                                  aaCAGGTTCAGCGCATTAAAAACCCGCACAACGTGATAAACAAAACTATTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTa                                                   <  1:482802/150‑1 (MQ=255)                                                                                                       ggccAGCGAATCAAGAATCTGCGCAAGGTGATGAGAAAAACGGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTc                                              <  1:457473/147‑1 (MQ=255)                                                                                                       ggacAGCGCATTCAGAACCTGCGCAAGGGGATGAGCAAAACTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACttt                                               <  1:539491/146‑1 (MQ=255)                                                                                                                 ttCAAAATCTGCGAAACGAGATAAAAAAAACTGAAAAAC‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcaccac                                         <  2:105600/145‑1 (MQ=255)                                                                                                                  taaGAATCTGCACAAGGTAAAAAAAAAAACTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGa                                      <  2:762402/145‑1 (MQ=255)                                                                                                                   cAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAAt                                 >  2:880649/1‑151 (MQ=255)                                                                                                                    aaattcggcgcAACGTAATGAAAAAAACTGTAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATc                                <  1:805101/144‑1 (MQ=255)                                                                                                                    aaaaTCGGCACAACGTGAAAAAAAAAAC‑GGTAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGa                                      <  2:422251/143‑1 (MQ=255)                                                                                                                    aGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATc                                >  2:127681/1‑151 (MQ=255)                                                                                                                     gAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCt                               >  1:298466/1‑151 (MQ=255)                                                                                                                       aTCGGCACAACGTGAGAAAAAAAACGGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGt                             <  2:513506/151‑1 (MQ=255)                                                                                                                           gcgcAACGTTATGAAAAAAATTGTAAACG‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGTAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGt                             <  2:132988/147‑1 (MQ=255)                                                                                                                               aaCGTGATGAAAAAAGGGGGTAAAG‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTAt                          <  1:335424/146‑1 (MQ=255)                                                                                                                                aCGTGATGAGAAAAACTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCg                        <  1:320739/147‑1 (MQ=255)                                                                                                                                  gggatgaCAAAAATTTGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACa                  <  2:79472/149‑1 (MQ=255)                                                                                                                                  gTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAAc                   >  1:63623/1‑150 (MQ=255)                                                                                                                                     atgaGAAAACCTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAAc                   <  2:146500/147‑1 (MQ=255)                                                                                                                                      tgaGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCg               >  1:806210/1‑150 (MQ=255)                                                                                                                                           aaacCTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAAt           <  1:252990/149‑1 (MQ=255)                                                                                                                                            aaaCTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGt        >  2:675351/1‑151 (MQ=255)                                                                                                                                                   tAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGt  >  2:660387/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                       ATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGT  >  NC_000913/1755725‑1756020

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