Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,410,053 (T)9→8 100% intergenic (+47/‑80) rlmB → / → yjfI 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,410,0450T.87.5% 212.3 / 31.0 64intergenic (+39/‑88)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (5/3);  new base . (30/26);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATG  >  NC_000913/4409947‑4410140
                                                                                                  |                                                                                               
cTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtt                                                                                                >  2:996786/1‑100 (MQ=255)
  gAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTTTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttt                                                                                               <  2:154496/99‑1 (MQ=255)
   aaCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttt                                                                                             <  2:113963/100‑1 (MQ=255)
     cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                           >  2:2161120/1‑100 (MQ=255)
     cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                           <  1:1224386/100‑1 (MQ=255)
     cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                           <  1:968212/100‑1 (MQ=255)
      gTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcat                                                                                          <  1:1082384/100‑4 (MQ=37)
       tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                         <  2:835641/100‑2 (MQ=255)
       tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                         <  2:1862157/100‑2 (MQ=255)
       tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                         >  1:1974726/1‑99 (MQ=255)
       tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                         <  1:29678/100‑2 (MQ=255)
         tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                       <  1:996791/100‑4 (MQ=255)
         tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                       >  2:1647639/1‑97 (MQ=255)
         tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                       <  1:1168882/100‑4 (MQ=255)
          cGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtct                                                                                      >  1:1439028/1‑96 (MQ=255)
            gTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCtat                                                                                   >  2:1140919/1‑100 (MQ=255)
            gTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCtat                                                                                   <  1:473886/100‑1 (MQ=255)
                cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                               <  1:2252807/100‑1 (MQ=255)
                cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                               >  1:2023903/1‑100 (MQ=255)
                   ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTTTATAGc                                                                                >  2:1663448/1‑97 (MQ=255)
                    cGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTCTATAGCTGTTa                                                                           >  2:1031499/1‑100 (MQ=255)
                     ggAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTaa                                                                          >  2:211948/1‑100 (MQ=255)
                           tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGCTGGTTTTTTTTGTATATCGCTGGAAAGCTaa                                                                     >  1:2199898/1‑99 (MQ=255)
                             cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTCTATATTTTGTAAGTTAAtt                                                                   >  2:1818481/1‑99 (MQ=255)
                             cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                  >  1:1989252/1‑100 (MQ=255)
                             cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                  >  1:1953219/1‑100 (MQ=255)
                             cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                  <  1:790103/100‑1 (MQ=255)
                              ttATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTAc                                                                 <  2:894584/100‑1 (MQ=255)
                               tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                                >  2:665094/1‑100 (MQ=255)
                               tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                                <  1:1030836/100‑1 (MQ=255)
                               tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCAGGTAAGATAATTACg                                                                >  1:2321603/1‑100 (MQ=255)
                                   tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTCTATAGCTGGTAAGATAATTACGTAtt                                                            >  2:1752567/1‑100 (MQ=255)
                                   tGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAtt                                                            >  2:2023681/1‑100 (MQ=255)
                                    gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATTCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTg                                                           <  2:268821/100‑1 (MQ=255)
                                    gAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTa                                                              >  2:159145/1‑97 (MQ=255)
                                              gcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCa                                                 <  2:499652/100‑1 (MQ=255)
                                                cAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                               >  1:1307035/1‑100 (MQ=255)
                                                cAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAtt                                               <  2:1485140/100‑1 (MQ=255)
                                                 aGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCAttt                                              >  1:1483595/1‑100 (MQ=255)
                                                    gcAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttat                                           <  2:901518/100‑1 (MQ=255)
                                                     cAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTttatt                                          <  1:1293245/100‑1 (MQ=255)
                                                       gCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTaa                                        >  2:213046/1‑100 (MQ=255)
                                                        cTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAat                                       >  1:2270717/1‑100 (MQ=255)
                                                             ttCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTa                                  <  1:729863/100‑1 (MQ=255)
                                                                tcATAATTATGAAAAGCCTTCCAGATTTGGATGGTTTTTTTTTTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                               <  2:616878/100‑1 (MQ=255)
                                                                 tatagacaaaaaaaACCATCCAAATCTGGATGGCTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                              >  1:2256657/11‑100 (MQ=21)
                                                                 cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATATCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                              >  2:939521/1‑100 (MQ=255)
                                                                 cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                              >  2:938334/1‑100 (MQ=255)
                                                                 cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                              >  2:708573/1‑100 (MQ=255)
                                                                 cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAAGATTTAAACt                              <  1:1928216/100‑1 (MQ=255)
                                                                   taaTTATGAAAAGCCATTCAGATTTGGATAGTTATTTTTGTTTATATATGTTAAGAGGATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                            <  2:153879/99‑1 (MQ=255)
                                                                   gAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTa                            >  2:1282225/1‑100 (MQ=255)
                                                                      ttATGAAAGGCCATCCAGATTTGGATGGTTTTTTTTGTTTATAGCTGTGAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                         <  2:968734/100‑1 (MQ=255)
                                                                      ttATGAAAAGCCCTCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                         <  1:2072305/100‑1 (MQ=255)
                                                                      ttATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                         >  2:2334335/1‑100 (MQ=255)
                                                                      ttATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                         <  2:34124/100‑1 (MQ=255)
                                                                       tATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCt                        <  2:1659032/100‑1 (MQ=255)
                                                                       caaaaaaaACCATCCAAATCTGGATGGCTTTTTTTGTCTATAGCTGGTAAAATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATc                         >  2:1404101/5‑99 (MQ=25)
                                                                        aTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTTCTATAGTTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                       <  2:2206951/100‑1 (MQ=255)
                                                                        aTGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                       <  1:1677765/100‑1 (MQ=255)
                                                                         aaaaaaACCATCCAAATCTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtat                      >  1:510897/3‑100 (MQ=38)
                                                                          gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                     >  2:1193973/1‑100 (MQ=255)
                                                                           aaaaGCCATCCAGATTTGGATGGTTTTTTTTTTTTTTAGTTTTTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAATTTTATTTTatat                    >  2:1076260/1‑100 (MQ=255)
                                                                              aGCTATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                   <  2:62138/98‑1 (MQ=255)
                                                                               gCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa                >  1:481238/1‑100 (MQ=255)
                                                                               gCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaa                >  1:1289274/1‑100 (MQ=255)
                                                                                ccATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat               <  1:1647650/100‑1 (MQ=255)
                                                                                ccATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaat               <  2:1113855/100‑1 (MQ=255)
                                                                                 cATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaata              >  1:440757/1‑100 (MQ=255)
                                                                                   tCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGaataat            >  1:2306343/1‑100 (MQ=255)
                                                                                    ccAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata           >  1:1652220/1‑100 (MQ=255)
                                                                                    ccAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata           >  2:1501351/1‑100 (MQ=255)
                                                                                      aGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata         <  2:630071/100‑1 (MQ=255)
                                                                                      aGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata         <  2:764075/100‑1 (MQ=255)
                                                                                       gATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata           <  2:2281010/97‑1 (MQ=255)
                                                                                          tttcctagTTTTTTTTTTTATTTCGTTGGTAAGATAATTAATTATTCTAAACATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg       <  2:2120201/92‑1 (MQ=255)
                                                                                          ttGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGAc     <  1:944928/100‑1 (MQ=255)
                                                                                          ttGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGAc     <  1:71703/100‑1 (MQ=255)
                                                                                             gATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATg  >  2:1331100/1‑100 (MQ=255)
                                                                                                  |                                                                                               
CTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACATG  >  NC_000913/4409947‑4410140

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