Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,430,159 C→T 8.7% S126L (TCG→TTG) ‡ 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,1590CT 8.7% 127.3 / 2.4 46S126L (TCG→TTG) ‡PP_3045ATP‑dependent protease ClpP
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base T (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.91e-01

AACACAAACGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAG  >  NC_002947/3430011‑3430307
                                                                                                                                                    |                                                                                                                                                     
aaCACAAACGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATc                                                                                                                                                       <  6:153216/149‑1 (MQ=255)
  cacaAACGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGa                                                                                                                                                     <  7:40556/149‑1 (MQ=255)
   acaAACGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGaa                                                                                                                                                    >  7:64514/1‑149 (MQ=255)
    caAACGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAAt                                                                                                                                                   <  2:254383/149‑1 (MQ=255)
     aaaCGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATg                                                                                                                                                  >  2:127628/1‑149 (MQ=255)
     aaaCGGAGGTGACCTGTTCGAGGGGCTTGCAATCAACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGACGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTATGGTCATGGCATCGAAt                                                                                                                                                   >  8:117086/1‑148 (MQ=255)
            ggTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATTCGATGCCATgac                                                                                                                                           <  6:182411/149‑3 (MQ=255)
              tGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCtga                                                                                                                                         <  8:208392/149‑1 (MQ=255)
                  cTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATc                                                                                                                                                       >  7:11773/1‑131 (MQ=255)
                  cTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATc                                                                                                                                                       <  8:11773/131‑1 (MQ=255)
                  cTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTgatgat                                                                                                                                     <  1:332457/149‑1 (MQ=255)
                       cGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAt                                                                                                                                 >  5:298203/1‑148 (MQ=255)
                                  ccATCCACATCGCGTTGAGCCGCTTGTGTGAATGTTGCACCGGACGAGTCGATGCCTTGGCCGCAAGTTCAGGCAGCGTGGCAGTCTGCGGAGCTCACCGTGTGGTCATTGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTacac                                                                                                                     <  7:319891/148‑1 (MQ=255)
                                    aTCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATTCGATGCCAt                                                                                                                                              <  8:217494/122‑1 (MQ=255)
                                    aTCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATTCGATGCCAt                                                                                                                                              >  7:217494/1‑122 (MQ=255)
                                    aTCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACACTTACACCTg                                                                                                                  >  5:283837/1‑149 (MQ=255)
                                       cacaACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGt                                                                                                               >  7:52824/1‑149 (MQ=255)
                                        acaACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGc                                                                                                                                            <  1:224864/120‑1 (MQ=255)
                                        acaACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGc                                                                                                                                            >  2:224864/1‑120 (MQ=255)
                                          aaCGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGATgg                                                                                                            >  2:28815/1‑149 (MQ=255)
                                          aaCGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGAGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGACATGGCATCGAAAGCCAAGCTGAAGAACCAA‑AACCCTAACACCTGGCTTgg                                                                                                            >  2:94584/1‑149 (MQ=255)
                                                cTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGAAGCGTGGCAGTATGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTTGGTTGGCGGTGa                                                                                                      <  7:117086/149‑1 (MQ=255)
                                                 ggAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGAt                                                                                                     <  6:226034/148‑1 (MQ=255)
                                                   aGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGAt                                                                                                     <  6:221152/147‑1 (MQ=255)
                                                                   gTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCgg                                                                                   >  3:12741/1‑149 (MQ=255)
                                                                    ttGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCggg                                                                                  >  3:73693/1‑149 (MQ=255)
                                                                      gCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGATGGCGGTGATGCCGAGGACCTTCGTCGGGt                                                                                 <  1:28815/148‑1 (MQ=255)
                                                                         ccGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCg                                                                             >  6:316801/1‑149 (MQ=255)
                                                                            gCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCa                                                                                                                                  <  5:315437/94‑1 (MQ=255)
                                                                            gCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCa                                                                                                                                  >  6:315437/1‑94 (MQ=255)
                                                                               gAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGt                                                                        <  1:94584/149‑1 (MQ=255)
                                                                                aGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATg                                                                                                    >  1:335597/1‑119 (MQ=255)
                                                                                aGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATg                                                                                                    <  2:335597/119‑1 (MQ=255)
                                                                                aGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGc                                                                      <  2:321000/149‑1 (MQ=255)
                                                                                  tCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCtt                                                                    <  1:47248/149‑1 (MQ=255)
                                                                                     aTGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCAGCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGCACCGTCGTCGGGTGGCGGACGTGCTTGAt                                                                 >  8:165260/1‑149 (MQ=255)
                                                                                              ccGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTc                                                        <  8:64514/149‑1 (MQ=255)
                                                                                              ccGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTc                                                        >  5:164791/1‑149 (MQ=255)
                                                                                                    gTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATc                                                                >  4:151988/1‑135 (MQ=255)
                                                                                                    gTGCAGGCAGCGTGCCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATc                                                                <  3:151988/135‑1 (MQ=255)
                                                                                                        aggcagCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGa                                              >  5:126837/1‑149 (MQ=255)
                                                                                                                tGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGcc                                      >  2:12778/1‑149 (MQ=255)
                                                                                                                     gTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATaa                                 <  6:164791/149‑1 (MQ=255)
                                                                                                                             cgcTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAgg                         <  4:73693/149‑1 (MQ=255)
                                                                                                                                  aCCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTc                      <  1:127628/147‑1 (MQ=255)
                                                                                                                                   ccGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCg                    >  7:287619/1‑148 (MQ=255)
                                                                                                                                            cATGGCATTCGATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGAt           >  3:56525/1‑148 (MQ=255)
                                                                                                                                                gCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGAGGGGCTTGCTCAGGCGTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGAcgc      <  4:12741/149‑1 (MQ=255)
                                                                                                                                                 cATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAAAGATCGCACCCTATAAGGCCAAGGCTCCGAACAACGATGAcgcg     >  2:12668/1‑149 (MQ=255)
                                                                                                                                                    cGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAg  >  4:183228/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACACAAACGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCAT‑AACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAG  >  NC_002947/3430011‑3430307

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