Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 724,486 Δ1 bp 8.8% coding (292/960 nt) rluC → 23S rRNA pseudouridylate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported724,4860G. 8.8% 223.2 / 12.2 68coding (292/960 nt)rluC23S rRNA pseudouridylate synthase
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base . (2/4);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGC  >  Exported/724417‑724554
                                                                      |                                                                    
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:95637/64‑1 (MQ=255)
        agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGaa                                                                    <  1:2075886/64‑1 (MQ=255)
        agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAga                                                                  >  1:736010/1‑66 (MQ=255)
        agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATc                                                             <  1:652954/71‑1 (MQ=255)
        agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATc                                                             <  1:2351795/71‑1 (MQ=255)
        agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATc                                                             <  1:2876918/71‑1 (MQ=255)
        agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGGCATCTTATATGAAGATGATc                                                             <  1:1153155/71‑1 (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:2354177/1‑69 (MQ=255)
               gCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATc                                                        >  1:3578577/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:1761470/1‑67 (MQ=255)
                                 cAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAAc                                        <  1:223058/67‑1 (MQ=255)
                                 cAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTc                                    >  1:159988/1‑71 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTcc                                   >  1:2677790/1‑70 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                                  <  1:3425753/71‑1 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                                  <  1:3111083/71‑1 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                                  >  1:2039763/1‑71 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                                  <  1:1620963/71‑1 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                                  >  1:2381953/1‑71 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                                  >  1:704455/1‑71 (MQ=255)
                                    aaGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCgg                                 >  1:544815/1‑71 (MQ=255)
                                     aGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATca                                                            >  1:1730416/1‑43 (MQ=255)
                                        tGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTAcg                             >  1:1412349/1‑71 (MQ=255)
                                        tGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTAcg                             >  1:1448961/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:1100379/66‑1 (MQ=37)
                                                        gatCATCTTCATATAAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg                <  1:506781/66‑1 (MQ=37)
                                                        gatCATCTTCATATAAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg                <  1:1712763/66‑1 (MQ=37)
                                                        cGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg               >  1:1483304/1‑69 (MQ=255)
                                                         gTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAg            >  1:3051149/1‑71 (MQ=255)
                                                         gTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAg            <  1:3715022/71‑1 (MQ=255)
                                                         atCATCTTCATATAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg               >  1:2017459/2‑68 (MQ=38)
                                                         atCATCTTCATATAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg               >  1:864379/2‑68 (MQ=38)
                                                         atCATCTTCATATAAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg                <  1:634243/66‑1 (MQ=37)
                                                          tCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATgg                  <  1:2645892/64‑1 (MQ=255)
                                                          tCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATgg                  <  1:1492587/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)
                                                                  atatGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACAt                    >  1:1579559/1‑55 (MQ=255)
                                                                  atatGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACAt                    >  1:3704702/1‑55 (MQ=255)
                                                                  atatGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGttt      <  1:1421489/69‑1 (MQ=255)
                                                                  atatGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGtt       <  1:1154109/68‑1 (MQ=255)
                                                                  atatGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGt        <  1:213656/67‑1 (MQ=255)
                                                                   tatGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGtt       <  1:2114344/67‑1 (MQ=255)
                                                                   tatGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGtt       <  1:1874626/67‑1 (MQ=255)
                                                                    atGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGc  >  1:3681373/1‑71 (MQ=255)
                                                                     tGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGc  <  1:3383413/70‑1 (MQ=255)
                                                                      |                                                                    
GAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTTATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGC  >  Exported/724417‑724554

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: