Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,389,301 G→C 7.9% G40A (GGC→GCC)  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,3010GC 7.9% 148.8 / 5.1 52G40A (GGC→GCC) rlmKL23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base C (2/2);  total (29/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 = 2.63e-01

GGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGA  >  NC_002947/2389165‑2389448
                                                                                                                                        |                                                                                                                                                   
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                                                                                                                             >  1:20938/1‑136 (MQ=255)
                         gCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGa                                                                                                                             <  2:20938/136‑1 (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                                                                                                        <  1:559620/121‑1 (MQ=255)
                                                             tCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCgcgg                                                                                                        >  2:559620/1‑121 (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                                                                                   >  2:227817/1‑98 (MQ=255)
                                                                                                         cTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGc                                                                                   <  1:227817/98‑1 (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                 >  2:528620/2‑143 (MQ=255)
                                                                                                                              aagGCCGAGGCCCTTGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGAt                 <  1:528620/142‑1 (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)
                                                                                                                                        |                                                                                                                                                   
GGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGA  >  NC_002947/2389165‑2389448

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