Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,852,600 C→T 7.5% intergenic (‑271/+26) 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,6000CT 7.5% 152.8 / 2.2 54intergenic (‑271/+26)recR/ybaBrecombination protein RecR/nucleoid‑associated broad specificity regulator
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base T (2/2);  total (27/27)
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.60e-01

GCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTCAGCGCAGCGGCAATGAG  >  NC_002947/4852453‑4852742
                                                                                                                                                    |                                                                                                                                              
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                                                                            >  3:85800/1‑132 (MQ=255)
                                                                                    ccACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCcggcg                                                                            <  4:85800/132‑1 (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                                                                                     <  5:263204/108‑1 (MQ=255)
                                                                                                   ccGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGc                                                                                     >  6:263204/1‑108 (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                 <  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)
                                                                                                                                               ggcAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTTCTCTGCTCGACCTTGCGCACTGCGTCGTTCAGCGCAGCGGCAATGAg  <  7:73586/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGC‑AAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCATGCCGGCGGTCATGCCGCCCATCTTTTCCTGACTGCTCTGCTCGACCTTGCGCACTGCGTCGTTCAGCGCAGCGGCAATGAG  >  NC_002947/4852453‑4852742

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