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

TATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCA  >  NC_000913/2421287‑2421585                                                                                                                                                        |                                                                                                                                                     tattatTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGAt                                                                                                                                                      >  1:496128/1‑151 (MQ=255)  attattTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGAt                                                                                                                                                      <  1:221940/150‑1 (MQ=255)             cTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgc                                                                                                                                          >  2:751707/1‑151 (MQ=255)                     tacaacGGAGAACCTGCAATGGCAAAGACTGCCGCTATTATTCGTATAAAGAACCTTAGTTTGCGTACGTATATAGGCAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCTCGCGGATa                                                                                                                                     >  2:413209/2‑148 (MQ=255)                       caacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACACCAGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgc                                                                                                                                 >  2:618874/1‑150 (MQ=255)                               aaCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                          >  2:921690/1‑149 (MQ=255)                                  cTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                      >  2:406251/1‑150 (MQ=255)                                        ggccacaCCCTGCCGTTATTATTCATAAAAAAACCCTTCTTTTGCGTACTTTTATCGGAAT‑AAAGAAAAAAAAAATAAACACCCGTCAGGATTTTTTTATAAATGTGACTACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                               <  1:439800/148‑1 (MQ=255)                                              acccTGCCTTTATTTTTCAAAAAAAAACCTTTCTTTTGCGTACTTTTATCGAAAT‑AAAGAAAAAAAAAATAAACACCCGTCAGAATTTTGTTAAAAAGGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                             <  1:589389/145‑1 (MQ=255)                                               ccctcccttttttattCAAAAAAAAACCCTTCGTTTGCTTATGTTTAT‑GGAATAAAAGAAAGAAAAAATAAACACCCGCAAGGATTTTTTTAAAAAGGAGCCGCCCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATc                                                                                                              <  2:1040046/134‑1 (MQ=255)                                                ccTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                        >  2:464979/1‑150 (MQ=255)                                                           tattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGATGATCCACTATCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                                                                             >  2:136193/1‑150 (MQ=255)                                                           tattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                                                            >  1:904030/1‑151 (MQ=255)                                                           tattCGAAAAAAGACCCTTCTTTTGCGTACTTTTATGGAAAT‑AAAGAAAAAAAAAATAAACACCCGTCAGAATTTTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                                                                             <  1:733109/150‑1 (MQ=255)                                                             ttctaaaaaaaaaCCTTTTTTGGCTTACTTTTATAGAAAT‑AAAAAAAAAAAAAATAACCAACCTTAAGGATTTTTTAATAAAGGTAACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                              <  2:315775/140‑1 (MQ=255)                                                                  ataAAGAACCTTCGTTTGCGTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCAGATAAAGAGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                     >  2:475222/1‑151 (MQ=255)                                                                  aaaaaaacccTTCTTTGGCGTACTTTTATCGAAAT‑AAAGGAAAAAAAAATAAACACCCGTCAGAATTTTGTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                          <  1:31967/138‑1 (MQ=255)                                                                   taAAGAACCTTCGTTTGCGTTCGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGAACCGt                                                                                    >  2:719076/1‑151 (MQ=255)                                                                        acccTTCTTTTGCGTACGTTTATCGAAAT‑AAAGAAAAAAAAAATTAACACCCGTCAGAATTTTGTTATCAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                <  1:516358/148‑1 (MQ=255)                                                                        acccTTCGTTTGCGTACTTTTATCGAAAT‑TAAGGAAAAAAAAATTAACACCCGTAAGGATTTTTTTATCAATGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                <  1:872293/148‑1 (MQ=255)                                                                                tttGCTTGCTTTTATGAGAAT‑AAAGAAATAAAAAATAAACACCCGACAGGATTTTGTTAAAAAGGTGCCGCTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACat                                                                        <  1:783925/150‑1 (MQ=255)                                                                                ttgtcttgtttttATAGAAAT‑AAAGGAAAAAAAAATTAACACCTGTAGGGTTTTTTTTATAAATGTGACTACCCACCACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                             <  2:1042318/135‑1 (MQ=255)                                                                                     gTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                   >  2:513522/1‑150 (MQ=255)                                                                                      tACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAACAACATTATTc                                                                   >  2:44838/1‑149 (MQ=255)                                                                                          tttATCGGAAT‑GAAGGAAGAAAAAATAAACAGCCGTGGGGAGTTTTTTATCATGTTGTCAACCCTCTCCCCCGCCGATAGAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                              <  2:184335/150‑1 (MQ=255)                                                                                            tATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                            >  1:844741/1‑150 (MQ=255)                                                                                            tATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                            >  2:865855/1‑150 (MQ=255)                                                                                             aTCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                            >  2:528232/1‑149 (MQ=255)                                                                                              tCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                            >  2:521491/1‑148 (MQ=255)                                                                                                 gAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                       >  2:88315/1‑150 (MQ=255)                                                                                                  aaT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAAt                                                     >  1:50688/1‑151 (MQ=255)                                                                                                       aaGAAAAAAAAAATAAACACCCGTAGGTATTTTGTTATAAATGTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                    <  1:681772/148‑1 (MQ=255)                                                                                                        atgacgaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCAGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACTGTAACGAAAAACATTATTCATCATGTAGAGAATAACCg                                                >  2:688542/6‑151 (MQ=255)                                                                                                           aaaaaaaaaTAAACACCCGTCAGGATATTGTTATCAATGGGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCg                                                <  1:837843/142‑1 (MQ=255)                                                                                                           aaaaaaaaaTAAACAACCGTCAGGATTTTTTTATAAAGGTGACGACCCCCACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGt                                               <  1:722810/143‑1 (MQ=255)                                                                                                              aaaaaaTAACCAACCGACAGTATTTTTTTATCAATGGGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                           <  2:442181/147‑1 (MQ=255)                                                                                                               aaaaaTAAACACCCGTCAGAATTTTGTTATCAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                         <  1:95026/149‑1 (MQ=255)                                                                                                                 aaaTAAACACCCGTAAGGATTTTTTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                       <  1:70012/151‑1 (MQ=255)                                                                                                                  aaTAAACACCCGACAGAATATTTTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                      <  1:154392/151‑1 (MQ=255)                                                                                                                      aacaacCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                  >  1:390631/1‑151 (MQ=255)                                                                                                                       acaacCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                  >  1:415958/1‑150 (MQ=255)                                                                                                                         acccGTCAGAATTTTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGa                                 <  1:486615/147‑1 (MQ=255)                                                                                                                          acCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGACGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAACGTTTCTCTTTGCTGTaa                                >  2:134132/1‑149 (MQ=255)                                                                                                                          acCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                               >  2:291465/1‑150 (MQ=255)                                                                                                                            cGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACTTTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAt                            >  2:501742/1‑151 (MQ=255)                                                                                                                             gTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt                           >  1:169249/1‑151 (MQ=255)                                                                                                                                    tatTTTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTaa                         <  1:467417/146‑1 (MQ=255)                                                                                                                                       tGTTATCAATGTGACGATCCACTACCCCTCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAgg                   >  2:393375/1‑149 (MQ=255)                                                                                                                                       tGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAt                 >  2:207069/1‑151 (MQ=255)                                                                                                                                        gtaaaaaaTTTGACGACCCCCACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGa                  <  1:761622/143‑1 (MQ=255)                                                                                                                                           ataaaggtgACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgt               <  1:381474/143‑1 (MQ=255)                                                                                                                                              aaTGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg          >  1:714717/1‑151 (MQ=255)                                                                                                                                               aTGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg          >  2:1017060/1‑150 (MQ=255)                                                                                                                                                     cGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc   >  1:579641/1‑151 (MQ=255)                                                                                                                                                       aTCCACTACCCCGCCGATAAAGCGCGCAATAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGAAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCa  >  2:365948/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                     TATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAAT‑TAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCA  >  NC_000913/2421287‑2421585

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