Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,122,627 G→A 7.8% G41G (GGC→GGT)  PP_4508 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,122,6270GA 7.8% 146.8 / 2.5 51G41G (GGC→GGT) PP_4508AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base A (2/2);  total (24/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 = 9.87e-01

CGAGGACGAAACTGAGGTCGGAGCTTTCGGCGATGATGCGCTGGTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCTTGAACATGGCTCTGCAATCCAGCCCTTCCAGTT  >  NC_002947/5122479‑5122775
                                                                                                                                                    |                                                                                                                                                    
cGAGGACGAAACTGAGGTCGGAGCTTTCGGCGATGATGCGCTGGTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAgg                                                                                                                                                      >  8:58396/1‑149 (MQ=255)
            tGAGGTCGGAGCTTTCGGCGATGATGCGCTGGTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCa                                                                                                                                          >  4:104642/1‑149 (MQ=255)
                gTCGGAGCTTTCGGCGATGATGCGCTGGTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGc                                                                                                                                      <  3:249006/149‑1 (MQ=255)
                 tCGGAGCTTTCGGCGATGATGCGCTGGTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGcc                                                                                                                                     <  1:26518/149‑1 (MQ=255)
                    gAGCTTTCGGCGATGATGCGCTGGTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGa                                                                                                                                  <  5:249922/149‑1 (MQ=255)
                             gCGATGATGCGCTGGTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGc                                                                                                                                                     >  7:212429/1‑121 (MQ=255)
                             gCGATGATGCGCTGGTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGc                                                                                                                                                     <  8:212429/121‑1 (MQ=255)
                             gCGATGATGCGCTGGTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCAc                                                                                                                         >  5:106227/1‑149 (MQ=255)
                                      cgcTGGTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGtt                                                                                                                 <  8:183132/148‑1 (MQ=255)
                                       gcTGGTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATCGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGc                                                                                                               >  6:10770/1‑149 (MQ=255)
                                           gTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCa                                                                                                             <  3:104642/147‑1 (MQ=255)
                                                gcACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCgg                                                                                                      <  1:304141/149‑1 (MQ=255)
                                                      ggCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATg                                                                                                <  6:267613/149‑1 (MQ=255)
                                                           tCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTc                                                                                           >  6:296753/1‑149 (MQ=255)
                                                               aGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCt                                                                                                                                    >  3:48830/1‑104 (MQ=255)
                                                               aGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCt                                                                                                                                    <  4:48830/104‑1 (MQ=255)
                                                                        ccAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGt                                                                              <  7:58396/149‑1 (MQ=255)
                                                                          aGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGt                                                                              <  6:183475/147‑1 (MQ=255)
                                                                          aGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGt                                                                              >  5:183475/1‑147 (MQ=255)
                                                                                 ggCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGACCGATGGCCTGATTGCCTGACAGCTCCACCGCCAGTTGCCATAgg                                                                                                         <  8:201320/113‑1 (MQ=255)
                                                                                 ggCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGACCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAgg                                                                                                         >  7:197384/1‑113 (MQ=255)
                                                                                 ggCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGACCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAgg                                                                                                         >  7:201320/1‑113 (MQ=255)
                                                                                 ggCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGACCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAgg                                                                                                         <  8:197384/113‑1 (MQ=255)
                                                                                  gCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGAt                                                                    >  8:54519/1‑149 (MQ=255)
                                                                                   cTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGAtc                                                                   >  7:211285/1‑149 (MQ=255)
                                                                                    tGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGAtcg                                                                  >  6:108819/1‑149 (MQ=255)
                                                                                         aTCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGAAGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAg                                                             >  1:97106/1‑149 (MQ=255)
                                                                                            aaCGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGggc                                                          >  5:44208/1‑149 (MQ=255)
                                                                                            aaCGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGggc                                                          <  7:54519/149‑1 (MQ=255)
                                                                                             aCGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCa                                                                                                                    >  1:55408/1‑90 (MQ=255)
                                                                                             aCGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCa                                                                                                                    <  2:55408/90‑1 (MQ=255)
                                                                                              cGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGggcgg                                                        <  2:281480/149‑1 (MQ=255)
                                                                                                  tAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGAGGAAAACGTGCGTCGTGATCGTCGAGGGCGGCa                                                      >  8:46647/1‑147 (MQ=255)
                                                                                                     cccACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTc                                                 >  4:141244/1‑149 (MQ=255)
                                                                                                        accacATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGa                                                              >  1:178477/1‑133 (MQ=255)
                                                                                                        accacATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGa                                                              <  2:178477/133‑1 (MQ=255)
                                                                                                        accacATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCaa                                              >  1:198528/1‑149 (MQ=255)
                                                                                                        accacATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCaa                                              <  3:273490/149‑1 (MQ=255)
                                                                                                              tGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCaa                                        <  6:44208/149‑1 (MQ=255)
                                                                                                               ggAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAAt                                       >  4:248962/1‑149 (MQ=255)
                                                                                                                 aaCGACGCCGGGCGGACCTCCCGTGCCTTGTTCAGGCGGTTGGCCTCATTGCCTGAGAGGTCCACCGCGTGTTGCCTTTGGGGGGTCATGGAGTCTTTCGGAAAACGTGCGTCGGGATCGTCGa                                                              <  3:52102/124‑1 (MQ=255)
                                                                                                                    gacgCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCtt                                  <  2:198528/149‑1 (MQ=255)
                                                                                                                        ccGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCt                                   >  5:270339/1‑144 (MQ=255)
                                                                                                                        ccGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCt                                   <  6:270339/144‑1 (MQ=255)
                                                                                                                         cGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCTTg                                 <  6:178962/145‑1 (MQ=255)
                                                                                                                         cGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCTTg                                 >  5:178962/1‑145 (MQ=255)
                                                                                                                            gCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCTTGAACATgg                          <  3:248962/149‑1 (MQ=255)
                                                                                                                                accaccCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCTTGAACATGGctct                      <  5:10770/149‑1 (MQ=255)
                                                                                                                                        tGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCTTGAACATGGCTCTGCAATc                >  8:246797/1‑147 (MQ=255)
                                                                                                                                           cATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCTTGAACATGGCTCTGCAATCCAGcc           <  7:246797/149‑1 (MQ=255)
                                                                                                                                                  aGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCTTGAACATGGCTCTGCAATCCAGCCCTTCCAg    <  3:136372/149‑1 (MQ=255)
                                                                                                                                                  aGGCCGATGGCCTCATTGCCAGACAGCTCCACCCCCAGAAGCCAGAGGCGGGTCATGCAGTCTTGCGGAAAACCAGCGTCGGGATCGTCGAGGGCGGCAAAGTCCACCCCCAAGTGCTTGTAAATGGCTCTCGCATCCAGCCCTTcccg    >  7:85498/1‑147 (MQ=255)
                                                                                                                                                   ggCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCTTGAACATGGCTCTGCTATCCAGCCCTTCCAGt   <  8:211285/149‑1 (MQ=255)
                                                                                                                                                    gCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCTTGAACATGGCTCTGCAATCCAGCCCTTCCAGtt  >  6:308681/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGAGGACGAAACTGAGGTCGGAGCTTTCGGCGATGATGCGCTGGTAGCGCACCAGGCGCTCGAAGCCCTCGGCCAGGGTGCGGCTGGACATCAACGCATAGCCCACCACATGGAACGACGCCGGGCGGACCACCCGTGCCATGTTCAGGCCGATGGCCTCATTGCCTGACAGCTCCACCGCCAGTTGCCATAGGCGGGTCATGGAGTCTTGCGGAAAACGTGCGTCGGGATCGTCGAGGGCGGCAAAGTCCAACCCCAATTGCTTGAACATGGCTCTGCAATCCAGCCCTTCCAGTT  >  NC_002947/5122479‑5122775

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