Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,350,8530AC45.3% 81.9 / 13.9 54S725A (TCA→GCA) arcBaerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base A (26/3);  new base C (0/24);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 4.17e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC  >  NC_000913/3350705‑3351003
                                                                                                                                                       |                                                                                                                                                      
tttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcc                                                                                                                                                      >  1:189465/1‑151 (MQ=255)
       cccccGTTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACTTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgccg                                                                                                                                                <  1:916443/150‑1 (MQ=255)
         cacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCAc                                                                                                                                             <  2:868613/151‑1 (MQ=255)
              tttCAGCCTTTCTACTTCGTGCCGCCACTCTTCTTTCACCCCTCCAACCAATCCACCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGCTTTTGCAGATGGCTTACCCCCACTGCCCCCGCCGCACCt                                                                                                                                           <  1:42609/148‑1 (MQ=255)
                    cACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAtttt                                                                                                                                   >  1:169503/1‑150 (MQ=255)
                     aCTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTAt                                                                                                                                 >  1:50445/1‑151 (MQ=255)
                      cTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATg                                                                                                                                >  2:951689/1‑151 (MQ=255)
                       ttCTCCGTCGTGCCGCCCCTCTTCTTTCATCTCTACAACCCATCCACCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTAt                                                                                                                                 <  2:422777/149‑1 (MQ=255)
                            cgtcgtGCCGCCCCTCTTTTTTCATCTCTCCAACCCATCCACCGACTTTTTCTCCCCAGCCCGGAAGGTCAGGAGACTAATTTTGCTGACCCAGCTGTTGCAGATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                                           <  2:127801/150‑1 (MQ=255)
                               cgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTc                                                                                                                       >  2:620364/1‑151 (MQ=255)
                                gtGACCCCACTCTTCTTTCATCTCTTCAATCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTc                                                                                                                       <  2:76622/150‑1 (MQ=255)
                                    cccccctctTCTTTCATCCCTTAAACCCATCCCCCGACTTTTTCTCCCCAGGCCGGAAGGTCGGGAACCTGAATTTGCTCACCCAGTTTTTGCAGAGGGCTTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca                                                                                                                   <  1:90701/145‑1 (MQ=255)
                                    cGCCCCTCTTCTTTCATCTCTCCAACCCATCCCCCGCCGTTATCTCCCCAGCCCGGAAGGTCAGGAGACTGAATTTCCTGCCCCAGTTGTTGAAGATGGCTTACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGtccttcct                                                                                                                        <  2:632314/145‑1 (MQ=255)
                                      cccctctTTTTTCATCTCTTCAACCCATTCCCCGACGTTATCTCCCCAGCCCGAAAGGTCAGGAGACTGATTTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca                                                                                                                   <  2:184148/145‑1 (MQ=255)
                                       cACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACTCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc                                                                                                               >  1:933573/1‑151 (MQ=255)
                                        aCTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc                                                                                                              >  2:1028597/1‑151 (MQ=255)
                                            ttcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttt                                                                                                          >  1:343100/1‑151 (MQ=255)
                                              cttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                        >  1:589273/1‑151 (MQ=255)
                                               tttCTTCTCTAAAACCCATTCCCCGACTTTATCTCCCCAGGCCGAAAGGTCAGGAGCCTGATTTTCCTGCCCCATTTTTTGCAGGTGCCATACCCCCCTTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                         <  1:34310/149‑1 (MQ=255)
                                               ttccacctctCCAATCCATCCCCCGTCGTTTTCTCCCCAGGCCGGAGGGTCGGGAGACTAAATTTGCTCCCCCAGCTTTTGCAGATGGCGTACCCCCCCTGGCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                        <  1:309471/144‑1 (MQ=255)
                                                ttCATCTCTCCAACCCATCCCCCGACGTTATCTTCCCAGGCCGGAGGGTCAGGAGACTGATTTTGCTGCCCCATTTGTTGCAGATGGCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                         <  1:910047/148‑1 (MQ=255)
                                                   aTCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCt                                                                                                   >  2:435/1‑151 (MQ=255)
                                                     ctctTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                                                 >  1:254238/1‑151 (MQ=255)
                                                        tccAACCCATTCCCCGACTTTATCTCCCCAGGCGGGAAGGCCGGGAGACTAATTTTGCTGCCCCATTTTTTGCAGATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                                                 <  2:192250/146‑1 (MQ=255)
                                                             cccttcccccGTCGTTATCTCCCCAGCCCGGAGGGCCAGGAGACGGATTTTGCGGCCCCAGTTGTTGAAGATGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCa                                                                                          <  1:882442/142‑1 (MQ=255)
                                                              ccATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                           >  2:379301/1‑148 (MQ=255)
                                                                aTTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt                                                                                      >  1:446783/1‑151 (MQ=255)
                                                                  tCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                    >  1:782046/1‑151 (MQ=255)
                                                                   ccccGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                    <  2:203054/148‑1 (MQ=255)
                                                                   ccccGACGTTATCTCCCCAGGCCGGAAGGTCGGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGAGGGCGAAACCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                    <  2:589897/148‑1 (MQ=255)
                                                                   cACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                    >  1:570026/1‑150 (MQ=255)
                                                                   cACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                    >  2:319028/1‑150 (MQ=255)
                                                                    cccGACTTTCTCTCCCCAGCCGGAAGGGTCGGGAGACGAATTTTGCGCCCCCGGTTGTTGCAGAGGCCTAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAt                                                                                       <  1:982174/145‑1 (MQ=255)
                                                                     ccGACGTTATCTCCCCAGGCGGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGTTTTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACt                                                                                 <  2:199181/151‑1 (MQ=255)
                                                                          gTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                                            >  2:650930/1‑151 (MQ=255)
                                                                                 tCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg                                                                      >  2:1013841/1‑150 (MQ=255)
                                                                                  cccAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGAc                                                                    >  2:500394/1‑151 (MQ=255)
                                                                                                     gagaCTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtctt                                                  >  1:178557/1‑151 (MQ=255)
                                                                                                       gaCTGATTTCGCGGCCCCAGCTGTTCCAGAGGGCTTAACCCCCCTGCCCCCGCCGCTCCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttc                                                 <  2:859447/150‑1 (MQ=255)
                                                                                                            aaTTTGCTGCCCCAGCTTTTGCAGATGGCGTAACCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAac                                           <  2:237677/151‑1 (MQ=255)
                                                                                                              tttGCTGACCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAaca                                          <  2:159911/150‑1 (MQ=255)
                                                                                                                 gCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGc                                      >  1:843364/1‑151 (MQ=255)
                                                                                                                   tGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCt                                     >  1:88565/1‑150 (MQ=255)
                                                                                                                         aGCTGTTCCGGGGGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGt                              <  1:577446/151‑1 (MQ=255)
                                                                                                                              ttGCAGGTGGCAAACCCCCCTGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCgg                           <  1:929739/149‑1 (MQ=255)
                                                                                                                              ttGCAGATGGCATACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGt                          <  1:411447/150‑1 (MQ=255)
                                                                                                                                  aGATGGCTAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGa                        <  1:1028576/148‑1 (MQ=255)
                                                                                                                                        cGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg               >  1:66581/1‑151 (MQ=255)
                                                                                                                                         gTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCgg              >  2:212595/1‑151 (MQ=255)
                                                                                                                                         cttacccccACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCtt                 <  2:696320/144‑1 (MQ=255)
                                                                                                                                              ccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCt          <  1:275122/151‑1 (MQ=255)
                                                                                                                                                ccACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTAc        <  2:797817/151‑1 (MQ=255)
                                                                                                                                                    tGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGt    >  1:569984/1‑151 (MQ=255)
                                                                                                                                                    tGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGt    >  1:167332/1‑151 (MQ=255)
                                                                                                                                                       aCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTc  >  2:514347/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
TTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC  >  NC_000913/3350705‑3351003

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