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

AGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCAC  >  NC_000913/2421302‑2421586                                                                                                                                              |                                                                                                                                                aGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGTTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGCTCAACTACCCCGCCGa                                                                                                                                        <  2:313367/151‑1 (MQ=255)    aaaaCAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATa                                                                                                                                      >  2:444809/1‑150 (MQ=255)        caacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcg                                                                                                                                 >  1:5090/1‑151 (MQ=255)          acGGAGAACCTGCAATCGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcg                                                                                                                                 >  2:791575/1‑149 (MQ=255)           cGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgcgc                                                                                                                              >  1:226742/1‑151 (MQ=255)            ggAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgcgc                                                                                                                              >  1:673468/1‑150 (MQ=255)            ggAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgcgc                                                                                                                              >  1:492408/1‑150 (MQ=255)                aaCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                           >  1:185596/1‑149 (MQ=255)                       aagggCACACCCTGCCGTTATTATTCATAAAAAGACCCTTCTTTTGCGTACGTTTATCGAAATAAAGGAAAAAAAAATAAACACCCGTAAGGATTTTGTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                                                                  <  2:3148/148‑1 (MQ=255)                        atgccacccccTGCCGTTATTATTCTAAAAAAAACCCTTTTTTTGCGTGTTTTTACCGGAATAAAGAAAAAAAAAATTAACAACCGTAGGGATATTTTTACAAATGGGCCGATCCCCCACCCCGCCGATAAAGCGCGCACTAGCGAAg                                                                                                                    <  2:559697/139‑1 (MQ=255)                                 cctcccttttttattCTTAAAAAAACCTTTCTTTTGCGTACTTTTATCGAAATTAAGAAAAAAAAAATAACCACCCGTCAGAATATTTTTACAAATGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                         <  1:964858/140‑1 (MQ=255)                                 ccgccctttattttTCATAAAAAAACCCTTCTTTTGCGTACTTTTTTGGAAATAAAGAAAAAAAAAATAAACACCCGTCAGAATATTTTTATAAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                           <  2:997143/140‑1 (MQ=255)                                 ccTGCCGTTATTATTCTAAAAAATACCGTTCGTTTGCTTCTGTTTATAGAAATGAAGAAAAAAAAAATAAAAAACGGTAAGGAAGTTATAATCAGTGAGACATCCCAAACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATc                                                                                                               <  2:802760/144‑1 (MQ=255)                                 ccTGCCGTTATTATTCGAAAAAAGACCCTTCGTTTGCGTACTTTTATCGAAATAAAGGAAGAAGAAATTAACACCCGTCAGAATTTTTTTATAAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                         <  1:729747/150‑1 (MQ=255)                                         tattatTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                                                                                >  1:97952/1‑151 (MQ=255)                                                   ataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                      >  1:406699/1‑151 (MQ=255)                                                          aCCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                >  1:53826/1‑150 (MQ=255)                                                              tcttttGCGTACTTTTATCGAAATTAAGGAAGAAGAAATAAACACCCGTCAGGATATTGTTACAAATGTGACCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaaa                                                                            <  1:248596/147‑1 (MQ=255)                                                                 tttGCGTACTTTTATGGAAATTAAGAAAAAAAAAATAAACACCCGTCAGGATATTGTTAACAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                             <  1:782967/146‑1 (MQ=255)                                                                 tttGCGTACTTTTATCGAAATAAAGGAAGAAAAAATAAACACCCGTCAGGATTTTGTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACa                                                                          <  1:829843/149‑1 (MQ=255)                                                                 tttGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                        >  1:911511/1‑151 (MQ=255)                                                                    gcgtttttttATCGAAATAAAGAAAGAAAAAATAAAAACCCGTAAGGATATTTTTAAAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                       <  2:865576/142‑1 (MQ=255)                                                                       tttgTTTATCGAAATTAAGAAAAAAAAAATAAACACCCGTGAGGATATTGTTATAAATGGCACGCTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCa                                                                   <  1:832114/147‑1 (MQ=255)                                                                       tACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                    >  2:554627/1‑149 (MQ=255)                                                                       tACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                    >  1:906423/1‑149 (MQ=255)                                                                           tttATGGAAATAAAGAAAAAAAAAATAAACAACCGTCAGGATATTTTAATCAATGTGCCAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                 <  2:267772/148‑1 (MQ=255)                                                                           tttATCGAAATAAAGAAAAAAGAAATAAACACCCGTCAGGATTTTTTTATAAATGTGACGATCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                               <  2:432716/150‑1 (MQ=255)                                                                            ttATCGAAATTAAGAAAAAAAAAATTAACACCCGTCAGGATATTTTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg                                                              <  1:931402/150‑1 (MQ=255)                                                                            ttATCGAAAAAAAGGAAAAAAAAATAACCACCCGTCAGGATATTTTTACCAAGGTGACTATCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg                                                              <  2:663808/150‑1 (MQ=255)                                                                               tCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                             >  1:92179/1‑148 (MQ=255)                                                                                 gaaaTAAAGAAAAAAAAAATAAACACCCGTCAGGATATTGTTATAAATGTGCCGACCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                        <  1:412658/149‑1 (MQ=255)                                                                                   aaTTAAGAAAGAAAAAATAAACAACCGTCAGGTTTTTGTTACAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                       <  1:185193/150‑1 (MQ=255)                                                                                   aaTTAAGAAAAAAAAAATAAACACCCGTAAGGATATTGTTACAAATGTGACGCTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                       <  2:923640/150‑1 (MQ=255)                                                                                   aaTTAAGAAAAAAAAAATAAACACCCGTAAGGATATTGTTAAAAATGGGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga                                                          <  1:458839/147‑1 (MQ=255)                                                                                      tAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                   >  2:989593/1‑151 (MQ=255)                                                                                      tAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                   >  1:391129/1‑151 (MQ=255)                                                                                       aaGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                  >  1:612654/1‑151 (MQ=255)                                                                                       aaGAAAAAAAAAATAAACACCCGAAAGGATTTTTTTACAAATGTAACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                  <  2:82744/151‑1 (MQ=255)                                                                                           aagaagaaATTAACACCCGCCATGTTATTTTTATCAATGTGACCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                              <  2:317714/151‑1 (MQ=255)                                                                                            agaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                              >  1:426693/1‑150 (MQ=255)                                                                                             gaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                            >  1:764816/1‑151 (MQ=255)                                                                                              aagaaATAAACACCCGTCAGGATTTTTTTATCAATGTGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtctc                                           <  2:428695/151‑1 (MQ=255)                                                                                                aaaaTAAACACCCTTAAGGATTTTTTTATCAATGTACCAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                         <  2:711820/150‑1 (MQ=255)                                                                                                 aaaTAAACAACCGTAAGGATTTTTTTATAAATGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                        <  2:856520/151‑1 (MQ=255)                                                                                                  aaTTAACAACCGCCGGGATATTTTTATCAATGCGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                       <  1:263100/151‑1 (MQ=255)                                                                                                    ttAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCt                                     >  1:855311/1‑151 (MQ=255)                                                                                                    ttAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCt                                     >  2:267165/1‑151 (MQ=255)                                                                                                      aacaacCGTCAGGATATTGTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                   <  1:510855/151‑1 (MQ=255)                                                                                                        caacCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGa                                  >  2:986143/1‑150 (MQ=255)                                                                                                         acccgaaaGGATATTTTTAAAAAGGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCt                                     <  2:723013/139‑1 (MQ=255)                                                                                                            cGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAt                             >  2:884724/1‑151 (MQ=255)                                                                                                                       ttttaataaTGTTACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAg                     <  2:138371/141‑1 (MQ=255)                                                                                                                           ataaaTGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgt                <  1:903627/146‑1 (MQ=255)                                                                                                                           aTCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg               <  1:879206/150‑1 (MQ=255)                                                                                                                              aaTGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg           >  2:243238/1‑151 (MQ=255)                                                                                                                               aTGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg           >  2:418492/1‑150 (MQ=255)                                                                                                                                 gtgACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg           <  1:140738/148‑1 (MQ=255)                                                                                                                                       aTCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCAc  >  1:925895/1‑151 (MQ=255)                                                                                                                                         ccACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCAc  <  1:848795/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                AGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCAC  >  NC_000913/2421302‑2421586

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