Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 2,145,263 A→G 7.0% E297G (GAA→GGA)  rluC → Ribosomal large subunit pseudouridine synthaseC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,145,2630AG 7.0% 176.3 / 3.8 57E297G (GAA→GGA) rluCRibosomal large subunit pseudouridine synthaseC
Reads supporting (aligned to +/- strand):  ref base A (25/28);  new base G (2/2);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

ATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGCTACTGATCTTCGACTGGGACGGCACCCTAGCCGATTCCATTGGG  >  KT2440/2145131‑2145396
                                                                                                                                    |                                                                                                                                     
aTGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAAt                                                                                                                                     >  8:249969/1‑135 (MQ=255)
  gCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAAtt                                                                                                                                    >  8:261013/1‑134 (MQ=255)
      ggCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGc                                                                                                                               >  8:148230/1‑135 (MQ=255)
        cacaTGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTTACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCt                                                                                                                              >  2:265559/1‑134 (MQ=255)
        cacaTGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGGATTGAAGCt                                                                                                                              >  5:25712/1‑134 (MQ=255)
        cacaTGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGGATTGAAGCTc                                                                                                                             >  5:20924/1‑135 (MQ=255)
        cacaTGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTc                                                                                                                             >  2:83860/1‑135 (MQ=255)
          caTGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGa                                                                                                                           >  8:319370/1‑135 (MQ=255)
             gagcGCTGGAGACATCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGgc                                                                                                                        <  8:142322/132‑1 (MQ=255)
              aTCGCGGGCGACAGCAAATACGGCGACGACGACTTCAGCCGGGAAATTCGCGAGTTGGGCTGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTTAATTGAAGCTCGAGgcgc                                                                                                                      <  3:156893/136‑1 (MQ=255)
               tCGCTGTCTACAGCAAATACGGCGACGAAGAATTCACCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCAGGACGGCGGTGAATTGAAGCTCGAGgcgc                                                                                                                      <  2:179043/135‑1 (MQ=255)
               tCGCTGGCGACACCAATTACGGCGACGAAGACTTCAGCCGGTAACTTCGCTAGTTGGGCGGCAAACGCCGGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGgcgc                                                                                                                      <  1:191795/135‑1 (MQ=255)
                cGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTTACTGTACCGCTGCCGGACGGCGGTTAATTGAAGCTCGAGgcgc                                                                                                                      <  1:265559/134‑1 (MQ=255)
                     gCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTgg                                                                                                                >  5:115302/1‑135 (MQ=255)
                      cGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCACGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTgg                                                                                                                <  7:148230/134‑1 (MQ=255)
                        aCAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGa                                                                                                               >  7:368758/1‑133 (MQ=255)
                        aCAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGAGAGTTGAGCGGCAAACGCCTGTTCCTGCATGACTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATg                                                                                                             >  8:77188/1‑135 (MQ=255)
                          aGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGGCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAg                                                                                                           <  1:14849/135‑1 (MQ=255)
                          aGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAg                                                                                                           <  7:249969/135‑1 (MQ=255)
                           gCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGt                                                                                                         >  3:155848/1‑136 (MQ=255)
                                   ggCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGc                                                                                                   <  8:368758/134‑1 (MQ=255)
                                      gacgaAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCaaaa                                                                                               >  2:368412/1‑135 (MQ=255)
                                        cgaAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTCCCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAAcc                                                                                            <  3:166219/136‑1 (MQ=255)
                                           aGACTTCAGCCGGGACATTGGCGAGGTGGGCGGCTAACGCCTGTTCCTGCAGGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGt                                                                                          <  2:192759/135‑1 (MQ=255)
                                           aGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGGATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGt                                                                                          <  6:25712/135‑1 (MQ=255)
                                           aGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGt                                                                                          <  2:280372/135‑1 (MQ=255)
                                            gACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGGATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGt                                                                                          <  6:20924/134‑1 (MQ=255)
                                             aCTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGTATGAGGTCTGGGCAAAAACCGtt                                                                                         >  1:365030/1‑134 (MQ=255)
                                                  aGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAAc                                                                                     <  1:368412/133‑1 (MQ=255)
                                                  aGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGATCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAAc                                                                                     <  1:40388/133‑1 (MQ=255)
                                                    ccGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTc                                                                                  >  5:195904/1‑134 (MQ=255)
                                                    ccGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCt                                                                                 >  6:224762/1‑135 (MQ=255)
                                                     cGGGAAATGCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTCCCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTg                                                                                <  5:239307/135‑1 (MQ=255)
                                                       ggAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAg                                                                              >  5:337720/1‑135 (MQ=255)
                                                          aaTTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTgc                                                                           >  1:100161/1‑135 (MQ=255)
                                                          aaTTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTgc                                                                           >  2:234027/1‑135 (MQ=255)
                                                                  aGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGAt                                                                   <  1:83860/135‑1 (MQ=255)
                                                                        gcggcAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGCGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGa                                                             <  5:201198/135‑1 (MQ=255)
                                                                        gcggcAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGa                                                             <  6:115302/135‑1 (MQ=255)
                                                                                cGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAgg                                                      >  2:126200/1‑134 (MQ=255)
                                                                                cGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCt                                                    <  4:6565/136‑1 (MQ=255)
                                                                                 gCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCt                                                    <  7:216119/135‑1 (MQ=255)
                                                                                  ccTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTAt                                                  <  4:7975/136‑1 (MQ=255)
                                                                                       tCCTGCATGCCTACGCGTTGACTGTCCCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGc                                              <  7:319370/135‑1 (MQ=255)
                                                                                              tGCCTACGCGTTGACTGTACCGCTGCCGGACGGTGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGCTACTGAt                                       >  6:121582/1‑135 (MQ=255)
                                                                                               gCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTg                                                                                <  6:326016/93‑1 (MQ=255)
                                                                                               gCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTg                                                                                >  5:326016/1‑93 (MQ=255)
                                                                                                   aCGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGCTACTGATCTTCg                                  >  7:79892/1‑135 (MQ=255)
                                                                                                   aCGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGCTACTGATCTTCg                                  >  5:123244/1‑135 (MQ=255)
                                                                                                   aCGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGCTACTGATCTTCGa                                 >  4:317407/1‑136 (MQ=255)
                                                                                                   aCGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTATGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGCTACTGATCTTCGa                                 >  4:212354/1‑136 (MQ=255)
                                                                                                          gACTGTACCGCTGCCGGACGGCGGTGAATTGACGCTCGCGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGCTACTGATCTTCGACTGGGa                           <  4:198534/135‑1 (MQ=255)
                                                                                                           aCTGTACCGCTGCCGGACGGCGGTGAATTGCAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGCTACTGATCTTCGACTGGGa                           <  8:227427/134‑1 (MQ=255)
                                                                                                                aacGCTGCCGGACGTCGTTTAATTGAAGCTAGAGGCTAATGTTGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCAATTATCCATGAAAAAAGGCTATGAGCTACTGATCTTCGACTGGGACGGc                       <  5:161660/131‑1 (MQ=255)
                                                                                                                  cGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGCTACTGATCTTCGACTGGGACGGCAcc                    >  1:91684/1‑134 (MQ=255)
                                                                                                                            cggcggTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGCTACTGATCTTCGACTGGGACGGCACCCTAGCCGATTc         <  1:45631/135‑1 (MQ=255)
                                                                                                                                   gAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGCTACTGATCTTCGACTGGGACGGCACCCTAGCCGATTCCATTggg  <  5:121582/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                     
ATGCCGGGCACATGATCGCTGGCGACAGCAAATACGGCGACGAAGACTTCAGCCGGGAAATTCGCGAGTTGGGCGGCAAACGCCTGTTCCTGCATGCCTACGCGTTGACTGTACCGCTGCCGGACGGCGGTGAATTGAAGCTCGAGGCGCCTGTGGATGAGGTCTGGGCAAAAACCGTTGAACGTCTGAGTGCGCCCTGATCCATGAAAAAAGGCTATGAGCTACTGATCTTCGACTGGGACGGCACCCTAGCCGATTCCATTGGG  >  KT2440/2145131‑2145396

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