Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,389,305 C→T 7.3% L41L (CTC→CTT)  rlmKL → 23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,389,3050CT 7.3% 161.3 / 2.9 55L41L (CTC→CTT) rlmKL23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (29/22);  new base T (2/2);  total (31/24)
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.33e-01

GGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTA  >  NC_002947/2389165‑2389454
                                                                                                                                            |                                                                                                                                                     
ggCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGAc                                                                                                                                              >  2:510531/1‑150 (MQ=255)
ggCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGAc                                                                                                                                              >  1:513731/1‑150 (MQ=255)
 gCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACg                                                                                                                                             >  1:127353/1‑150 (MQ=255)
 gCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACg                                                                                                                                             >  1:783580/1‑150 (MQ=255)
            ggggCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAg                                                                                                                                  <  2:127353/150‑1 (MQ=255)
                   tGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGcacaaca                                                                                                                           <  2:783580/150‑4 (MQ=255)
                   tGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCt                                                                                                                           >  1:244154/1‑150 (MQ=255)
                      gTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGCGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCgg                                                                                                                        <  1:54240/150‑1 (MQ=255)
                      gTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCgg                                                                                                                        >  2:252510/1‑150 (MQ=255)
                         gCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGa                                                                                                                                   <  2:20938/136‑1 (MQ=255)
                         gCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGa                                                                                                                                   >  1:20938/1‑136 (MQ=255)
                           cccGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAAGACCTCGGCCTTGACGAGGTTCGCGAGCACACCTCGGCCAtt                                                                                                                   >  2:779339/1‑150 (MQ=255)
                                  gtcaAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCgcggcg                                                                                                            >  2:588201/1‑150 (MQ=255)
                                      aaTAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGcgccgc                                                                                                        <  2:513731/150‑1 (MQ=255)
                                        tAGCCCAAGGTGCGAGAATACTCATGTCGGACCGTCTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGccgccg                                                                                                      <  2:191541/150‑1 (MQ=255)
                                         aGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGa                                                                                                     <  1:779339/150‑1 (MQ=255)
                                            ccAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACAt                                                                                                  >  2:288749/1‑150 (MQ=255)
                                             cAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATg                                                                                                 >  1:288807/1‑150 (MQ=255)
                                              aaGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATgg                                                                                                <  1:483550/150‑1 (MQ=255)
                                                   gcTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCGCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACt                                                                                           <  2:244154/150‑1 (MQ=255)
                                                         aTACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTAc                                                                                     <  1:252510/150‑1 (MQ=255)
                                                             tCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCgcgg                                                                                                              >  2:559620/1‑121 (MQ=255)
                                                             tCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCgcgg                                                                                                              <  1:559620/121‑1 (MQ=255)
                                                                    ggACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTgtgtgt                                                                          >  2:835894/1‑150 (MQ=255)
                                                                    ggACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTgtgtgt                                                                          >  1:460786/1‑150 (MQ=255)
                                                                     gACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTc                                                                         >  2:410648/1‑150 (MQ=255)
                                                                      aCCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCtt                                                                                       <  2:810560/135‑1 (MQ=255)
                                                                      aCCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCtt                                                                                       >  1:810560/1‑135 (MQ=255)
                                                                          tttCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTc                                                                    >  2:428339/1‑150 (MQ=255)
                                                                            tCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCAc                                                                  <  2:288807/150‑1 (MQ=255)
                                                                              gAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGt                                                                >  1:32774/1‑150 (MQ=255)
                                                                                       cTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAAc                                                       >  2:546632/1‑150 (MQ=255)
                                                                                          aCCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGAc                                                                                                                                              >  2:618625/1‑60 (MQ=255)
                                                                                          aCCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGAc                                                                                                                                              <  1:618625/60‑1 (MQ=255)
                                                                                          aCCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACcgc                                                    <  2:364766/150‑1 (MQ=255)
                                                                                            cTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTg                                                                       >  1:430208/1‑129 (MQ=255)
                                                                                            cTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTg                                                                       <  2:430208/129‑1 (MQ=255)
                                                                                                cccAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTg                                              >  1:25465/1‑150 (MQ=255)
                                                                                                     aGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGt                                         <  1:546632/150‑1 (MQ=255)
                                                                                                         cTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGc                                                                                         <  1:227817/98‑1 (MQ=255)
                                                                                                         cTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGc                                                                                         >  2:227817/1‑98 (MQ=255)
                                                                                                            gAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAg                                  <  1:288749/150‑1 (MQ=255)
                                                                                                                     cTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCg                         >  2:810879/1‑150 (MQ=255)
                                                                                                                            tcaagGCCGAGGCCCTTGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTg                                                                       <  1:123247/94‑1 (MQ=255)
                                                                                                                            tcaagGCCGAGGCCCTTGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTg                                                                       >  2:123247/4‑97 (MQ=255)
                                                                                                                              aagGCCGAGGCCCTTGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGAt                       <  1:528620/142‑1 (MQ=255)
                                                                                                                              aagGCCGAGGCCCTTGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGAt                       >  2:528620/2‑143 (MQ=255)
                                                                                                                                   cAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCcga           <  1:510531/150‑1 (MQ=255)
                                                                                                                                    aaGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCcgac          >  2:59877/1‑150 (MQ=255)
                                                                                                                                    aaGGGCCTCGGCCTTGACGAGGTGCGCGAGAACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGCTCACTTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCcgac          >  1:355355/1‑150 (MQ=255)
                                                                                                                                     aGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTgg                                          >  1:695302/1‑117 (MQ=255)
                                                                                                                                     aGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTgg                                          <  2:695302/117‑1 (MQ=255)
                                                                                                                                      gggCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCcgacga        >  2:31717/1‑150 (MQ=255)
                                                                                                                                          cTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATgaag                   <  2:481430/135‑1 (MQ=255)
                                                                                                                                          cTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATgaag                   >  1:481430/1‑135 (MQ=255)
                                                                                                                                            cGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTa  <  1:410648/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
GGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTA  >  NC_002947/2389165‑2389454

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