Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,547,431 +G 11.8% coding (1498/1542 nt) ppx → exopolyphosphatase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,547,4311.→G11.8% 14.2 ‑2.1 ‑0.534coding (1498/1542 nt)ppxexopolyphosphatase
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (16/14):  total (20/14)
Fisher's exact test strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.54e-01

AGGAGCAAGAATACTGGGAAGGCGTGGCTGGCTGGC.GGTTGAAAATTGAAGAAGAAAGTACACCAGAAATC  >  REL606/2547396‑2547466
                                    |                                   
aGGAGCAAGAATACTGGGAAGTCGtggctggctggc                                      <  1:3583202/36‑1 (MQ=255)
aGGAGCAAGAATACTGGGAAGGCGtggctggctggc                                      <  1:6579497/36‑1 (MQ=255)
 ggAGCAAGAATACTGGGAAGGCGTGGCTGGCTGGC.g                                    >  1:1788298/1‑36 (MQ=255)
 ggAGCAAGAATACTGGGAAGGCGTGGCTGGCTGGC.g                                    >  1:1956094/1‑36 (MQ=255)
  gAGCAAGAATACTGGGAAGGCGTGGCTGGCTGGC.gg                                   <  1:970691/36‑1 (MQ=255)
   aGCAAGAATACTGGGAAGGCGTGGCTGGCTGGC.GGt                                  <  1:2350643/36‑1 (MQ=255)
   aGCAAGAATACTGGGAAGGCGTGGCTGGCTGGC.GGt                                  <  1:6213791/36‑1 (MQ=255)
    gCAAGAATACTGGGAAGGCGTGGCTGGCTGGC.GGtt                                 <  1:6553234/36‑1 (MQ=255)
      aaGAATACTGGGAAGGCGTGGCTGGCTGGC.GGTTgg                               >  1:3199794/1‑35 (MQ=255)
        gAATACTGGGAAGGCGTGGCTGGCTGGCGGGTTGaa                              >  1:4698548/1‑36 (MQ=37)
        gAATACTGGGAAGGCGTGGCTGGCTGGC.GGTTGtaa                             >  1:7925967/1‑36 (MQ=255)
         aaTACTGGGAAGGCGTGGCTGGCTGGCGGGTTGaaa                             >  1:986544/1‑36 (MQ=37)
         aaTACTGGGAAGGCGTGGCTGGCTGGC.GGTTGaaaa                            <  1:5242699/36‑1 (MQ=255)
         aaTACTGGGAAGGCGTGGCTGGCTGGC.GGTTGaaaa                            <  1:66838/36‑1 (MQ=255)
           tACTGGGAAGGCGTGGCTGGCTGGC.GGTTGGAAAtt                          >  1:6685121/1‑36 (MQ=255)
           tACTGGGAAGGCGTGGCTGGCTGGC.GGTTGAAAAtt                          <  1:3203957/36‑1 (MQ=255)
           gaCTGGGAAGGCGTGGCTGGCTGGC.GGTTGAAAAtt                          <  1:1421905/35‑1 (MQ=255)
            aCTGGGAAGGCGTGGCTGGCTGGG.GGTTTGAAAtta                         >  1:1863295/1‑35 (MQ=21)
            aCTGGGAAGGCGTGGCTGGCTGGC.GGTTGGAAAttt                         >  1:8417857/1‑35 (MQ=255)
            aCTGGGAAGGCGTGGCTGGCTGGC.GGGTTGAAATTg                         >  1:5794551/1‑36 (MQ=25)
             cTGGGAAGGCGTGGCTGGCTGGC.GGTTGAAAATTga                        <  1:7625710/36‑1 (MQ=255)
              tGGGAAGGCGTGGCTGGCTGGC.GGTTGGAAAtttga                       >  1:7795595/1‑33 (MQ=37)
              tGGGAAGGCGTGGCTGGCTGGC.GGTTGAAAATTtaa                       >  1:1045203/1‑36 (MQ=255)
              tGGGAAGGCGTGGCTGGCTGGC.GGTTGAAAATTgaa                       <  1:7534914/36‑1 (MQ=255)
              tGGGAAGGCGTGGCTGGCTGGC.GGTTGAAAATTgaa                       <  1:2404195/36‑1 (MQ=255)
               gggAAGGCGTGGCTGGGTGGCGGGTTGAAATttttg                       >  1:7089097/1‑33 (MQ=16)
               gggAAGGCGTGGCTGGCTGGG.GGGTGAAAATTtaag                      >  1:1101001/1‑36 (MQ=25)
               gggAAGGCGTGGCTGGCTGGC.GGTTGAAAATTgaag                      <  1:78981/36‑1 (MQ=255)
               gggAAGGCGTGGCTGGCTGGC.GGTTGAAAATTgaag                      <  1:818855/36‑1 (MQ=255)
               gggAAGGCGTGGCTGGCTGGC.GGTTGAAAAATtaag                      >  1:7857592/1‑36 (MQ=37)
                ggAAGGCGTGGCTGGCTGGCGGGGTAAAAATTgaga                      >  1:6747164/1‑34 (MQ=16)
                 gAAGGCGTGGCTGGCTGGC.GGTTGAAAATTgaagaa                    <  1:7124474/36‑1 (MQ=255)
                 gAAGGCGTGGCTGGCTGGC.GGGTGAAAATTtaagga                    >  1:7379928/1‑34 (MQ=25)
                  aaGGCGTGGCTGGCTGGC.GGGTGGAAAATaaagaag                   >  1:4926780/1‑36 (MQ=16)
                   aGGCGTGGCTGGCTGGC.GGGTGGAAATTttagaagg                  >  1:3041009/1‑35 (MQ=16)
                    ggCGTGGCTGGCTGGC.GGTTGGAAAATttagaagga                 >  1:5818154/1‑34 (MQ=16)
                    ggCGTGGCTGGCTGGC.GGTTGAAAATTgaagaagaa                 <  1:3528998/36‑1 (MQ=255)
                     gCGTGGCTGGCTGGC.GGTTGAAAATTTGAGAAGaaa                >  1:3433762/1‑36 (MQ=25)
                       gTGGCTGGCTGGC.GGTTTGAAATTTAAGAAGGAAga              >  1:782604/1‑35 (MQ=16)
                        tggctggctggc.GGTTGAAAATTGAAGAAGAAAGTa             <  1:8298294/36‑1 (MQ=255)
                          gctggctggc.GGTTTAAAATTGAAGAAGAAAGTaca           >  1:7849211/1‑36 (MQ=255)
                          gctggctggc.GGTTGAAAATTGAAGAAGAAAGTaca           <  1:7452477/36‑1 (MQ=255)
                              gctggc.GGTTGAAAATTGAAGAAGAAAGTACACCAg       <  1:6265226/36‑1 (MQ=255)
                               ctggc.GGTTGAAAATTGAAGAAGAAAGTACACCAGa      <  1:3674541/36‑1 (MQ=255)
                                 ggc.GGTTGAAAATTGAAGAAGAAAGTACACCAGaaa    <  1:7220415/36‑1 (MQ=255)
                                   c.GGTTGAAAATTGAAGAAGAAAGTACACCAGAAATc  <  1:7358494/36‑1 (MQ=255)
                                   c.GGTTGAAAATTGAAGAAGAAAGTACACCAGAAATc  <  1:4219968/36‑1 (MQ=255)
                                    |                                   
AGGAGCAAGAATACTGGGAAGGCGTGGCTGGCTGGC.GGTTGAAAATTGAAGAAGAAAGTACACCAGAAATC  >  REL606/2547396‑2547466

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 17 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG