Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,430,203 G→C 17.0% S126S (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,2030GC17.0% 51.2 / 3.9 24S126S (TCG→TCC) ‡PP_3045ATP‑dependent protease ClpP
Reads supporting (aligned to +/- strand):  ref base G (7/12);  new base C (2/2);  total (9/15)
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.48e-01

CGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCT  >  NC_002947/3430062‑3430352
                                                                                                                                              |                                                                                                                                                      
cGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCAt                                                                                                                                                <  1:7037/150‑1 (MQ=255)
cGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCAt                                                                                                                                                <  1:548484/150‑1 (MQ=255)
  gAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATTCGATGCCATga                                                                                                                                              <  2:507922/150‑2 (MQ=255)
  gAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATTCGATGCCATga                                                                                                                                              <  2:235117/150‑2 (MQ=255)
     gTGCCCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCtga                                                                                                                                           <  1:490310/150‑1 (MQ=255)
          cTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATTCGATGCCATgac                                                                                                                                             >  1:105418/1‑141 (MQ=255)
          cTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATTCGATGCCATgac                                                                                                                                             <  2:105418/143‑3 (MQ=255)
          cTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATc                                                                                                                                      <  2:412237/150‑1 (MQ=255)
             ttCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAGTGCCATGCTGATGATc                                                                                                                                      >  2:405040/1‑147 (MQ=255)
             ttCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATc                                                                                                                                      <  1:405040/147‑1 (MQ=255)
                                  aaCGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGc                                                                                                              <  1:186420/150‑1 (MQ=255)
                                       gcTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGa                                                                                                         <  1:310244/150‑1 (MQ=255)
                                                        aaCGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACGCATCCACCTGGGTTGGCGGTGATGCCGAGGACCTtcgtc                                                                                        >  1:32310/1‑150 (MQ=255)
                                                        aaCGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTtcgtc                                                                                        >  1:313358/1‑150 (MQ=255)
                                                                                  ttGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGcc                                                              >  1:13288/1‑150 (MQ=255)
                                                                                                    cgcttGGCACTCTGCGGCGCTCACCGTGTGCTCATGGCCTCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGTCTTCGAGGCAATGATCGCa                                            <  2:263402/146‑1 (MQ=255)
                                                                                                                    gcgcTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCCAGGGCCTTCGTCGGGTGGCGGACGTGCTTGAGCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAgg                           <  1:531963/150‑1 (MQ=255)
                                                                                                                     gggtCACCGTGTGGTCAGGGCATCGAATGCCACGCTGATGATCCATAACCCTTACCCCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGCTCGCAGCCTATAAGGCCAAGGc                          <  2:32310/147‑1 (MQ=255)
                                                                                                                        tCACCCTTTGGTCATGGCATCTAATGCCATGCTGATGATCCATAACCCTTACACCTGAGTTGGCGGTGATGCCGAGTACCTTCTTCGGGTGGCGGAGCTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTcc                       <  1:400515/150‑1 (MQ=255)
                                                                                                                                tgGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCg               >  2:601714/1‑150 (MQ=255)
                                                                                                                                   gcATGGCATTCGATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCgg                                                                                      >  2:95400/2‑77 (MQ=255)
                                                                                                                                   gcATGGCATTCGATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCgg                                                                                      <  1:95400/76‑1 (MQ=255)
                                                                                                                                   gcATGGCATTCGATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATg            >  2:549663/2‑150 (MQ=255)
                                                                                                                                     aTGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAAAGGGCATGCGTGATAAGGCCAAGGCTCCGAACATCGATGAc          >  2:511697/1‑150 (MQ=255)
                                                                                                                                              gAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCt  <  2:313358/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
CGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCT  >  NC_002947/3430062‑3430352

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: