Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,434,323 T→A 15.8% intergenic (+75/‑61) 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,3230TA15.8% 209.6 / 29.2 82intergenic (+75/‑61)rhlB/gppATP‑dependent RNA helicase/guanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base T (38/31);  new base A (6/7);  total (44/38)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCA  >  Exported/2434256‑2434393
                                                                   |                                                                      
aTGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtag                                                                      >  1:2023369/1‑70 (MQ=255)
   tCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtagt                                                                     <  1:1230803/68‑1 (MQ=255)
    taTACACAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGt                                                                 >  1:4293028/2‑71 (MQ=9)
     aTACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGt                                                                 >  1:5601956/1‑70 (MQ=255)
      tACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGa                                                               >  1:3953088/1‑71 (MQ=255)
      tACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGa                                                               >  1:195662/1‑71 (MQ=255)
      tACACAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGTGa                                                               <  1:4024673/71‑1 (MQ=9)
         acAAAATCATTCAAGTTGCATCAAGGCGGCAAGTGAGCGAATCCCGATGAGCTTACTCAGGTAAGTGATTc                                                            >  1:2925141/2‑71 (MQ=12)
           aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa                                                           <  1:4193246/70‑1 (MQ=255)
           aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCa                                                           >  1:3893752/1‑70 (MQ=255)
           aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAg                                                          >  1:5458564/1‑71 (MQ=255)
           aaaaTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAg                                                          >  1:2517386/1‑71 (MQ=255)
              aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGAtt                                                             >  1:4116376/1‑65 (MQ=255)
              aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGt                                                        >  1:192398/1‑70 (MQ=255)
              aTCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGt                                                        >  1:5154368/1‑70 (MQ=255)
                 aTTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTg                                                                >  1:5533870/1‑59 (MQ=255)
                 aTTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTg                                                                >  1:2951237/1‑59 (MQ=255)
                 aTTCAAACTGTATCAAGATGGCAAGAGAATGAATCACGATGAGCTTACTGTAGTAAGTg                                                                >  1:3781978/1‑59 (MQ=255)
                    cAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtag                                                                      >  1:5087141/1‑50 (MQ=255)
                    cAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTgtag                                                                      >  1:1442162/1‑50 (MQ=255)
                      aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                                                            <  1:5226499/58‑1 (MQ=255)
                      aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                                                            <  1:2491927/58‑1 (MQ=255)
                      aaCTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTc                                                            <  1:4275211/58‑1 (MQ=255)
                           tATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGAtt                                                             <  1:55517/52‑1 (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:4304180/71‑1 (MQ=255)
                                        agagAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAg                             <  1:2893665/71‑1 (MQ=255)
                                         gagaATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCaa                                        >  1:279177/1‑59 (MQ=255)
                                           gaatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc                                          >  1:4161099/1‑55 (MQ=255)
                                            aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATg                                               <  1:2961155/49‑1 (MQ=255)
                                            aatgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATg                                               <  1:4740567/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:3701032/1‑52 (MQ=255)
                                              tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc                                          >  1:4582472/1‑52 (MQ=255)
                                              tgaatCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGcc                                          >  1:724454/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:672959/1‑56 (MQ=255)
                                                 atCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACag                                   >  1:5572557/1‑56 (MQ=255)
                                                   cccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGa                    >  1:1279910/1‑69 (MQ=255)
                                                    ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                  <  1:3958831/70‑1 (MQ=255)
                                                    ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                  <  1:121317/70‑1 (MQ=255)
                                                    ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                  >  1:594413/1‑70 (MQ=255)
                                                    ccGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATg                  >  1:3322455/1‑70 (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:3698710/1‑45 (MQ=255)
                                                               aCTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAgaga                                >  1:4518654/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)
                                                                   |                                                                      
ATGTCATACTCAAAATCATTCAAACTGTATCAAGATGGCAAGAGAATGAATCCCGATGAGCTTACTGTAGTAAGTGATTCAGGTGAAAGAATGCAGCCAACACAGAGACAGATTGAAGGATGAAGAGTATGGGTTCCA  >  Exported/2434256‑2434393

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: