Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 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
*Exported724,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

CTGAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTT  >  Exported/724415‑724549
                                                                   |                                                                    
cTGAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTtata                                                                   <  1:2703634/70‑1 (MQ=255)
cTGAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTatat                                                                  >  1:1062477/1‑71 (MQ=255)
  gAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATg                                                                 >  1:3487787/1‑70 (MQ=255)
  gAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATg                                                                 >  1:1433102/1‑70 (MQ=255)
   aGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATg                                                                 >  1:1121806/1‑69 (MQ=255)
      gggAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAgat                                                            <  1:2564688/71‑1 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGaa                                                               <  1:2075886/64‑1 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGaa                                                               <  1:95637/64‑1 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAga                                                             >  1:736010/1‑66 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATc                                                        <  1:652954/71‑1 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATc                                                        <  1:2876918/71‑1 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATc                                                        <  1:2351795/71‑1 (MQ=255)
          agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGGCATCTTATATGAAGATGATc                                                        <  1:1153155/71‑1 (MQ=255)
              gagGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTtata                                                                   >  1:2688952/1‑56 (MQ=255)
              gagGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTtata                                                                   >  1:2609385/1‑56 (MQ=255)
              gagGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTtata                                                                   >  1:2281560/1‑56 (MQ=255)
               agGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATcac                                                      >  1:3728835/1‑68 (MQ=255)
               agGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATcac                                                      >  1:3675809/1‑68 (MQ=255)
                ggCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATcc                                                  >  1:510292/1‑71 (MQ=255)
                 gCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATc                                                   >  1:3578577/1‑69 (MQ=255)
                 gCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATc                                                   >  1:2354177/1‑69 (MQ=255)
                   ggTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTgg                                               >  1:2423243/1‑71 (MQ=255)
                        cGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCt                                           <  1:1946543/70‑1 (MQ=255)
                         gCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTg                                             <  1:3399475/67‑1 (MQ=255)
                         gCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTg                                          <  1:192704/70‑1 (MQ=255)
                          ccACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCAAATCCTGGTGCTGa                                         >  1:361702/1‑70 (MQ=255)
                           cacaTCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAAt                                       <  1:2924118/71‑1 (MQ=255)
                             caTCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATc                                                        <  1:1821595/52‑1 (MQ=255)
                               tCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTgg                                               <  1:302095/59‑1 (MQ=255)
                               tCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAAc                                   >  1:252207/1‑71 (MQ=255)
                                   cAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGt                                              >  1:274405/1‑56 (MQ=255)
                                   cAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAAc                                   <  1:223058/67‑1 (MQ=255)
                                   cAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAAc                                   >  1:1761470/1‑67 (MQ=255)
                                   cAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTc                               >  1:159988/1‑71 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTcc                              >  1:2677790/1‑70 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             >  1:2381953/1‑71 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             <  1:3425753/71‑1 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             <  1:3111083/71‑1 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             >  1:704455/1‑71 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             >  1:2039763/1‑71 (MQ=255)
                                     aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             <  1:1620963/71‑1 (MQ=255)
                                      aaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCgg                            >  1:544815/1‑71 (MQ=255)
                                       aGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATca                                                       >  1:1730416/1‑43 (MQ=255)
                                          tGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTAcg                        >  1:1448961/1‑71 (MQ=255)
                                          tGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTAcg                        >  1:1412349/1‑71 (MQ=255)
                                           ggcggcGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATa                                      <  1:2571812/56‑1 (MQ=255)
                                            gcggcgCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTAcggc                      >  1:2116782/1‑71 (MQ=255)
                                                   tGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACAt               >  1:476555/1‑71 (MQ=255)
                                                     gCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACAt               <  1:2002386/69‑1 (MQ=255)
                                                     gCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATgg             <  1:1449965/71‑1 (MQ=255)
                                                          gatCATCTTCATATAAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg           <  1:506781/66‑1 (MQ=37)
                                                          gatCATCTTCATATAAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg           <  1:1100379/66‑1 (MQ=37)
                                                          gatCATCTTCATATAAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg           <  1:1712763/66‑1 (MQ=37)
                                                          cGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg          >  1:1483304/1‑69 (MQ=255)
                                                           gTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAg       <  1:3715022/71‑1 (MQ=255)
                                                           gTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAg       >  1:3051149/1‑71 (MQ=255)
                                                           atCATCTTCATATAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg          >  1:864379/2‑68 (MQ=38)
                                                           atCATCTTCATATAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg          >  1:2017459/2‑68 (MQ=38)
                                                           atCATCTTCATATAAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg           <  1:634243/66‑1 (MQ=37)
                                                            tCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATgg             <  1:1492587/64‑1 (MQ=255)
                                                            tCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATgg             <  1:2645892/64‑1 (MQ=255)
                                                             cATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                             >  1:2077664/1‑47 (MQ=255)
                                                             cATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACa                >  1:404687/1‑60 (MQ=255)
                                                               tCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg          >  1:1870086/1‑64 (MQ=255)
                                                               tCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGt   >  1:1575678/1‑71 (MQ=255)
                                                                 ttATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGtt  >  1:1208846/1‑70 (MQ=255)
                                                                 ttATAAGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTAc                 >  1:1608227/1‑55 (MQ=255)
                                                                   |                                                                    
CTGAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTT  >  Exported/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: