Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,350,8440AC45.6% 86.8 / 18.1 57L728V (TTA→GTA) arcBaerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base A (27/4);  new base C (0/26);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 3.11e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGG  >  NC_000913/3350701‑3350978
                                                                                                                                                |                                                                                                                                      
ggCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTGCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACt                                                                                                                                 <  2:851748/151‑1 (MQ=255)
          acccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgccg                                                                                                                       >  2:402693/1‑151 (MQ=255)
                 cTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCttt                                                                                                                >  1:597846/1‑151 (MQ=255)
                     cAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAtt                                                                                                            >  1:784430/1‑151 (MQ=255)
                         aCTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTAt                                                                                                        >  1:707433/1‑151 (MQ=255)
                         aCTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTAt                                                                                                        >  2:388290/1‑151 (MQ=255)
                           ttCTACGTCGTGCCCCCCCTCTTCTTTCATCTCTTCAACCCATCCCCCGCCTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGTTTTTGCAGATGCCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTAt                                                                                                        <  1:358961/149‑1 (MQ=255)
                               aCGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                  >  1:359814/1‑151 (MQ=255)
                                cgtcgtGACGCCACTCTTCTTTCATCTCTCAAACCCATCCCCCGACGTTATCTCCCCAGGCGGGAAGGTCAGGAGACTAAATTGGCTGCCCCAGTTGTTGAAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtccttc                                                                                                 <  2:204441/151‑1 (MQ=255)
                                 gtcgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcc                                                                                                >  1:814445/1‑151 (MQ=255)
                                 gtcgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcc                                                                                                >  1:814431/1‑151 (MQ=255)
                                  tcgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTCAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtccttcct                                                                                               >  2:321383/1‑151 (MQ=255)
                                    gtGACCCCCCTCTTTTTTCATCTCTTCAACCCATCCCCCGCCGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtccttcct                                                                                               <  2:187836/149‑1 (MQ=255)
                                     tGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTca                                                                                             >  1:208182/1‑150 (MQ=255)
                                        cccccctctTCTTCCATCTCTTAAACCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCGGCCCCATTTGTTGAAGATGGCATACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa                                                                                         <  1:717187/146‑1 (MQ=255)
                                         ccccctctTTTTTCATCTTTTCAATCCATTCCCCGACTTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGTTGTTGAAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAAt                                                                                        <  1:316948/147‑1 (MQ=255)
                                          cctctctTCTTTCATCTCTGCAATCCATCCCCCGAGTTTATCCCCCCAGGCCGGAAGGTCAGGAGATTAATTTTCCTACCCCATTTGTTGGAGATGGCATACCCCCACAGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAAt                                                                                        <  2:158235/147‑1 (MQ=255)
                                          cccttcttctTTCACCTCTACAACCAATCCCCCACCTTTATCTCCCCAGGCGGGAAGGTGAGGAGACTAATTTTGCTGACCTAGTTGTTGCAGAGGGCATACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                       <  1:366559/147‑1 (MQ=255)
                                          cccctctTCTTTAATCTCTCCAACCCATCCCCCGACTTTATCTCCCCGGGCGGGAAGGTCGGGAGACTAAATTTGCTGACCCGCCTGTTGCAGAGGGCGTACCCCCCCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTca                                                                                             <  2:380875/142‑1 (MQ=255)
                                                ttcttTCATCTCTTCAACCCATTCCCCGCCTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTAAATTTGCTCCCCAAGTTTTTGCAGATGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                  <  1:406783/150‑1 (MQ=255)
                                                         ctctCCAACCCATCCCCCGACTTTATCTCCCCAGGCCGGAGGGTCAGGAGACTAAATTTGCTGCCCCAGTTGTTGCAGATGCCGAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATcc                                                                           <  2:589672/148‑1 (MQ=255)
                                                           ctcaaacccATTCCCCGACGTTATCTCCCCGGGCCGGAAGGTCAGGAGACTCAATTTGCTCACCCAGCTGTTGCAGATGGCGTACCCCCCCTCACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                        <  1:321338/142‑1 (MQ=255)
                                                           ctTCACCCCATCCCCCGACTTTATCTCCCCAGCCGGGAAGTTCAGGAGACTGATTTTCCTCACCCTTCTGTTGCAGTGGGCTTAACCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTAt                                                                             <  1:455618/144‑1 (MQ=255)
                                                           ctTAAACCCATCCACCGACGTTTTCTTCCCAGGCGGAAAGGTCAGGAGCTTGATTTTGCTGACCCAGTTGTTGAATGTGGCGAACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCt                                                                          <  2:775446/147‑1 (MQ=255)
                                                               aacccATCCCCCGCCTTTTTCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                         <  2:329860/141‑1 (MQ=255)
                                                               aaTCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                  >  1:710202/1‑151 (MQ=255)
                                                               aaTCCATCCCCCGACGTTTTCTCCCCAGGCGGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                     <  2:749972/148‑1 (MQ=255)
                                                                  ccATCCCCCGACGTTATCTCCCCGGGCGGGAAGGTAAGGAGACTGATTTTGCTGCCCCAGCTTTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGt                                                                   <  2:176372/147‑1 (MQ=255)
                                                                    aTTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt                                                             >  1:593425/1‑151 (MQ=255)
                                                                    aTTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt                                                             >  2:704060/1‑151 (MQ=255)
                                                                       ccccgcctttATCTCCCCAGGCCGGAAGGTCAGGAGACTAAATTTGCTGCCCCAGTTTTTGCAGATGGCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGa                                                               <  2:177756/138‑1 (MQ=255)
                                                                       ccccGACTTTATCTCCCCAGGCGGGAAGGTCGGGAAACTAAATTCGCCACCCCAGCTGTTGCAGATGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGa                                                               <  2:452812/144‑1 (MQ=255)
                                                                       ccccGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTCACCCAGTTGTTGAAGATGGCGTAACCCCCCGCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                           <  1:47640/148‑1 (MQ=255)
                                                                        aCCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGAc                                                         >  2:526526/1‑151 (MQ=255)
                                                                         ccGGCGTTTTCTCCCCGGGCCGGAAGGTCAGAAAACTAATTTGGCTCCCCTAGCTTTTCAGGATGGTGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGAc                                                         <  2:637857/150‑1 (MQ=255)
                                                                           gACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                      >  1:179473/1‑151 (MQ=255)
                                                                           gACGTTATCTCCCCAGGCGGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                       <  2:483999/150‑1 (MQ=255)
                                                                              gTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                   >  2:145401/1‑151 (MQ=255)
                                                                                tATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                 >  2:480644/1‑151 (MQ=255)
                                                                                       ccAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGAc                                           >  1:523185/1‑150 (MQ=255)
                                                                                         aGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                         >  2:380832/1‑150 (MQ=255)
                                                                                         aGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                         >  2:79295/1‑150 (MQ=255)
                                                                                          gcgcGGAAGGCCAGGACACTAATTTGCTTCCCCAAGTTTTTGCAAGTGGCTTACCCCCCTCCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                         <  1:155335/146‑1 (MQ=255)
                                                                                            ccGGAAGGTCGGGAGACTGATTTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATa                                        <  2:623400/148‑1 (MQ=255)
                                                                                            ccGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACACTGACATAGcc                                     >  2:344160/1‑151 (MQ=255)
                                                                                                aaGGTCAGGAGTCTGAATTGGCTGCCCCATCTGTTGCAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCgg                                  <  1:371127/150‑1 (MQ=255)
                                                                                                 aGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                 >  2:224779/1‑150 (MQ=255)
                                                                                                   gTCAGGAGACGGATTTTCCTGCCCCTGTTTTTGCAGATGGCCTACCCCCCCGCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                 <  1:567921/148‑1 (MQ=255)
                                                                                                     cAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatcat                            >  2:87837/1‑151 (MQ=255)
                                                                                                        gagaCTGATTTTGCTGACCCAGCTGTTGCAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtctt                         <  2:870259/151‑1 (MQ=255)
                                                                                                              gAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaa                    >  2:549749/1‑150 (MQ=255)
                                                                                                               attttGCTGCCCCAGTTGTTGCAGATGGCATACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCa                     <  1:622702/146‑1 (MQ=255)
                                                                                                                   tGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg              >  2:590360/1‑151 (MQ=255)
                                                                                                                      tGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCt            >  2:874320/1‑150 (MQ=255)
                                                                                                                         cccAGCTGTTGCAGAGGGCGAACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTaa          <  1:310633/149‑1 (MQ=255)
                                                                                                                                tGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCgg  >  1:626888/1‑151 (MQ=255)
                                                                                                                                tGTTGCAGAGGGCATACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCg   <  1:425239/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                      
GGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGG  >  NC_000913/3350701‑3350978

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: