Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,852,594 T→C 7.8% intergenic (‑265/+32) 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,5940TC 7.8% 141.3 / 2.8 51intergenic (‑265/+32)recR/ybaBrecombination protein RecR/nucleoid‑associated broad specificity regulator
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base C (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 = 8.49e-01

GCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTCAGCGCAGCGGCAAT  >  NC_002947/4852453‑4852739
                                                                                                                                              |                                                                                                                                                 
gCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAgcg                                                                                                                                            >  5:175265/1‑149 (MQ=255)
        gcccccccGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGCCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCCCTGTTGTCAGCGGCGGtttt                                                                                                                                    <  6:347601/145‑1 (MQ=255)
              ccGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTAt                                                                                                                              >  5:244155/1‑149 (MQ=255)
              ccGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTAt                                                                                                                              >  3:68660/1‑149 (MQ=255)
                      atgagctgcGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGCACGCTCCAACAAAAAAGCACACAGAAACTGGGGTGAAATCAACGC‑AAAAAAAACCGCCACAGTTGTCAGAGGCGGTTATGAGTATCGCAAtcc                                                                                                                      >  7:232508/5‑147 (MQ=25)
                              cGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAg                                                                                                                <  2:238555/147‑1 (MQ=255)
                                 aGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACgg                                                                                                           >  7:20831/1‑149 (MQ=255)
                                    tGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCAt                                                                                                        >  8:5535/1‑149 (MQ=255)
                                     agagttccGGGGCTGCTTCGCAGCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATc                                                                                                       <  5:83179/143‑1 (MQ=37)
                                      aaaTCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCt                                                                                                      >  1:184893/1‑149 (MQ=255)
                                        aTCCCGGGGCCGCTTCGCGCCCCTTTCCCGACACAACGCCCCTCCCACAAAACAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTg                                                                                                    <  4:271741/149‑1 (MQ=255)
                                            cGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGc                                                                                                >  1:67608/1‑149 (MQ=255)
                                                    cTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCggcggc                                                                                        <  6:175265/149‑1 (MQ=255)
                                                        gcgcCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCt                                                                                    <  2:184893/149‑1 (MQ=255)
                                                          gcCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCt                                                                                    >  3:234245/1‑147 (MQ=255)
                                                          gcCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCt                                                                                    <  4:234245/147‑1 (MQ=255)
                                                           ccccTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGc                                                                                  >  6:288985/1‑148 (MQ=255)
                                                           ccccTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGc                                                                                  <  5:288985/148‑1 (MQ=255)
                                                           ccccTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCa                                                                                 >  1:82047/1‑149 (MQ=255)
                                                             ccTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCAt                                                                                >  4:291562/1‑149 (MQ=255)
                                                             ccTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAATGGCATCTTGAAGCCCGGCGGCAGCTGCATg                                                                               >  8:245637/1‑149 (MQ=255)
                                                               tttCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGcc                                                                             >  2:261281/1‑149 (MQ=255)
                                                                    cgACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAATGGCATCTTGAAGCCCGGCGGCAGCTGCATGCcggcgg                                                                        <  7:245637/149‑1 (MQ=255)
                                                                      acacAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTc                                                                      <  4:68660/149‑1 (MQ=255)
                                                                        acaAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCAt                                                                    <  2:67608/149‑1 (MQ=255)
                                                                          aaGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATgc                                                                  >  1:130/1‑149 (MQ=255)
                                                                            ggCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATgccg                                                                <  5:218318/149‑1 (MQ=255)
                                                                             gCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTGTTGAGTATCGCAATGCCATTAGAACGGCAGCTTGAAGCCCGGCGGCAGATGCATGCCGGCGTGCATgacgc                                                               >  7:328099/1‑149 (MQ=255)
                                                                                   cccACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATc                                                          <  3:291562/149‑1 (MQ=255)
                                                                                    ccACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCcggcg                                                                         <  4:85800/132‑1 (MQ=255)
                                                                                    ccACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCcggcg                                                                         >  3:85800/1‑132 (MQ=255)
                                                                                       caAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTc                                                     <  6:244155/149‑1 (MQ=255)
                                                                                        aaaaaaGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGaa                                                                                                  <  5:303884/103‑1 (MQ=255)
                                                                                        aaaaaaGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGaa                                                                                                  >  6:303884/1‑103 (MQ=255)
                                                                                        aaaaaaGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTcc                                                    >  7:339273/1‑149 (MQ=255)
                                                                                         caacagcagcCCGCAAGTGGGTTGAAATCTCCGC‑AAAAAAATCCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCt                                                   <  8:232508/145‑1 (MQ=255)
                                                                                           aaaGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGa                                                 <  1:261281/149‑1 (MQ=255)
                                                                                                   ccGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGc                                                                                  >  6:263204/1‑108 (MQ=255)
                                                                                                   ccGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGc                                                                                  <  5:263204/108‑1 (MQ=255)
                                                                                                      cAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACtgctctgctc                                      <  7:5535/149‑1 (MQ=255)
                                                                                                         cTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGAc                                   >  8:28425/1‑149 (MQ=255)
                                                                                                          tGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCAGGCGGCAGCTGCATGCCGGCGGTCATGCAGCCCATCATTTACTGACTGCTCTGAGCGAcc                                  >  8:73586/1‑149 (MQ=255)
                                                                                                            ggTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGAcc                                  <  2:316825/147‑1 (MQ=255)
                                                                                                                         gC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGcgtcg                   >  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              >  4:267245/1‑148 (MQ=255)
                                                                                                                                aaaaCCGCCGCTGACAACAGTGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTcag              <  3:267245/148‑1 (MQ=255)
                                                                                                                                aaaaCCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTcagc             <  8:339273/149‑1 (MQ=255)
                                                                                                                                          cTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCAACTTGAAGCCCGGCGGCAGCTGCATGCCTGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCGCTGCGTCGTTCAGCGCAGCGGCaa   >  1:80293/1‑149 (MQ=255)
                                                                                                                                           tgttgtCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTAGACCTTGCGCACTGCGTCGTTCAGCGCAGCGGCAAt  >  6:79703/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTCAGCGCAGCGGCAAT  >  NC_002947/4852453‑4852739

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: ‑