Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 60,316 C→T E11K (GAA→AAA)  rluA ← dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091360,3160CT98.5% 225.0 / ‑4.6 68E11K (GAA→AAA) rluAdual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (33/34);  total (33/35)
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.65e-01

TGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCAACCTGCC  >  NC_000913/60220‑60409
                                                                                                |                                                                                             
tgtgCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGTATAACCAACCAGGGTTTCTg                                                                                            <  1:304274/100‑1 (MQ=255)
 gtgCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGc                                                                                           >  1:1157378/1‑100 (MQ=255)
 gtgCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGc                                                                                           <  2:459032/100‑1 (MQ=255)
  tgCTCTTCCAGACTACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCg                                                                                          <  1:1682005/100‑1 (MQ=255)
   gCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCgg                                                                                         >  2:199586/1‑100 (MQ=255)
     tctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTg                                                                                       <  1:548848/100‑1 (MQ=255)
     tctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTg                                                                                       >  2:1089106/1‑100 (MQ=255)
      ctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTgg                                                                                      <  2:1103787/100‑1 (MQ=255)
      ctTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGAGGTgg                                                                                      <  1:779572/100‑1 (MQ=255)
        tCCAGACGACCCGGCACTGACAACAAGCCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGAt                                                                                    >  2:1108432/1‑100 (MQ=255)
         ccAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCAGTGGAtt                                                                                   >  2:1587395/1‑100 (MQ=255)
          cAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTg                                                                                  <  2:316431/100‑1 (MQ=255)
           aGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGt                                                                                 <  2:107740/100‑1 (MQ=255)
             acgacCCGGTACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAg                                                                               <  2:1432152/100‑1 (MQ=255)
                 cccGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGtt                                                                             <  1:1698815/98‑1 (MQ=255)
                    ggCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCa                                                                        <  1:469837/100‑1 (MQ=255)
                       aCTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATcc                                                                     <  1:591317/100‑1 (MQ=255)
                       aCTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATcc                                                                     <  1:972910/100‑1 (MQ=255)
                         tGACAACAAACCGCTCGGCTTGATGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATcctc                                                                   >  1:1056983/1‑98 (MQ=255)
                          gACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCa                                                                  >  1:1360352/1‑100 (MQ=255)
                          gACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCa                                                                  >  2:86618/1‑100 (MQ=255)
                           acaacaAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCAt                                                                 >  1:1500111/1‑100 (MQ=255)
                            caacaaACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCAtt                                                                >  2:1150561/1‑100 (MQ=255)
                            caacaaACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTACGGTGGATTGTAGTTTTCCATCCCCAtt                                                                >  2:1014897/1‑100 (MQ=255)
                                aaaCCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCg                                                            >  1:202460/1‑100 (MQ=255)
                                aaaCCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCg                                                            >  1:670603/1‑100 (MQ=255)
                                    cGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCt                                                         <  1:548077/99‑1 (MQ=255)
                                    cGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTc                                                        >  2:200409/1‑100 (MQ=255)
                                    cGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTc                                                        >  1:637468/1‑100 (MQ=255)
                                    cGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTc                                                        <  1:881265/100‑1 (MQ=255)
                                     gCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCAATCCCCATTTTCGGCTcc                                                       >  2:429674/1‑100 (MQ=255)
                                      cTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGc                                                          <  1:892135/96‑1 (MQ=255)
                                      cTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCg                                                      >  2:950980/1‑100 (MQ=255)
                                       tCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGt                                                     >  2:1802883/1‑100 (MQ=255)
                                        cGGCTTGTAGACCACCATAATATGGTCATCCTGATAAAGGATAACCAAACAGAGGTTCGACGACGGGTTGTAGCATTCCATCCCCAATCTCGGCTACggt                                                    >  1:1401161/1‑98 (MQ=255)
                                          gCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTAc                                                  <  2:1185265/100‑1 (MQ=255)
                                            ttgttgACCACCATAATATGTTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTg                                                <  2:1693470/100‑1 (MQ=255)
                                            ttgttgACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTg                                                <  2:405351/100‑1 (MQ=255)
                                              gttgACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGAt                                              <  2:476658/100‑1 (MQ=255)
                                              gttgACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCAACTCAGTTACTGa                                               >  2:278695/1‑99 (MQ=255)
                                                ngACCACCATAATATGGTCATCCAGATACAGGATAACCAACCAGGGTTTCTGCGGTGGCTAGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGAt                                              >  1:393462/2‑98 (MQ=255)
                                                 gACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCg                                           <  2:480110/100‑1 (MQ=255)
                                                   ccaccaTAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGtt                                         >  1:240077/1‑100 (MQ=255)
                                                     accaTAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGTGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTAc                                       >  1:504810/1‑100 (MQ=255)
                                                     accaTAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTAc                                       >  1:610418/1‑100 (MQ=255)
                                                        ataataTGGTCATCCTGTTACAGGATAACCAACCAGGGTTCCTGCGGTGGAGTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACaa                                     <  2:1148593/99‑1 (MQ=255)
                                                         taataTGTTCATCCTGATACATGATAACCAACCAGTTTGTCTTCGGTTGATTGTATTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAAc                                    <  2:1058397/99‑1 (MQ=255)
                                                            tatGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATc                                <  1:950973/100‑1 (MQ=255)
                                                             atGGTTATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTAGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGaac                                <  1:312025/99‑3 (MQ=255)
                                                             atGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATc                                <  2:702284/99‑1 (MQ=255)
                                                              tGGTCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCaa                              >  2:1296836/1‑100 (MQ=255)
                                                                 tCATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACg                           <  1:133841/100‑1 (MQ=255)
                                                                  cATCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGc                          >  2:867652/1‑100 (MQ=255)
                                                                    tCCTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAg                        <  1:421872/100‑1 (MQ=255)
                                                                      cTGATACAGTATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAggg                      <  2:528593/100‑1 (MQ=255)
                                                                      cTGATACAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAggg                      >  2:1351414/1‑100 (MQ=255)
                                                                        gATACAGGATAACCAACCTGGGTTTCTCCCGTCAATTGTAATTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAggg                      <  1:1142651/98‑1 (MQ=255)
                                                                            cAGGATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCa                <  1:1218895/100‑1 (MQ=255)
                                                                              ggATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGa              >  1:846481/1‑100 (MQ=255)
                                                                              ggATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGa              >  2:1506207/1‑100 (MQ=255)
                                                                               gATAACCAACCAGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGAc             >  2:1050948/1‑100 (MQ=255)
                                                                                    ccaaccaGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCaa        >  1:340604/1‑100 (MQ=255)
                                                                                    ccaaccaGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGAGCATACAGACGCCaa        >  2:1665173/1‑100 (MQ=255)
                                                                                     caaccaGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCAAc       <  1:1782148/100‑1 (MQ=255)
                                                                                     caaccaGCTTTTCTGCGGTGGATTCTAGTGTTCCATCCCCATTTGCGGCTCCGTTACTGATGCTTTACAACGATCAAACGCAGGGCATCCAGACGCCa         <  1:1658501/98‑1 (MQ=255)
                                                                                       accaGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCAAcct     <  2:1413976/100‑1 (MQ=255)
                                                                                        ccaGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCAAcctg    >  2:420542/1‑100 (MQ=255)
                                                                                          aGGGTTTCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCAAcctgcc  <  1:1405766/100‑1 (MQ=255)
                                                                                                |                                                                                             
TGTGCTCTTCCAGACGACCCGGCACTGACAACAAACCGCTCGGCTTGTTGACCACCATAATATGGTCATCCTGATACAGGATAACCAACCAGGGTTCCTGCGGTGGATTGTAGTTTTCCATCCCCATTTTCGGCTCCGTTACTGATGCGTTACAACGATCAAACGCAGGGCATCCAGACGCCAACCTGCC  >  NC_000913/60220‑60409

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