Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,327,261 A→G E137G (GAA→GGA)  rluB → 23S rRNA pseudouridine(2605) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,327,2610AG98.6% 248.3 / ‑4.8 70E137G (GAA→GGA) rluB23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (33/36);  total (33/37)
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.67e-01

TTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCA  >  NC_000913/1327126‑1327396
                                                                                                                                       |                                                                                                                                       
ttGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACg                                                                                                                                      <  2:1557835/139‑1 (MQ=255)
       tCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAatat                                                                                                                               >  2:1625759/1‑139 (MQ=255)
         tGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGc                                                                                                                             >  1:1222699/1‑139 (MQ=255)
          gCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGcc                                                                                                                            <  2:146879/139‑1 (MQ=255)
                    cGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTAt                                                                                                                  >  1:391711/1‑139 (MQ=255)
                         cgcACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGt                                                                                                             <  1:357498/139‑1 (MQ=255)
                                       ccGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGAcgc                                                                                               >  2:763033/1‑139 (MQ=255)
                                            ggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAAc                                                                                          >  1:12083/1‑139 (MQ=255)
                                            ggTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAAc                                                                                          >  2:708717/1‑139 (MQ=255)
                                               cGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGc                                                                                       <  1:708717/139‑1 (MQ=255)
                                                 cTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCcgtg                                                                                                                         <  2:533340/103‑1 (MQ=255)
                                                 cTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCcgtg                                                                                                                         >  1:533340/1‑103 (MQ=255)
                                                 cTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGt                                                                                     >  1:372776/1‑139 (MQ=255)
                                                            aTACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCg                                                                                      <  2:99614/127‑1 (MQ=255)
                                                            aTACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCg                                                                                      >  1:99614/1‑127 (MQ=255)
                                                                  gtgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGc                                                                                       <  1:1600772/120‑1 (MQ=255)
                                                                  gtgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGc                                                                                       >  2:1600772/1‑120 (MQ=255)
                                                                  gtgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGt                                                                    <  2:12083/139‑1 (MQ=255)
                                                                  gtgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGt                                                                    <  2:1222699/139‑1 (MQ=255)
                                                                  gtgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGt                                                                    <  1:250854/139‑1 (MQ=255)
                                                                   tgGTCTGCTGCTGTTGACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTg                                                                   <  1:249180/139‑1 (MQ=255)
                                                                   tgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTg                                                                   <  1:556791/139‑1 (MQ=255)
                                                                   tgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTg                                                                   <  2:566666/139‑1 (MQ=255)
                                                                   tgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTg                                                                   >  2:556791/1‑139 (MQ=255)
                                                                   tgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTg                                                                   <  2:372776/139‑1 (MQ=255)
                                                                   tgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTg                                                                   <  1:414808/139‑1 (MQ=255)
                                                                   tgGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTg                                                                   >  2:822/1‑139 (MQ=255)
                                                                    ggTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGcc                                                                                                                            <  1:75338/81‑1 (MQ=255)
                                                                    ggTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGcc                                                                                                                            >  2:75338/1‑81 (MQ=255)
                                                                     gTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAata                                                                                                                                <  1:1564999/76‑1 (MQ=255)
                                                                     gTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAata                                                                                                                                >  2:1564999/1‑76 (MQ=255)
                                                                        tgctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGtt                                                              <  1:1625759/139‑1 (MQ=255)
                                                                          ctgctgTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTgg                                                            >  2:376259/1‑139 (MQ=255)
                                                                                       ccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACg                                                                                                                                      >  1:60695/1‑52 (MQ=255)
                                                                                       ccGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACg                                                                                                                                      <  2:60695/52‑1 (MQ=255)
                                                                                            ggTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGc                                                                                                                             <  1:358411/56‑1 (MQ=255)
                                                                                            ggTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGc                                                                                                                             >  2:358411/1‑56 (MQ=255)
                                                                                            ggTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTCGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGTTCCGGCAGCtt                                          >  2:1091655/1‑139 (MQ=255)
                                                                                             gTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGt                                                                                                             <  1:56181/71‑1 (MQ=255)
                                                                                             gTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGt                                                                                                             >  2:56181/1‑71 (MQ=255)
                                                                                                 aCTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTaaa                                     >  1:1087350/1‑139 (MQ=255)
                                                                                                    ggCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGaaa                                                                                           <  1:41559/82‑1 (MQ=255)
                                                                                                    ggCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGaaa                                                                                           >  2:41559/1‑82 (MQ=255)
                                                                                                     gCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGa                                                                                                                                   >  2:794248/1‑41 (MQ=255)
                                                                                                     gCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGa                                                                                                                                   <  1:794248/41‑1 (MQ=255)
                                                                                                     gCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTc                                                                           <  2:957610/97‑1 (MQ=255)
                                                                                                     gCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTc                                                                           >  1:957610/1‑97 (MQ=255)
                                                                                                      cAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCgg                                                <  2:878795/123‑1 (MQ=255)
                                                                                                      cAAACCGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCgg                                                >  1:878795/1‑123 (MQ=255)
                                                                                                          ccGTTTAATGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAg                            >  2:675892/1‑139 (MQ=255)
                                                                                                                  tGCACCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAgcgg                    <  2:1087350/139‑1 (MQ=255)
                                                                                                                     aCCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGt                                                                                                             >  2:1257076/1‑47 (MQ=255)
                                                                                                                     aCCCAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGt                                                                                                             <  1:1257076/47‑1 (MQ=255)
                                                                                                                       ccAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGaaa                                                                                           >  2:1584921/1‑63 (MQ=255)
                                                                                                                       ccAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGaaa                                                                                           <  1:1584921/63‑1 (MQ=255)
                                                                                                                       ccAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCAt                                <  1:1621080/122‑1 (MQ=255)
                                                                                                                       ccAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCAt                                >  2:1621080/1‑122 (MQ=255)
                                                                                                                       ccAAGCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAgcggcggcg               >  2:29196/1‑139 (MQ=255)
                                                                                                                           gCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAg                            <  2:438853/122‑1 (MQ=255)
                                                                                                                           gCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAg                            >  1:438853/1‑122 (MQ=255)
                                                                                                                           gCCGTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAgg           <  1:675892/139‑1 (MQ=255)
                                                                                                                              gTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTc                                                                           >  2:946246/1‑72 (MQ=255)
                                                                                                                              gTGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTc                                                                           <  1:946246/72‑1 (MQ=255)
                                                                                                                               tGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGt                                                                    >  2:1319429/1‑78 (MQ=255)
                                                                                                                               tGAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGt                                                                    <  1:1319429/78‑1 (MQ=255)
                                                                                                                                gAAGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCa      <  2:391711/139‑1 (MQ=255)
                                                                                                                                 aaGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGc                                                                                       >  2:589081/1‑57 (MQ=255)
                                                                                                                                 aaGTTGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGc                                                                                       <  1:589081/57‑1 (MQ=255)
                                                                                                                                    ttGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCa  <  1:376259/139‑1 (MQ=255)
                                                                                                                                    ttGGACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCa  <  1:822/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TTGACCGTCTGCCAAAACTGCGTGGCGCACGCTGGATTGCCGTGGGTCGCTTGGACGTTAATACCTGTGGTCTGCTGCTGTTCACCACCGATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCA  >  NC_000913/1327126‑1327396

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑