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

CATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTG  >  NC_000913/2421295‑2421589                                                                                                                                                      |                                                                                                                                                   cATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTAcc                                                                                                                                                  >  2:413404/1‑151 (MQ=255) cATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTAc                                                                                                                                                   <  1:261601/150‑1 (MQ=255)   tCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTAcccc                                                                                                                                                >  2:976271/1‑151 (MQ=255)    ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCccg                                                                                                                                               <  2:883204/151‑1 (MQ=255)       aaGTAAAACAACGGAGACCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCccgccg                                                                                                                                            <  2:516037/151‑1 (MQ=255)          tAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGAt                                                                                                                                          >  2:77390/1‑150 (MQ=255)           aaaaCAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATa                                                                                                                                         >  1:1064171/1‑150 (MQ=255)                   ggAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAgcgcgc                                                                                                                                 >  2:306862/1‑150 (MQ=255)                        aCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGc                                                                                                                           >  2:23344/1‑151 (MQ=255)                         ccagCATTGCCACACCCTGCCGCTTTTTTTAGAAAAAAAACCCTTCGTTGGCGTACTTTTATAGAAATTAAGGAAGAAAAAATAACCCCCCGCGAGAATTTTGTTATAAAT‑GTGCCGAGCCCCCCCCCCGCCGATAAAGCGCGCACt                                                                                                                              <  1:212779/144‑1 (MQ=255)                         ccTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                          >  1:614769/1‑151 (MQ=255)                         ccTGAAAGGACACACCCTCCCTTTATTATTAGTAAAAAGACCCTTCTTTTGCTTATTTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGTAAGAATATTTTTATAAGG‑GTGACGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                          <  1:221443/151‑1 (MQ=255)                                   cacaACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                                 >  2:243942/1‑150 (MQ=255)                                             cgttattttTCGTAAAAAACCCCTTCGTTTGCTTACTTTTATCGAAATAAAGAAAAAAAAAATTAACACCCGTCAGGATTTTTTTATCAAT‑GTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                      <  2:793569/148‑1 (MQ=255)                                             cGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                      >  1:145894/1‑151 (MQ=255)                                                  ttattCGAATAAAGACCCTTTGTTTGCGTCCGTTTATGGAATTTAAGGAAAAAAAAATAAACACCCGTCAGTATTTTTTTATCATT‑GTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                                  <  2:515051/150‑1 (MQ=255)                                                     ttCTTAAAAAGACCCTTCTTTTGCGTACTTTTATAGAAATAAAGAAAGAAAAAATTAACACCCGTGGGGATTTTTTTATAAAT‑GTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                               <  2:602095/150‑1 (MQ=255)                                                      tCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                             >  2:501001/1‑151 (MQ=255)                                                          aaaaaGACCCTTCTTTTGCGTACGTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTCAGGATTTTTTTATCAAT‑GTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAcc                                                                                          <  2:236762/148‑1 (MQ=255)                                                            aaaGACCTTTCGTTTGCGTCTGTTTAGCGAAATAAAGAAAGAAAAAATAAACACCCGT‑AAGATTTTTTTAACAATGGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                       <  1:704147/151‑1 (MQ=255)                                                              aaacctttCTTTTGCGTACTTTTATGGAAATAAAGAAAAAAAAAATAACCACCCGTAGGAATTTTGTTATCAAT‑GTACCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                     <  2:324960/145‑1 (MQ=255)                                                                acccTTCTTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTCAGAATTTTTTTATAAAT‑GTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                    <  2:1040897/148‑1 (MQ=255)                                                                acccTTCGTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATTAACACCCGTCAGAATATTTTTAACAAT‑GTGACGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                    <  1:1025513/148‑1 (MQ=255)                                                                    ttCGTTTGCGTATTTTTATCGAATTAAAGAAAAAAAAAATAAACACCCGTCGGAATTTTGTTATAAAG‑GTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                                 <  2:202347/149‑1 (MQ=255)                                                                      cGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                              >  1:206412/1‑150 (MQ=255)                                                                        tttGCGTATTTTTATGGGAAAAAAGGAAAAAAAAATAAACACCCGACAGGATATTGTTATAAAT‑GTGACTCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                              <  1:1107254/148‑1 (MQ=255)                                                                          tGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                         >  2:162829/1‑151 (MQ=255)                                                                                  tttATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg                                                                 >  2:947444/1‑151 (MQ=255)                                                                                  tttATCGAAATAAAGAAAGAAAAAATAAACACCCGTAAGGATTTTTTTATCAAT‑GTGACTACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg                                                                     <  1:577102/147‑1 (MQ=255)                                                                                  ttatgcgaaaTAAAGAAAGAAAAAATAAACACCCGTCAGAATTTTGTTATAAAG‑TGGACCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                         <  1:281496/135‑1 (MQ=255)                                                                                      tCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                >  1:778596/1‑148 (MQ=255)                                                                                          aaTTAAGAAAGAAAAAATAAACACCCGTCAGAATTTTTTTATAAAT‑GTGCCAATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAAt                                                         <  2:1040625/151‑1 (MQ=255)                                                                                          aaTAAAGAAAAAAAAAATAAACCCCCGTAAGTATTTTTTTATAAAG‑GTGGCGCTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                          <  1:888427/150‑1 (MQ=255)                                                                                             tAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                      >  2:1124265/1‑151 (MQ=255)                                                                                              aaGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                     >  1:235092/1‑151 (MQ=255)                                                                                              aaGGAAAAAAAAATAAACACCCGTAAGGATATTGTTATCAAT‑GGGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                        <  2:973124/148‑1 (MQ=255)                                                                                              aaGAAAAAAGAAATTAACCACCGAGGGGATATTGTAGAGAAT‑GCAACAACCCAATCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                      <  2:228026/150‑1 (MQ=255)                                                                                                  aaaaaaaaaTAAACACCCGTCAGGATTTTTTTATAAAT‑GGGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                                 <  2:528794/145‑1 (MQ=255)                                                                                                     aagaaATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtctc                                              >  2:950556/1‑151 (MQ=255)                                                                                                          aTTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGc                                         >  2:922301/1‑151 (MQ=255)                                                                                                            tAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                       >  2:1041709/1‑151 (MQ=255)                                                                                                             aacaccCGTCAGTATTTTTTTATAATT‑GGGACCCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                      <  2:1083769/151‑1 (MQ=255)                                                                                                                acccctCAGTATTTTGTTATCTAT‑CTGCCGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCt                                        <  2:190483/141‑1 (MQ=255)                                                                                                                aacCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaa                                    >  2:1101015/1‑150 (MQ=255)                                                                                                                 acCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATCACCGTTTCTCTTTGCTGGaaa                                   >  1:925457/1‑150 (MQ=255)                                                                                                                  ccGTCAGAATATTTTTATAAAT‑GTGACGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGa                                     <  2:673291/147‑1 (MQ=255)                                                                                                                      aagaatttttttATCAAT‑GTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTaa                             <  1:77386/141‑1 (MQ=255)                                                                                                                          atatTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                          >  2:698599/1‑150 (MQ=255)                                                                                                                           tatTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                          >  1:391646/1‑149 (MQ=255)                                                                                                                             tttttATCAAT‑GTGCCACCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAgg                       <  1:158589/147‑1 (MQ=255)                                                                                                                             tttttATAAAT‑GTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAgg                       <  2:254718/147‑1 (MQ=255)                                                                                                                                     atgtttcaCGACCCTCTTCACCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAt                     <  1:703325/138‑1 (MQ=255)                                                                                                                                     aaT‑GTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg              <  1:998123/151‑1 (MQ=255)                                                                                                                                      aT‑GTGCAGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTc               <  2:1008753/149‑1 (MQ=255)                                                                                                                                              gATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc       >  2:611269/1‑150 (MQ=255)                                                                                                                                                 ccACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGt   >  1:235738/1‑151 (MQ=255)                                                                                                                                                  cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg  >  2:1134231/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTG  >  NC_000913/2421295‑2421589

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