Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,327,340 A→G E163E (GAA→GAG 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,3400AG98.2% 195.2 / ‑4.6 56E163E (GAA→GAGrluB23S rRNA pseudouridine(2605) synthase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (30/25);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCGGTTGGTGTGCAGGTGA  >  NC_000913/1327216‑1327473
                                                                                                                            |                                                                                                                                     
aTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCt                                                                                                                        >  2:9786/1‑140 (MQ=255)
aTGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCt                                                                                                                        >  1:831994/1‑140 (MQ=255)
    tGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTa                                                                                                                    >  2:482109/1‑140 (MQ=255)
     gAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTaa                                                                                                                   >  1:773361/1‑140 (MQ=255)
       aCTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTaaaa                                                                                                                 >  1:270757/1‑140 (MQ=255)
           gCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCAt                                                                                                             >  1:863971/1‑140 (MQ=255)
              aaCCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCaa                                                                                                          >  2:103122/1‑140 (MQ=255)
                    ttAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAg                                                                                                    <  1:621540/140‑1 (MQ=255)
                       aTGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAgcgg                                                                                                 <  2:773361/140‑1 (MQ=255)
                       aTGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAgcgg                                                                                                 <  1:482109/140‑1 (MQ=255)
                           aCCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAgcggcggc                                                                                             >  1:562829/1‑140 (MQ=255)
                           aCCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAgcggcggc                                                                                             >  2:67952/1‑140 (MQ=255)
                           aCCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAgcggcggc                                                                                             >  2:67951/1‑140 (MQ=255)
                            cccAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAgcggcggcg                                                                                            >  2:309087/1‑140 (MQ=255)
                             ccAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGa                                                                                           <  1:49951/140‑1 (MQ=255)
                                     tGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCa                                                                                   <  2:562829/140‑1 (MQ=255)
                                     tGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCa                                                                                   <  1:67951/140‑1 (MQ=255)
                                     tGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCa                                                                                   <  1:67952/140‑1 (MQ=255)
                                       aaGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAAc                                                                                 <  1:42469/140‑1 (MQ=255)
                                       aaGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAAc                                                                                 <  2:831994/140‑1 (MQ=255)
                                         gTTGAACGTGAATATGCTGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCa                                                                               <  2:863971/140‑1 (MQ=255)
                                            gAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTg                                                                            >  2:489627/1‑140 (MQ=255)
                                                   aaTATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGTGGCGGCGAAGGGATCAACCAGTGGTa                                                                         >  2:595679/1‑136 (MQ=255)
                                                   aaTATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTa                                                                         <  1:595679/136‑1 (MQ=255)
                                                        gCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGAc                                                                >  2:479930/1‑140 (MQ=255)
                                                        gCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGAc                                                                >  2:670963/1‑140 (MQ=255)
                                                          cgtgcgtgTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGActc                                                              >  2:278741/1‑140 (MQ=255)
                                                            tgcgtgTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTg                                                            >  1:9276/1‑140 (MQ=255)
                                                             gcgtgTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAGCGTGACTCTGa                                                           >  1:305036/1‑140 (MQ=255)
                                                                   aTTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAg                                                     >  1:470767/1‑140 (MQ=255)
                                                                        gTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGt                                                >  1:534151/1‑140 (MQ=255)
                                                                         tCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTa                                               >  1:243130/1‑140 (MQ=255)
                                                                             gTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCg                                           >  2:95034/1‑140 (MQ=255)
                                                                              ttGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGaa                                                      >  2:761537/1‑128 (MQ=255)
                                                                              ttGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGaa                                                      <  1:761537/128‑1 (MQ=255)
                                                                               tGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTg                                         <  1:479930/140‑1 (MQ=255)
                                                                                gacgacGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGa                                        <  1:103122/140‑1 (MQ=255)
                                                                                gacgacGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGa                                        <  1:95034/140‑1 (MQ=255)
                                                                                gacgacGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGa                                        <  2:431179/140‑1 (MQ=255)
                                                                                gacgacGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGa                                        <  1:309087/140‑1 (MQ=255)
                                                                                 acgacgCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGaa                                       >  2:392848/1‑140 (MQ=255)
                                                                                       cgAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTtcgt                                 >  1:259456/1‑140 (MQ=255)
                                                                                               cGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCtgtg                         >  1:114268/1‑140 (MQ=255)
                                                                                               cGTGATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCtgtg                         >  2:66383/1‑140 (MQ=255)
                                                                                                  gATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGa                      <  1:220349/140‑1 (MQ=255)
                                                                                                  gATTTGAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGa                      <  2:243130/140‑1 (MQ=255)
                                                                                                  gATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGa                      <  2:624876/140‑1 (MQ=255)
                                                                                                       gAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGGCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCgg                 >  2:888465/1‑140 (MQ=255)
                                                                                                       gAGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCgg                 >  1:657699/1‑140 (MQ=255)
                                                                                                        aGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCGGt                <  2:470767/140‑1 (MQ=255)
                                                                                                        aGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCGGt                <  2:270757/140‑1 (MQ=255)
                                                                                                        aGTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCGGt                <  1:9786/140‑1 (MQ=255)
                                                                                                         gTCGTGGCGTGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCGGtt               <  2:259456/140‑1 (MQ=255)
                                                                                                                  tGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCGGTTg              >  2:631653/1‑132 (MQ=255)
                                                                                                                  tGCAGTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCGGTTg              <  1:631653/132‑1 (MQ=255)
                                                                                                                      gTTGGAGGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCGGTTGGTGTGCAGGTGa  <  2:657699/140‑1 (MQ=255)
                                                                                                                            |                                                                                                                                     
ATGGTGAACTGGCAAACCGTTTAATGCACCCAAGCCGTGAAGTTGAACGTGAATATGCCGTGCGTGTATTTGGTCAGGTTGACGACGCGAAACTGCGTGATTTGAGTCGTGGCGTGCAGTTGGAAGATGGTCCGGCAGCTTTTAAAACCATCAAGTTCAGCGGCGGCGAAGGGATCAACCAGTGGTACAACGTGACTCTGACCGAAGGCCGTAACCGTGAAGTTCGTCGTCTGTGGGAAGCGGTTGGTGTGCAGGTGA  >  NC_000913/1327216‑1327473

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