Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,443	0	A	C	61.2%	29.9 / 62.2	50	Y40S (TAC→TCC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base A (18/1);  new base C (0/30);  total (18/32)
Fisher's exact test for biased strand distribution p-value = 1.64e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-09

ATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421296‑2421592                                                                                                                                                    |                                                                                                                                                      aTCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAccc                                                                                                                                                    >  2:430824/1‑151 (MQ=255)  tCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAcccc                                                                                                                                                   >  1:68947/1‑151 (MQ=255)  tCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAcccc                                                                                                                                                   >  2:80892/1‑151 (MQ=255)  tCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAcccc                                                                                                                                                   >  1:212157/1‑151 (MQ=255)       aGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgccg                                                                                                                                               >  1:373367/1‑150 (MQ=255)         tAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGAt                                                                                                                                             >  1:222645/1‑150 (MQ=255)                  ggAGACCCTGAAAGGGCACACCCTGCCGTTATTATTCGAAAAAAGACCCTTCGTTTGCGTACGTTTATCGGAATTAAGAAAAAAAAAATAAACACCCGTAAGGATTTTGTTATAAGTGTGACGATCCACTCCCCCGCCGATAAAgcgc                                                                                                                                      <  1:304540/148‑1 (MQ=255)                      acccTGAAAGGCCACACCCGGCCGTTATTATTCTTAAAAAGACCCTTCTTTTGCGTACGTTTATGGGAATTAAGGAAAAAAAAATAAACACCCGTCAGGATATTGTTATAAATGTGCCGACCCACTCCCCCGCCGATAAAGCGCGCACTAg                                                                                                                               <  2:291986/149‑1 (MQ=255)                      aaCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                                 >  1:434512/1‑149 (MQ=255)                        ccgacAATGCCCCACCCTCCCGTTTTTATTCATAAAAAAACCTTTCTTTTGCGTCCTTTTATCAAAATAAAGAAAGAAAAAATAAACACCCGACAGGATTTTTTTATAAATGTGACATCCCACTCCCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                             <  2:421340/147‑1 (MQ=255)                          tGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                            >  1:463656/1‑150 (MQ=255)                             aaTGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                            >  2:177391/1‑147 (MQ=255)                                 acacccccggCCGTTTTTTTTCATAAAAAGACCTTTGTTTTGTGTATTTTTATCGAATTAAAGAAAAAAAAAATAAACACCCGACGAGTTTTTTTTATAAATGGGCCGGCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATc                                                                                                                     <  2:542805/141‑1 (MQ=255)                                    caACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 >  1:391365/1‑151 (MQ=255)                                     acccggccttttttattCAAAAAAAAACCCTTCTTTTGCTTATTTTTATCGAAATAAAGAAAAAAAAAATAACCACCCGTCAGGATTTTGTTATAAATGAAACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                                      <  1:413317/133‑1 (MQ=255)                                       ctggCCGTTTTTTTTCTTAAAAAAACCTTTCTTTTGCGTACGTTTATAGAAATTAAGAAAAAAAAAATAAACACCGGTAAGGATATTTTTAAAAATGTGACGACCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                               <  1:164726/147‑1 (MQ=255)                                       cctccctttattttTCGTAAAAAAACCCTCTTTTTGCTTACGTTTAGCGAAATAAAGAAAAAAAAAATAAACACCCGTAAGGATTTTTTAAAAAGTGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                               <  2:70488/142‑1 (MQ=255)                                         ggcctttattatTCGTAAAAAAACCCTTCTTTTGCGTTCGTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGTAAGAATATTGTTATAAATGTAACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                            <  1:183133/145‑1 (MQ=255)                                            cGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                         >  2:183924/1‑151 (MQ=255)                                                 ttattaaaaaaaaGACCTTTCTTTGGCTTATTTTTATCGAAATAAAGAAAAAAAAAAACAACAACATTCAGAATATTGTTATAAATGTGACGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                                                                                      <  1:681851/139‑1 (MQ=255)                                                 ttatgCGTATAAAAACCCTTGTTTTGCGTACGTTTATCGAAATAAAAGAAGAAAAAATTAACACCCGTCAGGATTTTGTTAAAAATGGGACGACCCACAACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                                     <  1:500255/150‑1 (MQ=255)                                                         aaaaaaacccTTCTTTTGCGTCCGTTTATCGAAATAAAGAAAAAAAAAATAAACCCCCTCCAGAATATTTTTACCAATGTGACACCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                                                                                    <  2:448185/135‑1 (MQ=255)                                                         aaaaaGACCCTTCTTTTGCGTTCGTTTATCGAAATAAAGGAAGAAAAAATAAACACCCGTCAGTATTTTGTTATCAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                            <  1:58817/149‑1 (MQ=255)                                                           aaaaacctttCGTGTGCGTCGGTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGTAAGGATATTTTTATAATTGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                          <  2:566679/143‑1 (MQ=255)                                                                cctttCTTTTGCGTCCGTTTATGGAAATAAAGGAAAAAAAAATAAACACCCGTCAGGATATTGTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                      <  2:491561/147‑1 (MQ=255)                                                                 cctccttttGCGTCTTTTTATCGAAATTAAGAAAAAAAAAATAAACACCCGTAAGGATATTTTTAACAATGTGACACCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                      <  1:404277/143‑1 (MQ=255)                                                                   ttggTTTGCGTACTTTGAAGAGAATAAAGAAAAAAAAAATAAACAGCCGTGAGGATTTTGTTAAAAGTGAGACGAACCACTGCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                                   <  1:640735/147‑1 (MQ=255)                                                                          gCGTATTTTTATCGAAATTAAGAAAAAAAAAATAAACACCCGAAAGGATTTTGTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                             <  1:646614/149‑1 (MQ=255)                                                                             tACGTTTATCGAAATAAAGGAAAAAAAAATAAACACCCGTCAGAATTTTTTTATCAATGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                              <  2:359933/145‑1 (MQ=255)                                                                                ttttATCGAAAATAAGAAAAAAAAAATCAACACCCGGCAGGATTTTTTTATCAATGTGACCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          <  1:381230/145‑1 (MQ=255)                                                                                  ttatgaaaaTAAAGAAAAAAAAAATAAACACCCGTCAGGATATTTTTATAAATGTAACCACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg                                                                        <  1:116808/139‑1 (MQ=255)                                                                                     tCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                   >  2:339950/1‑148 (MQ=255)                                                                                      cgaaaTAAAGAAAAAAAAAATAAACACCCGTCAGGATATTTTTATCAATGTGACGCTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga                                                                <  1:120051/147‑1 (MQ=255)                                                                                                 aagaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                                    >  2:395575/1‑151 (MQ=255)                                                                                                            aacaccCGCCAGGATTTTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                         <  2:373455/151‑1 (MQ=255)                                                                                                            aacaacCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                         >  2:196604/1‑151 (MQ=255)                                                                                                             acaccCGTAAGGATTTTTTTATAAATGTCACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGa                                        <  2:434612/151‑1 (MQ=255)                                                                                                               aacCGTCAGGATATTGTTAAAATGGTGACCACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                      <  1:702478/151‑1 (MQ=255)                                                                                                                 cctacAGGATTTTTTTATAAATGTGACACCCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaaa                                    <  1:32316/147‑1 (MQ=255)                                                                                                                 ccGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaa                                     >  1:260234/1‑150 (MQ=255)                                                                                                                  cgtagggATTTTGTTAAAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGTTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAt                                   <  1:137872/146‑1 (MQ=255)                                                                                                                    tCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt                                  >  2:577860/1‑150 (MQ=255)                                                                                                                       ggATATTTTAAACAATGGGCCCATCCCCACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt                                  <  1:652658/147‑1 (MQ=255)                                                                                                                          tatTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAACACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                             >  1:168063/1‑149 (MQ=255)                                                                                                                            tttttATAAATGTAACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAg                           <  2:707281/146‑1 (MQ=255)                                                                                                                                    aaTGTGACGACCAACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg                     <  2:61260/147‑1 (MQ=255)                                                                                                                                     aTGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg                 >  2:147461/1‑150 (MQ=255)                                                                                                                                      tgagccgccccACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg                 <  2:344488/140‑1 (MQ=255)                                                                                                                                       gtgACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGata              <  2:237653/151‑1 (MQ=255)                                                                                                                                           cGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATc            <  1:536055/149‑1 (MQ=255)                                                                                                                                                    aCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  2:530571/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                      ATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421296‑2421592

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