Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,421,4020TA41.5% 62.2 / 17.2 41I26I (ATT→ATA) folXD‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base T (20/4);  new base A (0/17);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 4.70e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATAGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTT  >  NC_000913/2421259‑2421548
                                                                                                                                                |                                                                                                                                                  
ggATAGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAAcc                                                                                                                                             >  2:599401/1‑151 (MQ=255)
 gATAGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCg                                                                                                                                            <  2:210818/151‑1 (MQ=255)
 gATAGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCg                                                                                                                                            >  2:681901/1‑151 (MQ=255)
   tAGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTc                                                                                                                                          >  2:365759/1‑151 (MQ=255)
    aGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTc                                                                                                                                          >  2:584464/1‑150 (MQ=255)
       tAACAGAAACAGGTTCTCGTGTGTTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAgg                                                                                                                                       >  1:99267/1‑150 (MQ=255)
           aGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATAtt                                                                                                                                  <  2:556503/151‑1 (MQ=255)
           aGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATAtt                                                                                                                                  >  1:118101/1‑151 (MQ=255)
                      ctcGTGTATTATTTTATCCTAAGTAAAACAACGGAAAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAAAACCTTCGTTTGCGTACGTTTATCGGAATT‑TAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCaa                                                                                                                          >  1:547847/1‑148 (MQ=255)
                        cGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAAtgt                                                                                                                       <  1:512314/149‑1 (MQ=255)
                                 tttCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACt                                                                                                            >  1:5084/1‑151 (MQ=255)
                                     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                                                                     >  1:98467/1‑151 (MQ=255)
                                                                        caACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                     >  2:31031/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)
                                                                                                                                                |                                                                                                                                                  
GGATAGTTAACAGAAACAGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATT‑AAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTT  >  NC_000913/2421259‑2421548

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