Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,434,339 T→G 15.7% intergenic (+91/‑45) 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,3390TG15.7% 173.0 / 25.2 76intergenic (+91/‑45)rhlB/gppATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base G (4/8);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGC  >  Exported/2434270‑2434409
                                                                     |                                                                      
aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGt                                                                        >  1:192398/1‑70 (MQ=255)
aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGt                                                                        >  1:5154368/1‑70 (MQ=255)
             tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAg                                                            <  1:3800996/69‑1 (MQ=255)
             tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAg                                                            <  1:1632874/69‑1 (MQ=255)
             tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGc                                                           >  1:1956882/1‑70 (MQ=255)
               tCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCaa                                                        >  1:398427/1‑71 (MQ=255)
                 aaGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAaca                                                      >  1:2247963/1‑71 (MQ=255)
                 aaGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAac                                                       <  1:80674/70‑1 (MQ=255)
                  aGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAaca                                                      >  1:4882327/1‑70 (MQ=255)
                   gATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGaa                                                                 >  1:3420785/1‑58 (MQ=255)
                     tGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAAcaca                                                    >  1:4886431/1‑69 (MQ=255)
                     cggCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCaga                                                  <  1:3762691/70‑1 (MQ=2)
                       gCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAgaga                                                <  1:769099/71‑1 (MQ=255)
                          agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg                                             <  1:5298159/71‑1 (MQ=255)
                          agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg                                             <  1:2893665/71‑1 (MQ=255)
                          agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg                                             <  1:4304180/71‑1 (MQ=255)
                           gagaATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCaa                                                        >  1:279177/1‑59 (MQ=255)
                             gaatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc                                                          >  1:4161099/1‑55 (MQ=255)
                              aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATg                                                               <  1:4740567/49‑1 (MQ=255)
                              aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATg                                                               <  1:2961155/49‑1 (MQ=255)
                              aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGAtt                                          >  1:4042398/1‑70 (MQ=255)
                              aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTg                                         <  1:5089868/71‑1 (MQ=255)
                               atgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg                                             <  1:1837188/66‑1 (MQ=255)
                               atgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTg                                         <  1:1370517/70‑1 (MQ=255)
                                tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGaa                                                                 >  1:1316831/1‑45 (MQ=255)
                                tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc                                                          >  1:724454/1‑52 (MQ=255)
                                tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc                                                          >  1:3701032/1‑52 (MQ=255)
                                tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc                                                          >  1:4582472/1‑52 (MQ=255)
                                tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTg                                         <  1:955518/69‑1 (MQ=255)
                                 gaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAacac                                                     <  1:4405530/56‑1 (MQ=255)
                                  aatCCCGATGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAg                                      >  1:2506493/1‑70 (MQ=9)
                                  aatCCCGATGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAg                                      >  1:1530604/1‑70 (MQ=9)
                                   atCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACag                                                   >  1:3718940/1‑56 (MQ=255)
                                   atCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACag                                                   >  1:5572557/1‑56 (MQ=255)
                                   atCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACag                                                   >  1:672959/1‑56 (MQ=255)
                                     cccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGa                                    >  1:1279910/1‑69 (MQ=255)
                                      ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                                  <  1:3958831/70‑1 (MQ=255)
                                      ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                                  >  1:3322455/1‑70 (MQ=255)
                                      ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                                  >  1:594413/1‑70 (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:1367101/67‑1 (MQ=9)
                                         aTGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATg                                  <  1:128231/67‑1 (MQ=9)
                                         aTGAGCTTACTCAGGTAAGTGATTCGGGGGAGCGAACGCAGCCAAGGCAGAGGCGGCTTGAAGGATg                                  <  1:1081083/67‑1 (MQ=9)
                                          tGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAagag                             <  1:3058619/71‑1 (MQ=255)
                                          tGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAagag                             <  1:564370/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)
                                              cTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACaga                                                  <  1:3382326/46‑1 (MQ=255)
                                              cTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACaga                                                  <  1:567898/46‑1 (MQ=255)
                                               ttACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACaga                                                  <  1:1212478/45‑1 (MQ=255)
                                               ttACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACag                                                   <  1:5697365/44‑1 (MQ=255)
                                                tACTGTAGTAAGTGATTCAGGTGAAAGAATGAAGCCaa                                                        >  1:2602198/1‑38 (MQ=255)
                                                 aCTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCa                                                         <  1:4020884/36‑1 (MQ=255)
                                                 aCTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAgaga                                                >  1:4518654/1‑45 (MQ=255)
                                                 aCTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAgaga                                                >  1:3698710/1‑45 (MQ=255)
                                                 aCTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTa                           >  1:5067942/1‑66 (MQ=255)
                                                    gtagtaAGTGATTCACGTGAAAGAATGCAGCCAACACag                                                   >  1:3747104/1‑39 (MQ=255)
                                                     tagtaAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCa                  >  1:2404347/1‑71 (MQ=255)
                                                        taAGTGATTCAGGTGAAAGAATGCAGCCAACACaga                                                  <  1:618693/36‑1 (MQ=255)
                                                          aGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGa                                    >  1:1790166/1‑48 (MQ=255)
                                                            tGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCtcgtc           >  1:1550582/1‑71 (MQ=255)
                                                            tGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCtcgtc           >  1:2311966/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:3243562/47‑1 (MQ=255)
                                                               ttCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGaa                                <  1:1675654/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)
                                                                     |                                                                      
ATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCACCTCGTCGCTGTATGC  >  Exported/2434270‑2434409

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: