Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 724,481:1 +C 10.0% coding (287/960 nt) rluC → 23S rRNA pseudouridylate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE724,4811.C10.0% 204.0 / 25.7 60coding (287/960 nt)rluC23S rRNA pseudouridylate synthase
Reads supporting (aligned to +/- strand):  ref base . (31/23);  new base C (2/4);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CTGAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTT  >  minE/724415‑724549
                                                                   |                                                                    
cTGAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTt‑ata                                                                   <  1:2703634/70‑1 (MQ=255)
cTGAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑atat                                                                  >  1:1062477/1‑71 (MQ=255)
  gAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATg                                                                 >  1:3487787/1‑70 (MQ=255)
  gAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATg                                                                 >  1:1433102/1‑70 (MQ=255)
   aGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATg                                                                 >  1:1121806/1‑69 (MQ=255)
      gggAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAgat                                                            <  1:2564688/71‑1 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGaa                                                               <  1:2075886/64‑1 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGaa                                                               <  1:95637/64‑1 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAga                                                             >  1:736010/1‑66 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATc                                                        <  1:652954/71‑1 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATc                                                        <  1:2876918/71‑1 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATc                                                        <  1:2351795/71‑1 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGGCATCTT‑ATATGAAGATGATc                                                        <  1:1153155/71‑1 (MQ=255)
              gagGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTt‑ata                                                                   >  1:2688952/1‑56 (MQ=255)
              gagGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTt‑ata                                                                   >  1:2609385/1‑56 (MQ=255)
              gagGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTt‑ata                                                                   >  1:2281560/1‑56 (MQ=255)
               agGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATcac                                                      >  1:3728835/1‑68 (MQ=255)
               agGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATcac                                                      >  1:3675809/1‑68 (MQ=255)
                ggCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATcc                                                  >  1:510292/1‑71 (MQ=255)
                 gCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATc                                                   >  1:3578577/1‑69 (MQ=255)
                 gCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATc                                                   >  1:2354177/1‑69 (MQ=255)
                   ggTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTgg                                               >  1:2423243/1‑71 (MQ=255)
                        cGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCt                                           <  1:1946543/70‑1 (MQ=255)
                         gCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTg                                             <  1:3399475/67‑1 (MQ=255)
                         gCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTg                                          <  1:192704/70‑1 (MQ=255)
                          ccACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCAAATCCTGGTGCTGa                                         >  1:361702/1‑70 (MQ=255)
                           cacaTCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAAt                                       <  1:2924118/71‑1 (MQ=255)
                             caTCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATc                                                        <  1:1821595/52‑1 (MQ=255)
                               tCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTgg                                               <  1:302095/59‑1 (MQ=255)
                               tCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAAc                                   >  1:252207/1‑71 (MQ=255)
                                   cAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGt                                              >  1:274405/1‑56 (MQ=255)
                                   cAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAAc                                   <  1:223058/67‑1 (MQ=255)
                                   cAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAAc                                   >  1:1761470/1‑67 (MQ=255)
                                   cAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTc                               >  1:159988/1‑71 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTcc                              >  1:2677790/1‑70 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             >  1:2381953/1‑71 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             <  1:3425753/71‑1 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             <  1:3111083/71‑1 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             >  1:704455/1‑71 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             >  1:2039763/1‑71 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             <  1:1620963/71‑1 (MQ=255)
                                      aaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCgg                            >  1:544815/1‑71 (MQ=255)
                                       aGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATca                                                       >  1:1730416/1‑43 (MQ=255)
                                          tGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTAcg                        >  1:1448961/1‑71 (MQ=255)
                                          tGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTAcg                        >  1:1412349/1‑71 (MQ=255)
                                           ggcggcGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATa                                      <  1:2571812/56‑1 (MQ=255)
                                            gcggcgCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTAcggc                      >  1:2116782/1‑71 (MQ=255)
                                                   tGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACAt               >  1:476555/1‑71 (MQ=255)
                                                     gCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACAt               <  1:2002386/69‑1 (MQ=255)
                                                     gCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATgg             <  1:1449965/71‑1 (MQ=255)
                                                          gatCATCTTCATAT‑AAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg           <  1:506781/66‑1 (MQ=37)
                                                          gatCATCTTCATAT‑AAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg           <  1:1100379/66‑1 (MQ=37)
                                                          gatCATCTTCATAT‑AAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg           <  1:1712763/66‑1 (MQ=37)
                                                          cGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg          >  1:1483304/1‑69 (MQ=255)
                                                           gTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAg       <  1:3715022/71‑1 (MQ=255)
                                                           gTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAg       >  1:3051149/1‑71 (MQ=255)
                                                           atCATCTTCATAT‑AAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg          >  1:864379/2‑68 (MQ=38)
                                                           atCATCTTCATAT‑AAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg          >  1:2017459/2‑68 (MQ=38)
                                                           atCATCTTCATAT‑AAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg           <  1:634243/66‑1 (MQ=37)
                                                            tCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATgg             <  1:1492587/64‑1 (MQ=255)
                                                            tCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATgg             <  1:2645892/64‑1 (MQ=255)
                                                             cATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             >  1:2077664/1‑47 (MQ=255)
                                                             cATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACa                >  1:404687/1‑60 (MQ=255)
                                                               tCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg          >  1:1870086/1‑64 (MQ=255)
                                                               tCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGt   >  1:1575678/1‑71 (MQ=255)
                                                                 tt‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGtt  >  1:1208846/1‑70 (MQ=255)
                                                                 tt‑ATAAGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTAc                 >  1:1608227/1‑55 (MQ=255)
                                                                   |                                                                    
CTGAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTT  >  minE/724415‑724549

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