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

TAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCAC  >  NC_000913/2421300‑2421586                                                                                                                                                  |                                                                                                                                                tAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCccgc                                                                                                                                           >  2:169803/1‑150 (MQ=255)  aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCccgccg                                                                                                                                         >  2:79868/1‑151 (MQ=255)      aaaaCAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATa                                                                                                                                      >  2:40746/1‑150 (MQ=255)       aaaCAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTTGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATaa                                                                                                                                     >  1:8889/1‑150 (MQ=255)        aacaacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATaa                                                                                                                                     >  1:397621/1‑149 (MQ=255)            acGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAgcgcg                                                                                                                               >  1:527952/1‑151 (MQ=255)              ggAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAgcgcgc                                                                                                                              >  2:1030502/1‑150 (MQ=255)              ggAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAgcgcgc                                                                                                                              >  2:623401/1‑150 (MQ=255)              ggAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAgcgcgc                                                                                                                              >  2:496242/1‑150 (MQ=255)                         aagggCACACCCTGCCGTTATTATTCTTAAAAAGACCTTTCGTTTGCGTTCTTTTATCGAAATTAAGG‑AAGAAAAAATAAACACCCGTAAGAATATTGTTATCATTGTGACAA‑TCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAg                                                                                                                    <  2:283250/146‑1 (MQ=255)                                        cgttattttTCTTAAAAAAACCTTTCTTTTGCTTCCTTTTATCGGAATAAAGA‑AAAAAAAAATAAACACCCGTCAGTATTTTTTTATAAATGGGCCAC‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                   <  2:288940/148‑1 (MQ=255)                                            attattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                               >  2:225028/1‑151 (MQ=255)                                                     aaaaaaacccTTCTTTTGCGTACTTTTATCGAAATAAAGA‑AAAAAAAAATAAACACCCGTAAGAATTTTTTAACAAATGTGACGA‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAcc                                                                                       <  2:1115057/142‑1 (MQ=255)                                                       aaaaaccttttgTTTGCTTGCTTTTAGCGAAATAAAGA‑AAAAAAAAATAAACACCCGCCAGGATTTTTTAATCAATGTAACGA‑CCACCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                                                                        <  1:57025/136‑1 (MQ=255)                                                         aGACCCTTCTTTGGCGTACTTTTATCGAAATAAAGA‑AAAAAAAAATAAACACCCGTCAGGATTTTTTTACAAATGTGACGC‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                  <  2:505473/151‑1 (MQ=255)                                                               tctgTTTGGGTACTTTTAACGGAATTAAGA‑AAAAAAAAATCAACACCCTTAAGGATATTGTTCTAAATG‑GACCACCCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                              <  1:124659/146‑1 (MQ=255)                                                                   tttGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                        >  1:890364/1‑151 (MQ=255)                                                                      gCGTACTTTTATCGAAATAAAGA‑AAAAAAAAATAAACACCCGTAAGGATTTTTTTAACAATGTGACGA‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                       <  2:329107/149‑1 (MQ=255)                                                                        gttcGTTTATC‑GAATTAAGGAAAAAAAAAATAAACACCCGTCAGGATATTTTTATCAATGTGACGA‑CCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                    <  2:51652/147‑1 (MQ=255)                                                                         tACTTTTATGGAAATAAAGA‑AAAAAAAAATAAACAACCGTAAGAATATTTTAATCAATGGGCCCA‑TCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCa                                                                   <  1:705593/150‑1 (MQ=255)                                                                         tACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                    >  1:914371/1‑149 (MQ=255)                                                                           actttaagcgAAATAAGGAAAAAAAAAATAAACAACCGTGAGGTTATTGTTATAAATGAAACAC‑CCAACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattatt                                                                     <  2:636839/138‑1 (MQ=255)                                                                            ttttATAGAAATTAAGA‑AAAAAAAAATAAAAACCCATCAGTTTATTTTTATAAATGTGACTA‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCa                                                                <  1:722351/149‑1 (MQ=255)                                                                             tttATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg                                                              >  2:937398/1‑151 (MQ=255)                                                                             tttATCGAAATAAAGA‑AAAAAAAAATAACCACCCGTAAGGATTTTTTTATCAATGTGACGC‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg                                                                  <  1:825494/147‑1 (MQ=255)                                                                             tttAGCGAAATAAAGG‑AAAAAAAAATAAACACCCGACAGATTTTTGTTATAAATGTGCCGA‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                 <  1:1111340/148‑1 (MQ=255)                                                                              ttATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                             >  1:808491/1‑151 (MQ=255)                                                                                    gAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                        >  2:455312/1‑150 (MQ=255)                                                                                         aaGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                  >  1:729980/1‑151 (MQ=255)                                                                                         aaGA‑AAAAAAAAATAACCACCCGTAAGGATATTTTTATAAATGTGACAA‑CCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                  <  1:1033825/151‑1 (MQ=255)                                                                                               agaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtc                                             >  2:437866/1‑151 (MQ=255)                                                                                                 aaaaaaTAAACACCCGTAAGGTTTTTTTAATAAAGGTAACCA‑CCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtctc                                           <  1:628534/148‑1 (MQ=255)                                                                                                 aaaaaaTAAACACCCGTAAGGATTTTGTTATAAATGTGCCGA‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                            <  2:232806/147‑1 (MQ=255)                                                                                                  aaaaaTAAACACCCGCCAGAATATTGTTACCAATGTAACCA‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                          <  2:723575/149‑1 (MQ=255)                                                                                                   gaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                         >  2:1024678/1‑151 (MQ=255)                                                                                                     aaTTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                       >  2:15797/1‑151 (MQ=255)                                                                                                         aacaccGGTCAGTATATTGTTTTCAGTGAGACAT‑CCCACGCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                   <  2:25924/151‑1 (MQ=255)                                                                                                            acccGTAAGGTTATTTTTAAAATTGTGCCGA‑CCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaa                                 <  1:715072/148‑1 (MQ=255)                                                                                                            aacCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaa                                 >  1:606756/1‑150 (MQ=255)                                                                                                              ccGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaa                               >  1:460625/1‑150 (MQ=255)                                                                                                              ccGTAAGAATATTTTTATCAATGTGACGA‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                <  1:272175/149‑1 (MQ=255)                                                                                                                  aaGGATTTTGTTATCAATGTAAAGC‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTa                           <  1:344211/149‑1 (MQ=255)                                                                                                                   aGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAAc                         >  2:357528/1‑151 (MQ=255)                                                                                                                    ggATTTTGTTATCAATGTGCCGA‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTaa                          <  2:15128/149‑1 (MQ=255)                                                                                                                    ggATATTTTAATCAATGTGACCC‑TCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAt                             <  1:952728/146‑1 (MQ=255)                                                                                                                      atatTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACt                        >  2:310870/1‑149 (MQ=255)                                                                                                                      atatTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                       >  1:676847/1‑150 (MQ=255)                                                                                                                             tATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg               >  1:660053/1‑151 (MQ=255)                                                                                                                              ataaaTGTGACGA‑CCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg               <  1:339062/147‑1 (MQ=255)                                                                                                                                cAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTc            >  2:543661/1‑151 (MQ=255)                                                                                                                                    gtcccGA‑CCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgt                <  1:525438/139‑1 (MQ=255)                                                                                                                                          a‑TCCACTACCCCGCCGATAAAGCGCGCACTAGGGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCAc  >  1:227209/1‑151 (MQ=255)                                                                                                                                          a‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCAc  >  1:239242/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                TAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGG‑AAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGA‑TCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCAC  >  NC_000913/2421300‑2421586

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