Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,546,502 T→C 7.5% V10V (GTA→GTG)  rlmB ← 23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,546,5020TC 7.5% 156.1 / 2.2 53V10V (GTA→GTG) rlmB23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Reads supporting (aligned to +/- strand):  ref base T (21/28);  new base C (2/2);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TCTCACGGCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCACTGCTGC  >  NC_002947/5546354‑5546643
                                                                                                                                                    |                                                                                                                                              
tctcACGGCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTa                                                                                                                                               <  1:544642/150‑1 (MQ=255)
tctcACGGCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGACAGGGCTGCAGCAATGCCTGCACGGCGTGTa                                                                                                                                               <  1:937942/150‑1 (MQ=255)
       gcgcTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTa                                                                                                                                        <  2:730479/150‑1 (MQ=255)
             ggCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGAtttt                                                                                                                                  >  2:477743/1‑150 (MQ=255)
              gCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGAttttt                                                                                                                                 <  1:131012/150‑1 (MQ=255)
                   gCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCa                                                                                                                             >  2:60993/1‑149 (MQ=255)
                   gCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCa                                                                                                                             <  1:60993/149‑1 (MQ=255)
                    ccGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGc                                                                                                                           <  2:812173/150‑1 (MQ=255)
                      gACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTcc                                                                                                                              >  1:627650/1‑145 (MQ=255)
                      gACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTcc                                                                                                                              <  2:627650/145‑1 (MQ=255)
                        ccGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGAc                                                                                                                       <  1:355293/150‑1 (MQ=255)
                         cGGCACGCGGTTTTCCGCGGCCCGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgca                                                                                                                                      <  1:130335/134‑4 (MQ=255)
                         cGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgca                                                                                                                                      >  2:130335/1‑131 (MQ=255)
                               gcgGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGa                                                                                                                        <  2:42164/142‑1 (MQ=255)
                               gcgGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGa                                                                                                                        >  1:42164/1‑142 (MQ=255)
                                    tttCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGTCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATTGTGCTACAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGct                                                                                                          >  2:464363/1‑150 (MQ=255)
                                           ggCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGttt                                                                                                   <  1:236715/150‑1 (MQ=255)
                                               agccccagcAAGGTTTCCAGCCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGc                                                                                               <  1:412364/150‑1 (MQ=255)
                                                gcgccagcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCg                                                                                              >  2:797179/1‑150 (MQ=255)
                                                 cgccagcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGt                                                                                             >  2:884590/1‑150 (MQ=255)
                                                      gcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAgc                                                                                        >  1:163176/1‑150 (MQ=255)
                                                       cAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAgcg                                                                                       >  1:205390/1‑150 (MQ=255)
                                                         aGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTGTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGcc                                                                                     <  2:53270/150‑1 (MQ=255)
                                                          ggTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGc                                                                                               >  1:1112353/1‑139 (MQ=255)
                                                          ggTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGc                                                                                               <  2:1112353/139‑1 (MQ=255)
                                                           gTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGctc                                                                                                         >  1:375283/1‑128 (MQ=255)
                                                           gTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGctc                                                                                                         <  2:375283/128‑1 (MQ=255)
                                                                 aGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTcg                                                                             <  1:797179/150‑1 (MQ=255)
                                                                   gCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTcggc                                                                           >  1:73210/1‑150 (MQ=255)
                                                                         cTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCt                                                                     <  1:884590/150‑1 (MQ=255)
                                                                         cTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCt                                                                     >  2:934668/1‑150 (MQ=255)
                                                                          tCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCtt                                                                    <  2:163176/150‑1 (MQ=255)
                                                                            gCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTa                                                                  >  1:629886/1‑150 (MQ=255)
                                                                              tGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACg                                                                <  1:934668/150‑1 (MQ=255)
                                                                                cgcCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCCGCTCAC‑CATCCAATGGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGt                                                              <  2:897229/150‑1 (MQ=255)
                                                                                    cTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GACTTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCtt                                                          <  2:1008348/150‑1 (MQ=255)
                                                                                         gACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGttttt                                                                                                 >  2:993965/1‑106 (MQ=255)
                                                                                         gACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGttttt                                                                                                 <  1:993965/106‑1 (MQ=255)
                                                                                              ccAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTg                                                >  1:568523/1‑150 (MQ=255)
                                                                                                aGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGcc                                              <  1:477743/150‑1 (MQ=255)
                                                                                                 gATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCg                                             <  2:629886/150‑1 (MQ=255)
                                                                                                             ccGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGAc                                                                                                                       >  2:120450/1‑65 (MQ=255)
                                                                                                             ccGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGAc                                                                                                                       <  1:120450/65‑1 (MQ=255)
                                                                                                             ccGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCa                                 <  2:205390/150‑1 (MQ=255)
                                                                                                                 tttGGATGGTGCTGCAGCAATGCCTGCGCGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCtt                             <  1:318370/150‑1 (MQ=255)
                                                                                                                 tttGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCtt                             <  2:73210/150‑1 (MQ=255)
                                                                                                                     gATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGtt                                                                                                    >  1:68094/1‑75 (MQ=255)
                                                                                                                     gATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGtt                                                                                                    <  2:68094/75‑1 (MQ=255)
                                                                                                                                agcaATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTtgctgc              >  1:165478/1‑150 (MQ=255)
                                                                                                                                       ccTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCAct       >  1:369102/1‑150 (MQ=255)
                                                                                                                                         tctaCGGCGTACACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGcc                                              <  2:415230/106‑1 (MQ=255)
                                                                                                                                         tctaCGGCGTACACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGcc                                              >  1:415230/4‑109 (MQ=255)
                                                                                                                                            aCGGCGTGTACGCCGTAGATTTTTTCCAGCTGAATCAT‑GCAATGCTCTTGCTGTTGCGTGGAGAGCCGGACTTCGGCGGGCCTTTACAGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCACtgctgc  >  1:626931/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TCTCACGGCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAT‑GAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCACTGCTGC  >  NC_002947/5546354‑5546643

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