Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,443	0	A	C	44.6%	101.5 / 10.0	56	Y40S (TAC→TCC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base A (29/2);  new base C (0/25);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 6.30e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGA  >  NC_000913/2421305‑2421590                                                                                                                                              |                                                                                                                                                    aaaaCAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATa                                                                                                                                          >  2:810573/1‑150 (MQ=255)  aaaCAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATaa                                                                                                                                         >  1:555735/1‑150 (MQ=255)   aacaacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATaa                                                                                                                                         >  2:106880/1‑149 (MQ=255)   aacaacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATaa                                                                                                                                         >  1:517905/1‑149 (MQ=255)     caacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAgcg                                                                                                                                     >  2:464295/1‑151 (MQ=255)      aacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAgcgc                                                                                                                                    >  1:912278/1‑151 (MQ=255)            gaACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                               >  2:604538/1‑150 (MQ=255)                cTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                           >  1:562092/1‑150 (MQ=255)                cTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                           >  1:185233/1‑150 (MQ=255)                            aaCCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                               >  1:932278/1‑150 (MQ=255)                             aCCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                             >  2:807473/1‑151 (MQ=255)                              ccTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATg                                                                                                            >  1:494813/1‑151 (MQ=255)                                tGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                          >  1:701648/1‑151 (MQ=255)                                       attattCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                                   >  2:15877/1‑151 (MQ=255)                                        tttttagaaaaaaGACCC‑TTCGTTGGCGTACTTTTTTGGAAATTAAGA‑AAAAAAAAACAAACCACCGCCAGAATTTTGTTAAAAAGGTGC‑CGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGaa                                                                                                     <  2:701878/138‑1 (MQ=255)                                        ttattCGTAAAAAGACCC‑TTCGTTTGGTTATTTTTAT‑GGAAATAAGGAAAAAAAAAATAAACACCCGTAAGTATTTTGTTATAAATGTGC‑CGTCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                                   <  2:145991/150‑1 (MQ=255)                                           ttCGAAAAAAAAACC‑TTCTTTTGCGTATTTTTATCGAAATAAAGA‑AAAAAAAAATAAACAACCGTAAGAATTTTTTTATAAATGTGA‑CGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                                <  2:1016257/150‑1 (MQ=255)                                            ttattaaaaaaaCCTTTCTTTTGCTTTCTTTTATCGAAATAAAGA‑AAAAAAAAATAAACCCCCGTCAGAATTTTTTTATAAATGCGA‑CCCTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                                                                 <  2:374882/143‑1 (MQ=255)                                                ataAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                          >  2:663652/1‑151 (MQ=255)                                                ataAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                          >  1:7079/1‑151 (MQ=255)                                                aaaaaGAACC‑TTCTTTTGCGTACGTTTATCGGAATTAAGA‑AAGAAAAAATAAACACCCGGCAGGATATTGTTATAAATGTGA‑CGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                              <  1:317264/145‑1 (MQ=255)                                                   aaagaccttttttttGCTTCCTTTTATCGAAATTAAGA‑AAAAAAAAATTAACACCCGCCAGGATATTTTTATCAATGTGA‑CAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                          <  1:66064/137‑1 (MQ=255)                                                    aaaccc‑TTCTTTTGCGTAGTTTTATGGAAATAAAGA‑AAAAAAAAATAAACACCCGTAAGAATTTTTTTATCAATGTGA‑CACCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                      <  1:546057/147‑1 (MQ=255)                                                           ttCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                                 >  2:796741/1‑150 (MQ=255)                                                            tCGTTTGGGTATGTTTATCGGAATAAAGG‑AAGAAAAAATAAACACCCGTCAGAATTTTTTAATAAATGTGC‑CAATCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                     <  2:865234/145‑1 (MQ=255)                                                                  gcttgtttttATGGAAATTAAGG‑AAAAAAAAATTAACACCGGCGAGGATTTTGTAAT‑AATTTGAGCGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                          <  1:411227/143‑1 (MQ=255)                                                                  gCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattatt                                                                         >  2:124644/1‑151 (MQ=255)                                                                   cGTACTTTTATCGAAATAAAGA‑AAAAAAAAATTAACAACCGTCAGGATTTTGTTATCAATGTGA‑CGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                          <  1:254625/149‑1 (MQ=255)                                                                         tttATGGAAATAAGGG‑AAAAAAAAATAAACACCCGTCGGAATTTTTTTATCAATGTGA‑CGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCa                                                                    <  1:892412/149‑1 (MQ=255)                                                                         tttAGCGAAATAAAGA‑AAGAAAAAATTAACAACCGTCAGGATATTGTTATCAAGGTGA‑CGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                     <  1:614308/148‑1 (MQ=255)                                                                         tttAGCGAAATAAAGA‑AAAAAAAAATAAACACCCGTCAGAATTTTGTTATAAATGGGA‑CGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                   <  1:744989/150‑1 (MQ=255)                                                                          ttATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                 >  2:207088/1‑151 (MQ=255)                                                                             tCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                 >  2:1012545/1‑148 (MQ=255)                                                                               ggAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                               >  1:111695/1‑148 (MQ=255)                                                                                 aaTAAAGA‑AAAAAAAAATTAACACCCGTAAGAATATTTTTATAAATGTGC‑CAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag                                                             <  1:908850/148‑1 (MQ=255)                                                                                 aaTAAAGA‑AAAAAAAAATAAACACCCGTAAGGATTTTTTTATAAATGTGA‑CGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga                                                              <  2:659443/147‑1 (MQ=255)                                                                                    tAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                       >  2:987137/1‑151 (MQ=255)                                                                                     aaga‑aagaAGAAATAACCACCCGACAGAATATTTTAATCAATGTAC‑CAATCCCCCACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga                                                              <  2:185228/143‑1 (MQ=255)                                                                                     aaGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                      >  1:606485/1‑151 (MQ=255)                                                                                     aaGA‑AAAAAAAAATAAACACCCGTGAGATTTTTTTTAAAAATGTGA‑CGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATaa                                                        <  1:462187/149‑1 (MQ=255)                                                                                          aaaaaaaaaTAAACACCCGTAAGTATATTTTTATAAATGTGA‑CGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGt                                                    <  2:555719/143‑1 (MQ=255)                                                                                           agaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtc                                                 >  1:936598/1‑151 (MQ=255)                                                                                               gaaATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                             >  1:271819/1‑151 (MQ=255)                                                                                                   ttAACAACCGTCAGGATTTTTTTAAAAGTGTAA‑CGCCCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGc                                          <  1:439933/150‑1 (MQ=255)                                                                                                       cacccggaaGGATATTGTTATAAATGGGC‑CGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGa                                      <  1:90563/142‑1 (MQ=255)                                                                                                          ccgaaagaatatTTTTATAAATGTGC‑CAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                       <  2:574479/138‑1 (MQ=255)                                                                                                                 gATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                           >  1:843114/1‑151 (MQ=255)                                                                                                                   tattttaaaaattgtgA‑CACTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACt                            <  2:469900/136‑1 (MQ=255)                                                                                                                     ttttcagcattgtgA‑CCGTCCCCCACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACt                            <  1:486042/136‑1 (MQ=255)                                                                                                                        ttaaaaaTGTGA‑CCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtg                     <  2:822378/145‑1 (MQ=255)                                                                                                                           tCAATGTGA‑CGCTCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGc                  >  2:173590/1‑150 (MQ=255)                                                                                                                             aaTGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg               >  1:790274/1‑151 (MQ=255)                                                                                                                                     cGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGatat           <  1:560526/148‑1 (MQ=255)                                                                                                                                         ccACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGTGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGt    <  1:397124/151‑1 (MQ=255)                                                                                                                                         ccACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGAAc      <  1:340891/149‑1 (MQ=255)                                                                                                                                            cTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGa  >  1:1029179/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                    AAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGA‑CGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGA  >  NC_000913/2421305‑2421590

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