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

GAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGC  >  Exported/724417‑724554
                                                                      |                                                                    
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:95637/64‑1 (MQ=255)
        agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGaa                                                                    <  1:2075886/64‑1 (MQ=255)
        agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAga                                                                  >  1:736010/1‑66 (MQ=255)
        agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATc                                                             <  1:652954/71‑1 (MQ=255)
        agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATc                                                             <  1:2351795/71‑1 (MQ=255)
        agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATc                                                             <  1:2876918/71‑1 (MQ=255)
        agaagaGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGGCATCTT‑ATATGAAGATGATc                                                             <  1:1153155/71‑1 (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:2354177/1‑69 (MQ=255)
               gCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATc                                                        >  1:3578577/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:1761470/1‑67 (MQ=255)
                                 cAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAAc                                        <  1:223058/67‑1 (MQ=255)
                                 cAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTc                                    >  1:159988/1‑71 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTcc                                   >  1:2677790/1‑70 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                                  <  1:3425753/71‑1 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                                  <  1:3111083/71‑1 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                                  >  1:2039763/1‑71 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                                  <  1:1620963/71‑1 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                                  >  1:2381953/1‑71 (MQ=255)
                                   aaaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCg                                  >  1:704455/1‑71 (MQ=255)
                                    aaGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCgg                                 >  1:544815/1‑71 (MQ=255)
                                     aGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATca                                                            >  1:1730416/1‑43 (MQ=255)
                                        tGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTAcg                             >  1:1412349/1‑71 (MQ=255)
                                        tGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTAcg                             >  1:1448961/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:1100379/66‑1 (MQ=37)
                                                        gatCATCTTCATAT‑AAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg                <  1:506781/66‑1 (MQ=37)
                                                        gatCATCTTCATAT‑AAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg                <  1:1712763/66‑1 (MQ=37)
                                                        cGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg               >  1:1483304/1‑69 (MQ=255)
                                                         gTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAg            >  1:3051149/1‑71 (MQ=255)
                                                         gTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAg            <  1:3715022/71‑1 (MQ=255)
                                                         atCATCTTCATAT‑AAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg               >  1:2017459/2‑68 (MQ=38)
                                                         atCATCTTCATAT‑AAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcgg               >  1:864379/2‑68 (MQ=38)
                                                         atCATCTTCATAT‑AAGATGACCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATggcg                <  1:634243/66‑1 (MQ=37)
                                                          tCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATgg                  <  1:2645892/64‑1 (MQ=255)
                                                          tCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATgg                  <  1:1492587/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)
                                                                  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)
                                                                      |                                                                    
GAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGTCATCTT‑ATATGAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGC  >  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: ‑