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.7%	59.7 / 12.7	44	Y40S (TAC→TCC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base A (21/2);  new base C (0/21);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.26e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTG  >  NC_000913/2421301‑2421589                                                                                                                                                 |                                                                                                                                                   aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTCCCccgcc                                                                                                                                             <  2:166627/150‑1 (MQ=255)          aacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgc                                                                                                                                   >  1:38109/1‑151 (MQ=255)             ggAAAACCTCCAATGGCACAACCTGCCGTTATTATTCGTAAAAAGACCCTTCTTTTGCGTACTTTTATCGAAATA‑AAGAAAAAAAAAATA‑AACAACCGTCAGGATTTTTTTAAAATTGTGCCGACCCACTCCCCCGCCGATAAAGCGCGCa                                                                                                                                <  1:260097/151‑1 (MQ=255)                 acccTGCAAGGCCACACCCTCCCTTTATTATTCGTAAAAAGACCTTTCGTTTGCGTACGTTTATCGAAATA‑AAGAAAAAAAAAATA‑AACACCCGAAAGGATTTTTTTATAATTGTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTa                                                                                                                             <  1:776733/148‑1 (MQ=255)                      gCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                         >  2:817121/1‑149 (MQ=255)                        aaggcaccccccTCCCGTTATTATTCAAAAAAACACCCTTCTTTTGCGTATTTTTACCGAAATA‑AAGGAAAAAAAAAGA‑ACCACCCGAAAGCAATTTGTTAAAAATGTGCCCACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGaa                                                                                                                        <  2:70288/138‑1 (MQ=255)                                accttGCGGCTATTATTGTAATAAAAATCCTTCTTTTGGTTACTTTTATGGAAATT‑AAGAAAAAAAAAATA‑AACACCCGACAGGTTATTTTTACAAATGAGACCCCCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATc                                                                                                                  <  1:385654/142‑1 (MQ=255)                                    tGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                         >  2:411790/1‑151 (MQ=255)                                            ttattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                                  >  2:317561/1‑150 (MQ=255)                                                   tataAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAcc                                                                                          >  1:866569/1‑151 (MQ=255)                                                      aaaaacccTTCTTTTGCGTATTTTTATCGAAATA‑AAGAAAAAAAAAATA‑AACACCCGTCAGGATTTTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                       <  2:211574/145‑1 (MQ=255)                                                        aGAACCTTCTTTTGCGTACGTTTACCAAAATT‑AAAGAAGAAAAAATA‑AACCACCTTCAGGATATGGTTAGGAGGGTGGCGATCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                                                                                      <  1:317249/150‑1 (MQ=255)                                                             cTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                                 >  1:521999/1‑150 (MQ=255)                                                               tCGTTTGCGTTCTTTTATCGAAATA‑AAGAAAAAAAAAATA‑AACACCCGTAAGGATATTGTTATAAATGTGACAACCCCCACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                                <  2:203241/149‑1 (MQ=255)                                                                    tttcgtcattttatggAATTAAAGGAAAAAAAAATT‑AACCCCCTTCGGGATTTTTTAACCAATGTCACCACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                           <  1:817198/136‑1 (MQ=255)                                                                      cgtagttttagcaaaaTA‑AAGAAAAAAAAAATA‑AAAACCCGAGAGTATTTTTTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                           <  1:681850/133‑1 (MQ=255)                                                                        tACGTTTATCGGAATT‑AAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                       >  1:91724/1‑149 (MQ=255)                                                                         ccGTTTAT‑GGAAAT‑AAAGAAAAAAAAATTAAACACCCGTAAGGATTTTTTTAAAATTGTAACACCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                    <  1:374266/150‑1 (MQ=255)                                                                         acttttATGGAAATA‑AAGAAAAAAAAAAAA‑AAAACCCGACAGATTATTTTTATCAAGGTGACCACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCa                                                                      <  1:516179/146‑1 (MQ=255)                                                                              tagcgaaaTA‑AAGGAAGAAAAAATT‑AACACCCGTAAGAATTTTTTTACCAATGTGCCGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                <  1:362555/144‑1 (MQ=255)                                                                                   gAATT‑AAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                           >  1:8430/1‑150 (MQ=255)                                                                                    aaTT‑AAGAAAGAAAAATTT‑AACACCCTTCAGGTTATTGTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag                                                            <  1:724918/148‑1 (MQ=255)                                                                                         aaGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                     >  2:671985/1‑151 (MQ=255)                                                                                         aaGAAAAAAAAAATT‑AAAACCCGTCAGTATTTTTTTAACAAGGTGACTACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                      <  1:655121/150‑1 (MQ=255)                                                                                            gaagaagaaATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGtt                                                  >  2:399287/1‑151 (MQ=255)                                                                                             aaaaaaaaaTA‑AACACCCGTCAGGATATTTTTATAAATGGGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                                 <  1:324973/145‑1 (MQ=255)                                                                                             aaaaaaaaaTA‑AACACCCGTAAGATTTTTTTTATAAATGTAACCACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCg                                                    <  2:201347/142‑1 (MQ=255)                                                                                              agaagaaATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtc                                                >  2:387783/1‑151 (MQ=255)                                                                                                 aaaaaTT‑AACCCCCGTAAGTTTTTTTTTATAAATGGGCCGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                             <  2:219264/149‑1 (MQ=255)                                                                                                  gaaATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                            >  1:555717/1‑151 (MQ=255)                                                                                                    aaTT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                          >  2:358038/1‑151 (MQ=255)                                                                                                     aTT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGc                                         >  2:359230/1‑151 (MQ=255)                                                                                                         accaccggTCAGTATTTTTTTATCAATGTGACAACCCACACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                      <  2:869373/144‑1 (MQ=255)                                                                                                            aacctaaagtatatTGTTTTCAATGTAACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                            <  2:537863/132‑1 (MQ=255)                                                                                                            aacCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaa                                    >  1:106266/1‑150 (MQ=255)                                                                                                             acCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                   >  2:205267/1‑150 (MQ=255)                                                                                                              ccGTCAGTATATTGTTAACAATGTGACGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                   <  2:327808/149‑1 (MQ=255)                                                                                                                 gaagtattttGTAATAAAGGTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTa                              <  2:142197/143‑1 (MQ=255)                                                                                                                          tGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAt                     >  1:346836/1‑151 (MQ=255)                                                                                                                            ttAACAATGTGACAATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtg                    <  2:761144/150‑1 (MQ=255)                                                                                                                                        cGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc       >  1:694908/1‑151 (MQ=255)                                                                                                                                         gATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc       >  2:668563/1‑150 (MQ=255)                                                                                                                                         gATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc       >  1:643334/1‑150 (MQ=255)                                                                                                                                             cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg  >  2:155170/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                   AAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTG  >  NC_000913/2421301‑2421589

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