Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,430,349 G→C 7.7% S189S (TCG→TCC)  PP_3045 → ATP‑dependent protease ClpP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,430,3490GC 7.7% 159.6 / 4.3 52S189S (TCG→TCC) PP_3045ATP‑dependent protease ClpP
Reads supporting (aligned to +/- strand):  ref base G (25/23);  new base C (2/2);  total (27/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 = 5.34e-01

TGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGGTGAC  >  NC_002947/3430208‑3430495
                                                                                                                                             |                                                                                                                                                  
tGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCg                                                                                                                                             <  1:127836/149‑1 (MQ=255)
 gCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGt                                                                                                                                            >  7:10263/1‑149 (MQ=255)
  ccGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTg                                                                                                                                           <  5:127791/149‑1 (MQ=255)
  ccGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTg                                                                                                                                           <  5:205593/149‑1 (MQ=255)
   cGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTgg                                                                                                                                          >  3:223053/1‑149 (MQ=255)
   cGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTgg                                                                                                                                          >  3:116357/1‑149 (MQ=255)
    gAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCATCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCCGAGGCCGTgag                                                                                                                                         <  3:101263/149‑3 (MQ=255)
         ccTTCGTCGGGTGGCGGACGGGCTTGATCAGGCCTTCGAGGCAAGGCTCGCAGCCTATACGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGGCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTgg                                                                                                                                    <  3:42843/149‑1 (MQ=255)
          cTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTggg                                                                                                                                   <  8:10263/149‑1 (MQ=255)
          cTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATAGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTggg                                                                                                                                   >  4:178257/1‑149 (MQ=255)
           ttCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTgggg                                                                                                                                  >  5:19568/1‑149 (MQ=255)
                cgGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGc                                                                                                                             >  7:204608/1‑149 (MQ=255)
                       cGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAg                                                                                                                      <  2:154200/149‑1 (MQ=255)
                             tGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCGCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGt                                                                                                                <  4:158443/149‑1 (MQ=255)
                                  gATCAGGCCTTCGAGTCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATg                                                                                                            <  6:19568/148‑1 (MQ=255)
                                  gATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATg                                                                                                            <  5:187268/148‑1 (MQ=255)
                                    tCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATg                                                                                                            <  2:25842/146‑1 (MQ=255)
                                    tCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATg                                                                                                            >  1:25842/1‑146 (MQ=255)
                                     cAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGgtg                                                                                                         <  8:204608/148‑1 (MQ=255)
                                                      gATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGcc                                                                                       >  3:135725/1‑149 (MQ=255)
                                                          gCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTgg                                                                                   <  4:223053/149‑1 (MQ=255)
                                                          gCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTgg                                                                                   <  4:135725/149‑1 (MQ=255)
                                                             gCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGc                                                                                 >  5:275259/1‑148 (MQ=255)
                                                                 ataAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGgc                                                                            >  5:8306/1‑149 (MQ=255)
                                                                   aaGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGgcgc                                                                          >  6:79046/1‑149 (MQ=255)
                                                                     ggCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCgcgg                                                                        <  1:253353/149‑1 (MQ=255)
                                                                        cAAGGATCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTGAATGCCGAAAGTTGGCTCACGGACTCGGCGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGAGTCAAGGTTCAGGCGTg                                                                                         >  5:188973/1‑129 (MQ=255)
                                                                                         aTGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCaa                                                     >  2:65496/1‑148 (MQ=255)
                                                                                              gcgGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATAcccc                                               >  5:988/1‑149 (MQ=255)
                                                                                                  aGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGc                                             >  3:32193/1‑147 (MQ=255)
                                                                                                       agaCGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTAc                                                          <  6:70583/129‑1 (MQ=255)
                                                                                                       agaCGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTAc                                                          >  5:70583/1‑129 (MQ=255)
                                                                                                       agaCGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGcc                                        <  1:65496/147‑1 (MQ=255)
                                                                                                         aCGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAt                                    >  2:116641/1‑149 (MQ=255)
                                                                                                          cGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGTCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAtt                                   <  8:49116/149‑1 (MQ=255)
                                                                                                          cGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAtt                                   <  1:116641/149‑1 (MQ=255)
                                                                                                            tCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTg                                  >  5:274692/1‑148 (MQ=255)
                                                                                                              tGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGa                                >  4:144810/1‑148 (MQ=255)
                                                                                                              tGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGa                                >  6:9950/1‑148 (MQ=255)
                                                                                                              tGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGAt                               >  2:250859/1‑149 (MQ=255)
                                                                                                                gTCAATGCCGAAAGGTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGAt                               <  5:9950/147‑1 (MQ=255)
                                                                                                                  cAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCaa                                           >  5:198865/1‑133 (MQ=255)
                                                                                                                  cAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCaa                                           <  6:198865/133‑1 (MQ=255)
                                                                                                                       ccGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGTGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAg                      <  4:116357/149‑1 (MQ=255)
                                                                                                                               agcgccACGGCCTCCGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTc                         <  7:222248/133‑1 (MQ=255)
                                                                                                                               agcgccACGGCCTCCGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTc                         >  8:222248/6‑138 (MQ=255)
                                                                                                                                 cgccACGGCCTCCGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCGATTGGATCAGCTCAAGGCCAGTAggg            >  3:175873/4‑149 (MQ=255)
                                                                                                                                      cGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATggg       <  6:274692/149‑1 (MQ=255)
                                                                                                                                       ggCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGc                                         >  2:112174/1‑114 (MQ=255)
                                                                                                                                           tCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTa               >  2:55116/1‑136 (MQ=255)
                                                                                                                                           tCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTa               <  1:55116/136‑1 (MQ=255)
                                                                                                                                           tCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGGTGAc  >  6:162823/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
TGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGGTGAC  >  NC_002947/3430208‑3430495

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