Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,443	0	A	C	51.3%	43.2 / 57.7	39	Y40S (TAC→TCC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base A (17/2);  new base C (0/20);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 3.35e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-07

ATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421296‑2421592                                                                                                                                                     |                                                                                                                                                      aTCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAccc                                                                                                                                                    <  1:667645/151‑1 (MQ=255)    cTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgc                                                                                                                                                 >  1:708293/1‑151 (MQ=255)                       ccctaaaaggccacaCCCGGCCTTTATTATTCATAAAAAAACCCTTTGTTTGCGTACTTTTATCGAAATT‑AAGAAAAAAGAAATAAACACCCGTAAGGATATTGTAAAAAATGTGACGTCCCACCACCCCGCCGATAAAGCGCGCACTAg                                                                                                                               <  2:114649/139‑1 (MQ=255)                        ccTGCAAGGCCACACCCTGCCGTTTTTTTTCGAAAAAAAACCCTTCTTTTGCGTACTTTTATCGAAATA‑AAGAAAAAAAAAATAAACACCCGTCAGGATATTTTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGc                                                                                                                              <  2:109632/150‑1 (MQ=255)                        ccTACAATGCCACACCCCGCCTTTTTTATTCGAAAAAAGACCCTTCGTTTGCGTCTGTTTACGGGAATA‑AAGAAAGAAAAAATAAACACCCGGAAGGATTTTGTTATAAATGGGCCGACCCACTCCCCCGCCGATAAAGCGCGCACTAGc                                                                                                                              <  1:463700/150‑1 (MQ=255)                             aggggCACACCCTGCCGTTTTTATCCGTAAAAAGAACTTTCTTTGGGGTGTTTTTATGGAAATT‑AAGAAAAAAAAATTAAACACCCGCCAGGATATTTTAATCAATGGCACGACCCTCGCCCCCGCCGATAAAGCGCGCACTAGCGAAGa                                                                                                                         <  2:318229/147‑1 (MQ=255)                                    caACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 >  2:31031/1‑151 (MQ=255)                                    caACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 >  1:98467/1‑151 (MQ=255)                                     aaCCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 >  1:273052/1‑150 (MQ=255)                                     aaCCTGCCGATATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCAGTAGCGAAGATATCAACg                                                                                                                 >  1:452411/1‑150 (MQ=255)                                       cctccctttattatTTGTAAAAAAACCTTTTGTTTGCTTTC‑TTTTTTGGAATTAAAGAAAAAAAAAATAAACACCCGGAAGGTTTTTGTTATAATTGTGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                                                                                  <  1:405216/139‑1 (MQ=255)                                       ccTGCCGTTATTATTCGTAAAAAGACCCTTCGTTTGCGTACGTTTATCGAAATT‑AAGGAAGAAAAAATAAACACCCGTAAGGATATTTTTATCAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGa                                                                                                                <  2:161978/149‑1 (MQ=255)                                         tccctttattatTAGAAAAAAAACCCTTCGTTTGCTTCGGTTAATCGAAATA‑AAGAAAAAAAAAATAAACACCCGTAAGGATTTTTTAATAAAGGGGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                               <  1:513090/142‑1 (MQ=255)                                           cgttttttattCTTAAAAAAACCCTTTTTTGGCGTACTTTTAGCGAAATA‑AAGGAAAAAAAAATAAACACCCCCAAGGTTATTTTTACCAAGGGGACGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                             <  2:632877/142‑1 (MQ=255)                                              cTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGa                                                                                                        >  2:278029/1‑150 (MQ=255)                                                         aaaaagacctttgttttGCGTATTTTTAT‑GGAATTAAAGAAAAAAAAAATAAACACCCACAAGAATTTTTTAAAAAATGTGACGCCCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                                                                              <  2:286790/135‑1 (MQ=255)                                                           aaaGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                           >  2:302021/1‑150 (MQ=255)                                                            aaGAACCTTCGTTTGCGTATGTTTATCGGAATA‑AGGGAAGAAAAAATAAACAACCGTCGGGATATTTTTATAAAGGGGACGCTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                                                                                         <  1:361097/151‑1 (MQ=255)                                                                      gTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACat                                                                               >  1:458567/1‑151 (MQ=255)                                                                               cGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                       >  2:103774/1‑150 (MQ=255)                                                                                 tttATGGAAATA‑AAGAAAAAAAAAATAAACACCCGTCAGGATATTTTAATAAATGTGACAATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTAAg                                                                        <  2:315792/147‑1 (MQ=255)                                                                                 tttATCGGAATT‑AAGGAAGAAGAAATTAACAACCGACAGGATATTGTTATCAATGTGACGATCCAATAACCCGCCGATAAAGCGAGCACTAGCGAAGATATCACCGATGCGCTGAATTATCGCACAGTAACGAAAAACATTATTAAGCATg                                                                    >  1:666979/1‑151 (MQ=255)                                                                                       gaaaTA‑AAGAAAAAAAAAATAAACACCCGCAAGAATTTTTTTATCAATGTGCCGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                              <  2:115365/149‑1 (MQ=255)                                                                                        gAATT‑AAGGACGAAGAAATTAACCACCGTCAGGATATTGTTATCAATGTGACGATCCATTACCCCGCCGACAAGGCGCGCACTAGCGAAGACATCAACGATGCGCTGCATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag                                                               >  1:26429/1‑149 (MQ=255)                                                                                         aaTA‑AAGAAAAAAAAAATAAACACCCGTCAGTATTTTTTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga                                                                <  2:491869/147‑1 (MQ=255)                                                                                                  aagaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                                    >  2:198873/1‑151 (MQ=255)                                                                                                  aaaaaaaaaTAAACACCCGTCAGGATATTGTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                                    <  2:84004/145‑1 (MQ=255)                                                                                                  aaaaaaaaaTAAACACCCGTAAGAATTTTGTTATCAATATGACACTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                                    <  1:599109/145‑1 (MQ=255)                                                                                                        aaaTAAACACCCGTCAGAATTTTTTTATAAATGTGCCACCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                               <  2:548115/150‑1 (MQ=255)                                                                                                        aaaTAAACACCCGCCAGAATATTTTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                              <  2:99289/151‑1 (MQ=255)                                                                                                                        ggATATTTTTATAAATGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTaa                                <  2:55380/149‑1 (MQ=255)                                                                                                                          atatTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                             >  1:32487/1‑150 (MQ=255)                                                                                                                             tttttATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAg                           <  2:395494/146‑1 (MQ=255)                                                                                                                                ttAACAATGAGCCGCGCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtg                       <  2:708613/150‑1 (MQ=255)                                                                                                                                          gacactCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGa                <  1:1075/142‑1 (MQ=255)                                                                                                                                          gACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATc            >  1:238296/1‑151 (MQ=255)                                                                                                                                            cGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc          >  2:341484/1‑151 (MQ=255)                                                                                                                                              aTCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCAc        >  2:91752/1‑151 (MQ=255)                                                                                                                                                    tACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  1:330146/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421296‑2421592

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