Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,443	0	A	C	43.1%	82.5 / 11.8	51	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/22);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.47e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGT  >  NC_000913/2421296‑2421588                                                                                                                                                        |                                                                                                                                                  aTCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTAccc                                                                                                                                                <  1:119449/151‑1 (MQ=255)   ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCccg                                                                                                                                              <  2:344373/151‑1 (MQ=255)            aacaacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATaa                                                                                                                                       >  2:771846/1‑149 (MQ=255)                             aaggccacaCCCCGCCGTTTTTTTTCGTAAAAAGACCT‑TTCTTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATA‑AA‑CAACCGTCAGG‑ATTTTTTTATAAAGGTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAg                                                                                                                      <  1:159342/144‑1 (MQ=255)                                  cacaACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                                >  2:429032/1‑150 (MQ=255)                                  cacaACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCa                                                                                                                >  2:677378/1‑150 (MQ=255)                                      cccggCCGTTATTATTATAAAAAAAACCC‑TTCTTTTGCGTCCTTTTAGCGAAATAAAGAAAAAAAAAATA‑AA‑CACCCATCAGG‑ATATTTTTATCAATGTGCCGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                                                                              <  2:334434/144‑1 (MQ=255)                                       cgcccctcTATTATT‑GAAAAAAAAACCTTTGGTTTGGGTACGTTTAGGAGAATAAAGAAAAAAAAAAAA‑CA‑CAACCGGCAGG‑ATTTTTTTACAAAGGTGACGCCCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCaa                                                                                                               <  1:739278/139‑1 (MQ=255)                                       ccTGCCGTTATTATTCGTAAAAAGACCT‑TTCTTTTGCGTACGTTTA‑CGGAAATAAAGAAAAAAAAAAT‑AACCACCCGTAAGG‑TTTTTTTAATAAAGGGGACGTCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                                  <  1:201772/143‑1 (MQ=255)                                        cTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                         >  1:198613/1‑151 (MQ=255)                                        cTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                         >  2:412050/1‑151 (MQ=255)                                          gCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcgc                                                                                                       >  1:79291/1‑151 (MQ=255)                                             gCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                     >  2:548184/1‑150 (MQ=255)                                              cTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGa                                                                                                    >  1:527621/1‑150 (MQ=255)                                                  ttattaaaaaaaaaaCCTTTCGTTTGCGTATTTTTACCGAAATAAAGGAAAAAAAAATA‑AA‑C‑ACCGATAGGAATTTTTTTATAAATGTGACGCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGa                                                                                                    <  2:508873/138‑1 (MQ=255)                                                  tattCGAAAAAAAACCC‑TTCTTTTGCGTACTTTTACGGAAATAAAGAAAAAAAAAATA‑AA‑CAACCGTAAGA‑ATATTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                                                                                <  1:831032/150‑1 (MQ=255)                                                    ttCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                             >  2:501405/1‑151 (MQ=255)                                                      cGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                            >  1:619936/1‑150 (MQ=255)                                                         ataAAGAACC‑TTCGTTTGCGGACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAcc                                                                                         >  2:41313/1‑150 (MQ=255)                                                          aaaaaaaccttttttttGCGTCCTTTTATAGAAATAAAGAAAAAAAAAATA‑AA‑AACCCGTAAGA‑ATTTTGTAAAAAAGGTGCCTCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                                                           <  1:853867/134‑1 (MQ=255)                                                           aaaGAACT‑TTGTTTTGGTTCCTTTTATGGAAATAAAGGAAAAAAAAATA‑AC‑CACCCGTGAGT‑ATATTGTTAATAAGGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                      <  1:85230/151‑1 (MQ=255)                                                                      cGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                             >  2:780198/1‑150 (MQ=255)                                                                       ttttGCGTACTTTTATCGGAATAAAGGAAGAAAAAAAA‑AC‑CACCCGTAAGT‑ATATTTTAATCAATGTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACa                                                                            <  1:509585/149‑1 (MQ=255)                                                                        tttGCGTCCTTTTAGCGAAATAAAGAAAAAAAAAATT‑AA‑CACCCTTCAGA‑ATATTGTTATAAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACa                                                                            <  2:73151/149‑1 (MQ=255)                                                                             gTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                      >  2:248121/1‑150 (MQ=255)                                                                                 ttttATCGAAATAAAGAAAGAAAAAATC‑AA‑CACCCTTCAGT‑ATTTTGTTATAAATGTGACCCTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCa                                                                     <  2:435035/146‑1 (MQ=255)                                                                                  tttagcaaaaTAAAGGAAAAAAAAATA‑AA‑CACCCGTCAGG‑ATTTTTTTATAAATGTGACCACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                   <  2:324720/140‑1 (MQ=255)                                                                                  tttacggcaaTAAAGAAAGAAAAATTC‑AA‑CCACCCCAAGG‑ATTTTGTAATGAAGTTGACCACCACCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCa                                                                     <  1:607037/138‑1 (MQ=255)                                                                                   tttatggAAATAAAGAAAAAAAAATTAAA‑CACCCGTAAGG‑ATTTTTTTATAATTGTAACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatta                                                                         <  2:852196/137‑1 (MQ=255)                                                                                        ggAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATACACTACCCCACCGATAAAGCGAGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                             >  2:618480/1‑148 (MQ=255)                                                                                          aaTTAAGAAAGAAAAAATA‑AA‑AACCCGAAAGG‑ATTTTTTTAAAAATGTGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAAt                                                        <  1:455829/151‑1 (MQ=255)                                                                                          aaTTAAGAAAAAAAAAATA‑AA‑CACCCGTGAGA‑ATTTTTTTATAAATGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag                                                           <  1:756772/148‑1 (MQ=255)                                                                                          aaTAAAGAAAAAAAAAATA‑AA‑CACCCGTCAGG‑ATTTTTTTATCAATGTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                          <  2:694653/149‑1 (MQ=255)                                                                                               aGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCg                                                   >  2:733405/1‑151 (MQ=255)                                                                                                  aaaaagaaATA‑AA‑CCACCGCCGGA‑ATTTTATTAAAATTGTCCCGACCCCCACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGt                                                  <  1:384375/146‑1 (MQ=255)                                                                                                   agaagaaATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                                >  2:197314/1‑150 (MQ=255)                                                                                                      agaaATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                            >  1:567412/1‑151 (MQ=255)                                                                                                        aaaTT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                          >  1:868682/1‑151 (MQ=255)                                                                                                                  cacccgaaaggtttttGTAATCAATGTGCCAACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                     <  2:276045/135‑1 (MQ=255)                                                                                                                      gTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt                              >  2:542605/1‑151 (MQ=255)                                                                                                                       tCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt                              >  2:594883/1‑150 (MQ=255)                                                                                                                             atatTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                         >  1:327380/1‑150 (MQ=255)                                                                                                                              tatTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACt                          >  2:310162/1‑148 (MQ=255)                                                                                                                                tttttATAAATGTGACGACCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCa                        <  1:475582/145‑1 (MQ=255)                                                                                                                                ttttaataaaTGTGACGATCCACACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTaa                            <  2:471738/136‑1 (MQ=255)                                                                                                                                       cAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTc              >  2:6450/1‑151 (MQ=255)                                                                                                                                        aagttgACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgt                  <  1:298475/142‑1 (MQ=255)                                                                                                                                           gtgcccacccACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGa            <  1:694433/141‑1 (MQ=255)                                                                                                                                            tgACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGatat         >  2:496321/1‑151 (MQ=255)                                                                                                                                               cGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc      >  2:508508/1‑151 (MQ=255)                                                                                                                                                   ccACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGt  >  1:350126/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                  ATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACC‑TTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AA‑CAACCGTCAGG‑ATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGT  >  NC_000913/2421296‑2421588

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