Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,767,900 T→C 30.4% A266A (GCA→GCG rluD ← 23S rRNA pseudouridine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,767,9000TC30.4% 96.2 / 50.6 105A266A (GCA→GCGrluD23S rRNA pseudouridine synthase
Reads supporting (aligned to +/- strand):  ref base T (47/26);  new base C (26/6);  total (73/32)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

ATAAAGACGCAGCATGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATGA  >  minE/1767837‑1767966
                                                               |                                                                  
aTAAAGACGCAGCATGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCtt                                                                >  1:3181985/1‑68 (MQ=255)
aTAAAGACGCAGCATGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCtt                                                                >  1:2731356/1‑68 (MQ=255)
      aCGCAGCATGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGc                                                         <  1:385550/69‑1 (MQ=255)
              tGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTg                                                 >  1:3367383/1‑69 (MQ=255)
              tGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTg                                                 >  1:1875849/1‑69 (MQ=255)
              tGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTg                                                 >  1:205919/1‑69 (MQ=255)
              tGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTg                                                 >  1:2408783/1‑69 (MQ=255)
              tGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTg                                                 >  1:1585760/1‑69 (MQ=255)
              tGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTg                                                 >  1:326579/1‑69 (MQ=255)
              tGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTg                                                 >  1:1234761/1‑69 (MQ=255)
              tGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTg                                                 >  1:1169664/1‑69 (MQ=255)
                  tGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcgg                                             >  1:335392/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:2137089/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:1021003/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:2751178/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:2518702/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:3152469/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:2162216/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:2158391/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:210110/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:3142971/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:1857741/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:901622/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:1175286/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:1809948/1‑69 (MQ=255)
                        tAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:1200332/1‑69 (MQ=255)
                         aGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgga                                      >  1:2509351/1‑69 (MQ=255)
                         aGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgga                                      >  1:2922011/1‑69 (MQ=255)
                         aGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:626311/1‑68 (MQ=255)
                                      cAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGCCGCGGACGGCCACCATaaa                         >  1:1949561/1‑69 (MQ=255)
                                      cAAACATACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGgcggacg                                          >  1:2134087/1‑52 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCttt                                                   >  1:872583/1‑42 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcgg                                       >  1:875279/1‑54 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcg                                        >  1:2057640/1‑53 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcg                                        >  1:2612951/1‑53 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcg                                        >  1:2321250/1‑53 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcg                                        >  1:2205784/1‑53 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcg                                        >  1:1411284/1‑53 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcg                                        >  1:2165461/1‑53 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcg                                        >  1:1251030/1‑53 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcg                                        >  1:178447/1‑53 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcg                                        >  1:3159212/1‑53 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcg                                        >  1:518272/1‑53 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcg                                        >  1:2103319/1‑53 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgcg                                        >  1:1826473/1‑53 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGgcggacgc                                         >  1:2288670/1‑52 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAAc                        >  1:2724932/1‑69 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAAc                        >  1:1582851/1‑69 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAAc                        >  1:3136847/1‑69 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAAc                        >  1:2342165/1‑69 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAAc                        >  1:389419/1‑69 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAAc                        >  1:2179587/1‑69 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAAc                        >  1:586398/1‑69 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAAc                        >  1:2081271/1‑69 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAAc                        >  1:568832/1‑69 (MQ=255)
                                       aaaCTTACGCAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAAc                        >  1:58266/1‑69 (MQ=255)
                                          cTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCgg                     <  1:1571675/69‑1 (MQ=255)
                                          cTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCgg                     <  1:2484161/69‑1 (MQ=255)
                                          cTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCgg                     <  1:1492856/69‑1 (MQ=255)
                                          cTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCgg                     <  1:952107/69‑1 (MQ=255)
                                          cTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCgg                     >  1:417694/1‑69 (MQ=255)
                                               gCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCg                 <  1:350781/68‑1 (MQ=255)
                                               gCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCg                 <  1:368257/68‑1 (MQ=255)
                                               gCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCg                 <  1:2549569/68‑1 (MQ=255)
                                               gCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCg                 <  1:2049331/68‑1 (MQ=255)
                                               gCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCg                 <  1:2142536/68‑1 (MQ=255)
                                               gCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCg                 <  1:2498170/68‑1 (MQ=255)
                                               gCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCg                 <  1:2626500/68‑1 (MQ=255)
                                                cAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCg                 <  1:2029870/67‑1 (MQ=255)
                                                cAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGTGGACGCGGACGGCCACCATAAACCGGATCGcc               <  1:312718/69‑1 (MQ=255)
                                                cAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGTGGACGCGGACGGCCACCATAAACCGGATCGcc               <  1:896948/69‑1 (MQ=255)
                                                cAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGTGGACGCGGACGGCCACCATAAACCGGATCGcc               <  1:2931052/69‑1 (MQ=255)
                                                cAGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGTGGACGCAGACGGCCACCATAAACCGGATCGcc               <  1:2797685/69‑1 (MQ=255)
                                                 aGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGTGGACGCGGACGGCCACCATAAACCGGATCGcc               <  1:884003/68‑1 (MQ=255)
                                                 aGCGTGGAGATAAACGCTTCCGAAGCACCTTTTGGTGGACGCGGACGGCCACCATAAACCGGATCGcc               <  1:1775610/68‑1 (MQ=255)
                                                       gagaTAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCgg                     >  1:3150713/1‑56 (MQ=255)
                                                       gagaTAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCgg                     >  1:224301/1‑56 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGc                >  1:1596428/1‑56 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGc                >  1:928803/1‑56 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGc                >  1:1253057/1‑56 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGc                >  1:3111967/1‑56 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATg   <  1:1649592/69‑1 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATg   >  1:741887/1‑69 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATg   >  1:634153/1‑69 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATg   >  1:1658661/1‑69 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATg   <  1:95658/69‑1 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATg   >  1:1894020/1‑69 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATg   <  1:470617/69‑1 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATg   >  1:2566089/1‑69 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATg   >  1:2940314/1‑69 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACAATAAACCGGATCGCCCACCAGCGGATg   <  1:3151328/69‑1 (MQ=255)
                                                            aaaTGCTTCCGAAGCACCTTTTGGCGGACACGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATg   >  1:1452172/1‑69 (MQ=255)
                                                             aaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATg   <  1:2669353/68‑1 (MQ=255)
                                                             aaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATg   <  1:2980562/68‑1 (MQ=255)
                                                             aaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATGa  <  1:2263429/69‑1 (MQ=255)
                                                             aaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATGa  <  1:249417/69‑1 (MQ=255)
                                                             aaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATGa  <  1:3346905/69‑1 (MQ=255)
                                                             aaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATGa  <  1:2784975/69‑1 (MQ=255)
                                                             aaTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATGa  <  1:941972/69‑1 (MQ=255)
                                                              aTGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATGa  <  1:907053/68‑1 (MQ=255)
                                                               |                                                                  
ATAAAGACGCAGCATGGTTGCATGTAGCGCCTGGCGGTCAAACTTACGCAGCGTGGAGATAAATGCTTCCGAAGCACCTTTTGGCGGACGCGGACGGCCACCATAAACCGGATCGCCCACCAGCGGATGA  >  minE/1767837‑1767966

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