Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,437	0	T	C	48.1%	50.4 / 22.4	53	I38T (ATC→ACC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base T (21/6);  new base C (0/25);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 1.61e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATC  >  NC_000913/2421290‑2421582                                                                                                                                                    |                                                                                                                                                  tattTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATcc                                                                                                                                                 >  2:116834/1‑150 (MQ=255)  attTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCa                                                                                                                                                <  1:314046/150‑1 (MQ=255)    ttCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACt                                                                                                                                              >  1:386293/1‑150 (MQ=255)    ttCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACt                                                                                                                                              >  1:237602/1‑150 (MQ=255)       aTCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAccc                                                                                                                                          >  2:430824/1‑151 (MQ=255)        tCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAcccc                                                                                                                                         >  2:80892/1‑151 (MQ=255)        tCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAcccc                                                                                                                                         >  1:68947/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)                                                                                                                                             gtgACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGata    <  2:237653/151‑1 (MQ=255)                                                                                                                                                 cGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATc  <  1:536055/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATC  >  NC_000913/2421290‑2421582

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