Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,434,367 A→C 13.4% intergenic (+119/‑17) rhlB → / → gpp ATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,434,3670AC13.4% 154.9 / 19.7 67intergenic (+119/‑17)rhlB/gppATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base A (34/24);  new base C (4/5);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

AATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCATTGATCTCGGTTCGAATAGTTT  >  Exported/2434300‑2434436
                                                                   |                                                                     
aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGAtt                                                                     >  1:4042398/1‑70 (MQ=255)
aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTg                                                                    <  1:5089868/71‑1 (MQ=255)
 atgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTg                                                                    <  1:1370517/70‑1 (MQ=255)
  tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTg                                                                    <  1:955518/69‑1 (MQ=255)
    aatCCCGATGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAg                                                                 >  1:2506493/1‑70 (MQ=9)
    aatCCCGATGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAg                                                                 >  1:1530604/1‑70 (MQ=9)
       cccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGa                                                               >  1:1279910/1‑69 (MQ=255)
        ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                                                             >  1:3322455/1‑70 (MQ=255)
        ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                                                             >  1:594413/1‑70 (MQ=255)
        ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                                                             <  1:3958831/70‑1 (MQ=255)
        ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                                                             <  1:121317/70‑1 (MQ=255)
        ccGATGAGCTTACTGTAGTAAGTGATTCAGGGGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                                                             <  1:2059265/70‑1 (MQ=255)
        ccGATGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATg                                                             >  1:3565310/1‑70 (MQ=9)
        ccGATGAGCTCACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATg                                                             <  1:4671174/70‑1 (MQ=255)
         cGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa                                                           <  1:1773347/71‑1 (MQ=255)
           aTGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa                                                           >  1:5325034/1‑69 (MQ=255)
           aTGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATg                                                             <  1:1081083/67‑1 (MQ=9)
           aTGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATg                                                             <  1:128231/67‑1 (MQ=9)
           aTGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATg                                                             <  1:1367101/67‑1 (MQ=9)
            tGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAagag                                                        <  1:564370/71‑1 (MQ=255)
            tGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAagag                                                        <  1:3058619/71‑1 (MQ=255)
            tGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAagag                                                        <  1:5655139/71‑1 (MQ=255)
             gAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGCGACAGATTGAAGGATGAagag                                                        >  1:2168033/1‑70 (MQ=255)
             gAGCTTACTCAGGTAAGTGATTCGTGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATGAagtg                                                        >  1:3677193/1‑68 (MQ=255)
              aGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTa                                                      >  1:3899367/1‑71 (MQ=255)
              aGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATGAAgtgt                                                       <  1:2636176/70‑1 (MQ=9)
                   aCTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTa                                                      >  1:5067942/1‑66 (MQ=255)
                       tagtaAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCa                                             >  1:2404347/1‑71 (MQ=255)
                            aGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGa                                                               >  1:1790166/1‑48 (MQ=255)
                              tGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCtcgtc                                      >  1:2311966/1‑71 (MQ=255)
                              tGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCtcgtc                                      >  1:1550582/1‑71 (MQ=255)
                              tGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCtcgtc                                      >  1:2891341/1‑71 (MQ=255)
                              tGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATCAAGAGTATGGGTTCCACAtcgtc                                      >  1:5458551/1‑71 (MQ=255)
                               gATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCtcgtcg                                     >  1:4667734/1‑71 (MQ=255)
                                 ttCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa                                                           <  1:1675654/47‑1 (MQ=255)
                                 ttCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa                                                           <  1:3243562/47‑1 (MQ=255)
                                 ttCAGGGGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa                                                           <  1:3009514/47‑1 (MQ=255)
                                   cAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTg                                  <  1:5350950/70‑1 (MQ=255)
                                    aGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGt                                 <  1:2491598/70‑1 (MQ=255)
                                       tGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa                                                           >  1:3101460/1‑41 (MQ=255)
                                       tGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGc                             <  1:3294857/71‑1 (MQ=255)
                                        gAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                                                  >  1:3942289/1‑49 (MQ=255)
                                        gAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                                                  >  1:1985509/1‑49 (MQ=255)
                                         aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                                                  >  1:962305/1‑48 (MQ=255)
                                         aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                                                  >  1:1439269/1‑48 (MQ=255)
                                         aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                                                  >  1:922984/1‑48 (MQ=255)
                                         aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                                                  >  1:5610772/1‑48 (MQ=255)
                                         aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATggg                                                  >  1:4383222/1‑48 (MQ=255)
                                         aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATgg                                                   >  1:3802152/1‑47 (MQ=255)
                                         aaaGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATgg                                                   >  1:4693506/1‑47 (MQ=255)
                                            gAATGCAGCCAACACAGAGAGAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCa                        >  1:4839376/1‑71 (MQ=255)
                                            gAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCa                        >  1:1017520/1‑71 (MQ=255)
                                            gAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCa                        >  1:1944370/1‑71 (MQ=255)
                                             aaTGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCt                                          >  1:2340350/1‑52 (MQ=255)
                                             aaTGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCt                                          >  1:2115110/1‑52 (MQ=255)
                                             aaTGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCt                                          >  1:2739377/1‑52 (MQ=255)
                                               tGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCATTg                     <  1:2741606/71‑1 (MQ=255)
                                                gCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCATTg                     <  1:4640623/70‑1 (MQ=255)
                                                            gagaCAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCt                                   >  1:4109164/1‑44 (MQ=255)
                                                            gagaCAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCATTGATCTCGGTTCGaa        <  1:2745846/71‑1 (MQ=255)
                                                            gagaCAGATTGAAGGATGAAGAGTATGGGTTCCAACTCGTCGc                                    >  1:2279848/1‑43 (MQ=255)
                                                               aCAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCATTGATCTCGGTTCGAATAg     >  1:1756418/1‑71 (MQ=255)
                                                                  gATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCAt                       <  1:276422/50‑1 (MQ=255)
                                                                  gATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCAt                       <  1:2676765/50‑1 (MQ=255)
                                                                  gATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCAt                       <  1:5030767/50‑1 (MQ=255)
                                                                  gATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCAt                       <  1:3508596/50‑1 (MQ=255)
                                                                  gATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCATTGATCTCGGTTCGAATAGttt  <  1:1681249/71‑1 (MQ=255)
                                                                   aTTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCATTGATCTCGGTTCGAATAGttt  <  1:2790234/70‑1 (MQ=255)
                                                                   |                                                                     
AATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGCAGCCATTGATCTCGGTTCGAATAGTTT  >  Exported/2434300‑2434436

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: