Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,847,988 C→T 8.1% L54L (CTG→TTG)  rlmI → ribosomal RNA large subunit methyltransferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,847,9880CT 8.1% 137.7 / 2.6 49L54L (CTG→TTG) rlmIribosomal RNA large subunit methyltransferase I
Reads supporting (aligned to +/- strand):  ref base C (21/24);  new base T (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

GCCTTCGCCTCAAAGCCAATGCCGACCGCCGCCTGCGCGCCGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAGCGCCTGTTCGAGCAGCCG  >  NC_002947/5847842‑5848131
                                                                                                                                                  |                                                                                                                                               
gCCTTCGCCTCAAAGCCAATGCCGACCGCCGCCTGCGCGCCGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTg                                                                                                                                               <  5:185019/149‑1 (MQ=255)
gCCTTCGCCTCAAAGCCAATGCCGACCGCCGCCTGCGCGCCGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTg                                                                                                                                               <  4:97233/149‑1 (MQ=255)
      gCCTCAAAGCCAATGCCGACCGCGGCCTGCGCGCCGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGGCAACGGCAAGCCGCTGGGCAtc                                                                                                                                         <  6:59893/149‑1 (MQ=255)
           aaaGCCAATGCCGACCGCCGCCTGCGCGCCGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTcgc                                                                                                                                    <  8:53679/149‑1 (MQ=255)
              gCCAATGCCGACCGCCGCCTGCGCGCCGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCt                                                                                                                                 >  3:54564/1‑149 (MQ=255)
               ccAATGCCGACCGCCGCCTGCGCGCCGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCAtcgt                                                                                                                                       <  8:720/142‑1 (MQ=255)
               ccAATGCCGACCGCCGCCTGCGCGCCGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCAtcgt                                                                                                                                       >  7:720/1‑142 (MQ=255)
                 aaTGCCGACCGCCGCCTGCGCGCCGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAg                                                                                                                              <  1:57434/149‑1 (MQ=255)
                       gACCGCCGCCTGCGCGCCGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGaa                                                                                                                        <  2:267662/149‑1 (MQ=255)
                                    gcgcCGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTgc                                                                                                           <  7:120040/149‑1 (MQ=255)
                                       ccGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGa                                                                                                                         <  4:203858/132‑1 (MQ=255)
                                       ccGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGa                                                                                                                         >  3:203858/1‑132 (MQ=255)
                                            cACCTGTGGGTCTACAGCAACTACGTCGATGTCAGCGCGACCCCTTTGCAAGGCTGCCAGGCGGGCCAGCAGGCATTTCGCGAGGCAGCCAACGGCAAGCCGCTGGGATTCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCct                                                                                                   <  1:36084/149‑1 (MQ=255)
                                              ccTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCctgc                                                                                                 >  7:184547/1‑149 (MQ=255)
                                                    ggTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCa                                                                                            <  1:197062/148‑1 (MQ=255)
                                                      tCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAAGCTGATCTGCGCCCGCATGCTGTCAcgc                                                                                         >  2:211126/1‑149 (MQ=255)
                                                      tCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGcgcccg                                                                                                      >  3:11175/1‑136 (MQ=255)
                                                      tCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGcgcccg                                                                                                      <  4:11175/136‑1 (MQ=255)
                                                      tCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCAcgc                                                                                         <  4:54564/149‑1 (MQ=255)
                                                              aaCGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCAcgc                                                                                         <  6:103662/141‑1 (MQ=255)
                                                              aaCGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCAcgc                                                                                         >  5:103662/1‑141 (MQ=255)
                                                                  aaGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTg                                                                             >  5:253040/1‑149 (MQ=255)
                                                                   aGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTg                                                                                                                                <  4:147107/97‑1 (MQ=255)
                                                                   aGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTg                                                                                                                                >  3:147107/1‑97 (MQ=255)
                                                                       gATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGTGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGc                                                                         <  1:77726/148‑1 (MQ=255)
                                                                          gTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGa                                                                     >  2:259673/1‑149 (MQ=255)
                                                                                gcgACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTc                                                               >  5:150924/1‑149 (MQ=255)
                                                                                gcgACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTc                                                               >  6:176799/1‑149 (MQ=255)
                                                                                     cccTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCgctgc                                                          >  7:231975/1‑149 (MQ=255)
                                                                                          tGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGt                                                      <  6:150924/148‑1 (MQ=255)
                                                                                           gCAAGGCTTCCAGGCGGGCCAGAAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTACGAACAACCTGATCTGCGCACGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTa                                                     >  6:172856/1‑148 (MQ=255)
                                                                                                  ttCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGa                                                                                       <  3:101291/107‑1 (MQ=255)
                                                                                                  ttCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGa                                                                                       >  4:101291/1‑107 (MQ=255)
                                                                                                  ttCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCt                                             >  6:23682/1‑149 (MQ=255)
                                                                                                   tCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCt                                             >  3:226368/1‑148 (MQ=255)
                                                                                                                      cATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAAc                                         <  2:92298/133‑1 (MQ=255)
                                                                                                                      cATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTgc                         >  2:188323/1‑149 (MQ=255)
                                                                                                                             gAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAgcg                  >  3:20/1‑149 (MQ=255)
                                                                                                                             gAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAgc                   <  5:172856/148‑1 (MQ=255)
                                                                                                                             gAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGAACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAgcg                  >  2:158435/1‑149 (MQ=255)
                                                                                                                                   gCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAGCGCCTGtt            <  8:184547/149‑1 (MQ=255)
                                                                                                                                    ccAGCGGCTTGCCGTTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAgc                   <  3:117743/141‑1 (MQ=255)
                                                                                                                                    ccAGCGGCTTGCCGTTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAgc                   >  4:117743/1‑141 (MQ=255)
                                                                                                                                    ccAGCGGCTTGCCGTTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAGCGCCTGTTc           >  8:185605/1‑149 (MQ=255)
                                                                                                                                    ccAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAGCGCCTGTTc           <  4:226368/149‑1 (MQ=255)
                                                                                                                                        cggcttgCCGTTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAGCGCCTGTTCGagc       <  7:185605/143‑1 (MQ=255)
                                                                                                                                             aGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAGCGCCTGTTCGAGCAGCCg  <  7:252969/149‑1 (MQ=255)
                                                                                                                                                 gCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAGCGCCTGTTCGAGCAGCCg  <  6:236401/145‑1 (MQ=255)
                                                                                                                                                 gCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAGCGCCTGTTCGAGCAGCCg  >  5:236401/1‑145 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCCTTCGCCTCAAAGCCAATGCCGACCGCCGCCTGCGCGCCGGCCACCTGTGGGTCTACAGCAACGAAGTCGATGTCAGCGCGACCCCTTTGCAAGGCTTCCAGGCGGGCCAGCAGGCCATTCTCGAGGCAGCCAACGGCAAGCCGCTGGGCATCGTCGCGCTGAGTCCGAACAACCTGATCTGCGCCCGCCTGCTGTCACGCGACATCAAGCTGCCGCTGGACAAGTCGCTGCTGGTACACCGCCTGAACGTCGCCCTGTCGCTGCGCGAGCGCCTGTTCGAGCAGCCG  >  NC_002947/5847842‑5848131

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