Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,389,141 T→C 7.4% M1T (ATG→ACG) † 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,1410TC 7.4% 158.0 / 2.4 54M1T (ATG→ACG) †rlmKL23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base C (2/2);  total (26/28)
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.42e-01

ATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCG  >  NC_002947/2389005‑2389281
                                                                                                                                         |                                                                                                                                            
atgatCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgt                                                                                                                                  <  4:307915/149‑1 (MQ=255)
   atCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCATTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgt                                                                                                                                  >  1:84928/1‑147 (MQ=255)
   atCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCACGGGGCGACGCATCGCCCCgcga                                                                                                                                <  3:311491/148‑4 (MQ=255)
    tCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcagt                                                                                                                              >  3:138386/1‑149 (MQ=255)
     cTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcagt                                                                                                                              >  6:102243/1‑148 (MQ=255)
       aaGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcagtca                                                                                                                            >  8:242494/1‑148 (MQ=255)
        aGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCaaa                                                                                                                          <  5:149193/149‑1 (MQ=255)
               cTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGccc                                                                                                                    <  8:111985/148‑1 (MQ=255)
                         ggCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcag                                                                                                                               <  5:42692/127‑1 (MQ=255)
                         ggCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcag                                                                                                                               >  6:42692/1‑127 (MQ=255)
                         ggCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTg                                                                                                           <  4:138386/147‑1 (MQ=255)
                              ggCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCTAATAGCCCAAGGCGCTGGAATACt                                                                                                    <  2:84928/149‑1 (MQ=255)
                              ggCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACt                                                                                                    <  7:242494/149‑1 (MQ=255)
                              ggCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACt                                                                                                    <  4:96941/149‑1 (MQ=255)
                               gCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACt                                                                                                    >  3:79409/1‑148 (MQ=255)
                                        ggTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGAc                                                                                          <  5:102243/149‑1 (MQ=255)
                                         gTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGGCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGAcc                                                                                         >  1:98644/1‑149 (MQ=255)
                                          tCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGc                                                                                                                      <  4:18549/119‑1 (MQ=255)
                                            aGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGtt                                                                                      <  6:173085/149‑1 (MQ=255)
                                                  cccGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAAc                                                                                >  4:41232/1‑149 (MQ=255)
                                                   ccGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACt                                                                               >  5:67045/1‑149 (MQ=255)
                                                   ccGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACt                                                                               >  7:209906/1‑149 (MQ=255)
                                                    cGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACtt                                                                              >  6:289135/1‑149 (MQ=255)
                                                     ggCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACttt                                                                             <  3:41232/149‑1 (MQ=255)
                                                     ggCTTTTCTGCAGCAAGTTGTGTGTCGCAGAAGGCCACCAGCACGACTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACtt                                                                              >  4:273411/1‑148 (MQ=255)
                                                                cagcaCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTg                                                                   >  2:10353/1‑149 (MQ=255)
                                                                   caCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCt                                                                                                            <  2:83615/105‑1 (MQ=255)
                                                                   caCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTcc                                                                <  8:209906/149‑1 (MQ=255)
                                                                   caCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGAAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCt                                                                                                            >  1:83615/1‑105 (MQ=255)
                                                                      gTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCa                                                              >  5:80731/1‑148 (MQ=255)
                                                                       ttGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCAt                                                                                                 <  1:284045/112‑1 (MQ=255)
                                                                       ttGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCAt                                                                                                 >  2:284045/1‑112 (MQ=255)
                                                                              tCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCtt                                                      >  6:36529/1‑148 (MQ=255)
                                                                                  aGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCCCCTGTCCCAAAGGCCTTGAAAg                                                 >  6:14011/1‑149 (MQ=255)
                                                                                  aGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCCTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACCCATGTCGGACCGTTTAGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAg                                                 >  7:230810/1‑149 (MQ=255)
                                                                                       gCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGc                                            >  1:18115/1‑149 (MQ=255)
                                                                                            ccGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGcc                                       <  5:289135/149‑1 (MQ=255)
                                                                                                aCGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCg                                      >  4:288589/1‑146 (MQ=255)
                                                                                                          gAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                          <  5:14011/148‑1 (MQ=255)
                                                                                                          gAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                          <  2:98644/148‑1 (MQ=255)
                                                                                                          gAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                          <  6:67045/148‑1 (MQ=255)
                                                                                                          gAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAg                           <  8:21838/147‑1 (MQ=255)
                                                                                                               gtgtACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCg                    <  6:80731/149‑1 (MQ=255)
                                                                                                                  tACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                          <  3:134864/140‑1 (MQ=255)
                                                                                                                  tACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                          >  4:134864/1‑140 (MQ=255)
                                                                                                                      ggCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGa             <  3:273411/149‑1 (MQ=255)
                                                                                                                            cATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGt       <  3:288589/149‑1 (MQ=255)
                                                                                                                              tgacgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAagag                                 >  6:208320/4‑121 (MQ=255)
                                                                                                                              tgacgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAagag                                 <  5:208320/118‑1 (MQ=255)
                                                                                                                               gacgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTgcg    >  2:181998/3‑149 (MQ=255)
                                                                                                                                gcgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTgcgc   >  6:99058/1‑149 (MQ=255)
                                                                                                                                gcgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTgcgc   <  2:18115/149‑1 (MQ=255)
                                                                                                                                gcgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTgcgc   >  3:285611/1‑149 (MQ=255)
                                                                                                                                 cgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGcgcg  <  4:79409/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
ATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCG  >  NC_002947/2389005‑2389281

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