Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,852,596 T→A 7.9% intergenic (‑267/+30) recR ← / ← ybaB recombination protein RecR/nucleoid‑associated broad specificity regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,852,5960TA 7.9% 139.2 / 3.7 51intergenic (‑267/+30)recR/ybaBrecombination protein RecR/nucleoid‑associated broad specificity regulator
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base A (2/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

GCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTCAGCGCAGCGGCAATG  >  NC_002947/4852453‑4852740
                                                                                                                                                |                                                                                                                                                
gCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcg                                                                                                                                             >  5:175265/1‑149 (MQ=255)
        gcccccccGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGCCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCCCTGTTGTCAGCGGCGGtttt                                                                                                                                     <  6:347601/145‑1 (MQ=255)
              ccGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTAt                                                                                                                               >  5:244155/1‑149 (MQ=255)
              ccGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTAt                                                                                                                               >  3:68660/1‑149 (MQ=255)
                      atgagctgcGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGCACGCTCCAACAAAAAAGCACACAGAAACTGGGGTGAAATCAACGCAAAAAAAACCGCCACAGTTGTCAGAGGCGGTTATGAGTATCGCAAtcc                                                                                                                       >  7:232508/5‑147 (MQ=25)
                              cGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAg                                                                                                                 <  2:238555/147‑1 (MQ=255)
                                 aGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACgg                                                                                                            >  7:20831/1‑149 (MQ=255)
                                    tGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCAt                                                                                                         >  8:5535/1‑149 (MQ=255)
                                     agagttccGGGGCTGCTTCGCAGCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATc                                                                                                        <  5:83179/143‑1 (MQ=37)
                                      aaaTCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCt                                                                                                       >  1:184893/1‑149 (MQ=255)
                                        aTCCCGGGGCCGCTTCGCGCCCCTTTCCCGACACAACGCCCCTCCCACAAAACAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTg                                                                                                     <  4:271741/149‑1 (MQ=255)
                                            cGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGc                                                                                                 >  1:67608/1‑149 (MQ=255)
                                                    cTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCggcggc                                                                                         <  6:175265/149‑1 (MQ=255)
                                                        gcgcCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCt                                                                                     <  2:184893/149‑1 (MQ=255)
                                                          gcCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCt                                                                                     >  3:234245/1‑147 (MQ=255)
                                                          gcCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCt                                                                                     <  4:234245/147‑1 (MQ=255)
                                                           ccccTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGc                                                                                   >  6:288985/1‑148 (MQ=255)
                                                           ccccTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGc                                                                                   <  5:288985/148‑1 (MQ=255)
                                                           ccccTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCa                                                                                  >  1:82047/1‑149 (MQ=255)
                                                             ccTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCAt                                                                                 >  4:291562/1‑149 (MQ=255)
                                                             ccTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAATGGCATCTTGAAGCCCGGCGGCAGCTGCATg                                                                                >  8:245637/1‑149 (MQ=255)
                                                               tttCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGcc                                                                              >  2:261281/1‑149 (MQ=255)
                                                                    cgACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAATGGCATCTTGAAGCCCGGCGGCAGCTGCATGCcggcgg                                                                         <  7:245637/149‑1 (MQ=255)
                                                                      acacAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTc                                                                       <  4:68660/149‑1 (MQ=255)
                                                                        acaAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCAt                                                                     <  2:67608/149‑1 (MQ=255)
                                                                          aaGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATgc                                                                   >  1:130/1‑149 (MQ=255)
                                                                            ggCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATgccg                                                                 <  5:218318/149‑1 (MQ=255)
                                                                             gCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTGTTGAGTATCGCAATGCCATTAGAACGGCAGCTTGAAGCCCGGCGGCAGATGCATGCCGGCGTGCATgacgc                                                                >  7:328099/1‑149 (MQ=255)
                                                                                   cccACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATc                                                           <  3:291562/149‑1 (MQ=255)
                                                                                    ccACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCcggcg                                                                          <  4:85800/132‑1 (MQ=255)
                                                                                    ccACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCcggcg                                                                          >  3:85800/1‑132 (MQ=255)
                                                                                       caAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTc                                                      <  6:244155/149‑1 (MQ=255)
                                                                                        aaaaaaGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGaa                                                                                                   <  5:303884/103‑1 (MQ=255)
                                                                                        aaaaaaGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGaa                                                                                                   >  6:303884/1‑103 (MQ=255)
                                                                                        aaaaaaGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTcc                                                     >  7:339273/1‑149 (MQ=255)
                                                                                         caacagcagcCCGCAAGTGGGTTGAAATCTCCGCAAAAAAATCCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCt                                                    <  8:232508/145‑1 (MQ=255)
                                                                                           aaaGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGa                                                  <  1:261281/149‑1 (MQ=255)
                                                                                                   ccGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGc                                                                                   >  6:263204/1‑108 (MQ=255)
                                                                                                   ccGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGc                                                                                   <  5:263204/108‑1 (MQ=255)
                                                                                                      cAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACtgctctgctc                                       <  7:5535/149‑1 (MQ=255)
                                                                                                         cTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGAc                                    >  8:28425/1‑149 (MQ=255)
                                                                                                          tGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCAGGCGGCAGCTGCATGCCGGCGGTCATGCAGCCCATCATTTACTGACTGCTCTGAGCGAcc                                   >  8:73586/1‑149 (MQ=255)
                                                                                                            ggTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGAcc                                   <  2:316825/147‑1 (MQ=255)
                                                                                                                         gCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGcgtcg                    >  1:305096/1‑149 (MQ=255)
                                                                                                                             aaaaaaaCCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTc                 <  2:130/149‑1 (MQ=255)
                                                                                                                                aaaaCCGCCGCTGACAACAGTGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTgc                               <  6:371167/132‑1 (MQ=255)
                                                                                                                                aaaaCCGCCGCTGACAACAGTGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTgc                               >  5:371167/1‑132 (MQ=255)
                                                                                                                                aaaaCCGCCGCTGACAACAGTGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTcag               <  3:267245/148‑1 (MQ=255)
                                                                                                                                aaaaCCGCCGCTGACAACAGTGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTcag               >  4:267245/1‑148 (MQ=255)
                                                                                                                                aaaaCCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTcagc              <  8:339273/149‑1 (MQ=255)
                                                                                                                                          cTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCAACTTGAAGCCCGGCGGCAGCTGCATGCCTGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCGCTGCGTCGTTCAGCGCAGCGGCaa    >  1:80293/1‑149 (MQ=255)
                                                                                                                                           tgttgtCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTAGACCTTGCGCACTGCGTCGTTCAGCGCAGCGGCAAt   >  6:79703/1‑149 (MQ=255)
                                                                                                                                               gtCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGTCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTCAGCGCAGCGGCAATg  <  3:5569/146‑1 (MQ=255)
                                                                                                                                               gtCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGTCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTCAGCGCAGCGGCAATg  >  4:5569/1‑146 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTCAGCGCAGCGGCAATG  >  NC_002947/4852453‑4852740

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