Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,430,202 C→T 16.2% 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,2020CT16.2% 48.0 / 2.3 24S126L (TCG→TTG) ‡PP_3045ATP‑dependent protease ClpP
Reads supporting (aligned to +/- strand):  ref base C (7/13);  new base T (2/2);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGAC  >  NC_002947/3430062‑3430344
                                                                                                                                             |                                                                                                                                               
cGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCAt                                                                                                                                        <  1:548484/150‑1 (MQ=255)
cGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCAt                                                                                                                                        <  1:7037/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)
                                                                                                                                             |                                                                                                                                               
CGGAGGTGACCTGTTCGAGGGGCTTGCCATCCACAACGCGCTGAGCCGCTTGGGTGAACGTTGCACCGGCCGAGTCGATGCCTTGGCCGCAAGTGCAGGCAGCGTGGCAGTCTGCGGCGCTCACCGTGTGGTCATGGCATCGAATGCCATGCTGATGATCCATAACCCTTACACCTGGGTTGGCGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGAC  >  NC_002947/3430062‑3430344

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: