Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,947,7970A→G 5.4% 0.9 0.0 ‑0.25510aslAacrylsulfatase‑like enzyme
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (2/1):  ref base (35/17):  total (37/18)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.01e-01

GGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTC  >  REL606/3947762‑3947832
                                   |                                   
ggcgACAGCTACAACAAGAAGTTTGGGTGAAAACGa                                     >  1:3415449/1‑36 (MQ=255)
ggcgACAGCTACAACAAGAAGTTTGGGTGAAAACGa                                     >  1:2117166/1‑36 (MQ=255)
ggcgACAGCTACAACAAGAAGTTTGGGTGAAAACGa                                     <  1:5362845/36‑1 (MQ=255)
 gcgACAGCTACAACAAGAAGTTTGGGTGAAAACGaa                                    >  1:4578155/1‑36 (MQ=255)
 gcgACAGCTACAACAAGAAGTTTGGGTGAAAACGaa                                    >  1:5164657/1‑36 (MQ=255)
 gcgACAGCTACAACAAGAAGTTTGGGTGAAAACGaa                                    <  1:1598173/36‑1 (MQ=255)
  cgACAGCTACAACAAGAAGTTTGGGTGAAAACGaaa                                   <  1:5069976/36‑1 (MQ=255)
   gACAGCTACAACAAGAAGTTTGGGTGAAAACGaaaa                                  >  1:384351/1‑36 (MQ=255)
   gACAGCTACAACAAGAAGTTTGGGTGAAAACGaaaa                                  <  1:2474653/36‑1 (MQ=255)
   gACAGCTACAACAAGAAGTTTGGGTGAAAACGaaaa                                  <  1:5725955/36‑1 (MQ=255)
    aCAGCTACAACAAGAAGTTTGGGTGAAAACGAAAAt                                 <  1:4351662/36‑1 (MQ=255)
     cAGCTACAACAAGAAGTTTTGTTGAAAACGAAAAtt                                <  1:2693170/36‑1 (MQ=37)
     cAGCTACAACAAGAAGTTTGGGTGAAAACGAAAAtt                                <  1:3124154/36‑1 (MQ=255)
      aGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTc                               >  1:3652649/1‑36 (MQ=255)
      aGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTc                               >  1:3750030/1‑36 (MQ=255)
      aGCTACAACAAGAAGTTTGGGCGAAAACGAAAATTc                               >  1:3037292/1‑36 (MQ=255)
       gCTACAACAAGAAGTTTGGGTGAAAACGGAAATTcc                              <  1:1860854/36‑1 (MQ=255)
       gCTACAACAAGAAGTTTGGGTGAAAACGAAAATTcc                              >  1:5929357/1‑36 (MQ=255)
        cTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCa                             <  1:1657905/36‑1 (MQ=255)
         tACAACAAGAAGTTTGGGTGAAAACGAAAATTCCAt                            >  1:4303667/1‑36 (MQ=255)
         tACAACAAGAAGTTTGGGTGAAAACGAAAATTCCAt                            <  1:4469288/36‑1 (MQ=255)
           caacaaGAAGTTTGGGTGAAAACGAAAATTCCATGc                          <  1:1355878/36‑1 (MQ=255)
           caacaaGAAGTTTGGGTGAAAACGAAAATTCCATGc                          <  1:4222208/36‑1 (MQ=255)
            aacaaGAAGTTTGGGTGAAAACGAAAATTCCATGCa                         >  1:2199152/1‑36 (MQ=255)
            aacaaGAAGTTTGGGTGAAAACGAAAATTCCATGCa                         >  1:4655319/1‑36 (MQ=255)
             acaaGAAGTTTGGGTGAAAACGAAAATTCCATGCaa                        >  1:2585242/1‑36 (MQ=255)
              caaGAAGTTTGGGTGAAAACGAAAATTCCATGCaaa                       >  1:3007136/1‑36 (MQ=255)
               aagaagTTTGGGTGAAAACGAAAATTCCATGCaaaa                      >  1:2373749/1‑36 (MQ=255)
               aagaagTTTGGGTGAAAACGAAAATTCCATGCaaaa                      >  1:5906335/1‑36 (MQ=255)
               aagaagTTTGGGTGAAAACGAAAATTCCATGCaaaa                      <  1:3590860/36‑1 (MQ=255)
                agaagTTTGGGTGAAAACGAAAATTCCATGCAAAAt                     <  1:1827356/36‑1 (MQ=255)
                agaagTTTGGGTGAAAACGAAAATTCCATGCAAAAt                     <  1:1089842/36‑1 (MQ=255)
                 gaagTTTGGGTGAAAACGGAAATTCCATGCAAAATg                    >  1:1314963/1‑36 (MQ=255)
                 gaagTTTGGGTGAAAACGAAAATTCCATGCAAAATg                    >  1:3030158/1‑36 (MQ=255)
                 gaagTTTGGGTGAAAACGAAAATTCCATGCAAAATg                    >  1:1013028/1‑36 (MQ=255)
                   agTTTGGGTGAAAACGAAAATTCCATGCAAAATGCt                  >  1:979710/1‑36 (MQ=255)
                   agTTTGGGTGAAAACGAAAATTCCATGCAAAATGCt                  <  1:2659770/36‑1 (MQ=255)
                    gTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTc                 >  1:951514/1‑36 (MQ=255)
                    gTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTc                 >  1:5479067/1‑36 (MQ=255)
                    gTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTc                 >  1:5587428/1‑36 (MQ=255)
                    gTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTc                 >  1:1680273/1‑36 (MQ=255)
                     tttGGGTGAAAACGAAAATTCCATGCAAAATGCTcc                >  1:5505659/1‑36 (MQ=255)
                     tttGGGTGAAAACGAAAATTCCATGCAAAATGCTcc                >  1:1476258/1‑36 (MQ=255)
                     tttGGGTGAAAACGAAAATTCCATGCAAAATGCTcc                >  1:3811959/1‑36 (MQ=255)
                       tGGGTGAAAACGAAAATTCCATGCAAAATGCTCCgg              >  1:2732589/1‑36 (MQ=255)
                        gggTGAAAACGAAAATTCCATGAAAAATGCTCAGGt             >  1:3568690/1‑36 (MQ=37)
                         ggTGAAAACGAAAATTCCATGCAAAATGCTCCGGtt            <  1:5292355/36‑1 (MQ=255)
                         ggTGAAAACGAAAATTCCATGCAAAATGCTCCGGtt            >  1:665528/1‑36 (MQ=255)
                          gTGAAAACGAAAATTCCATGCAAAATGCTCCGGttt           <  1:3494267/36‑1 (MQ=255)
                          gTGAAAACGAAAATTCCATGCAAAATGCTCCGGttt           >  1:4709387/1‑36 (MQ=255)
                           tGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTc          >  1:4625037/1‑36 (MQ=255)
                            gAAAACGGAAATTCCATGCAAAATGCTCCGGTTTCa         >  1:3986076/1‑36 (MQ=255)
                            gAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCa         >  1:833949/1‑36 (MQ=255)
                            gAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCa         >  1:5261728/1‑36 (MQ=255)
                            gAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCa         >  1:4307663/1‑36 (MQ=255)
                             aaaaCGAAAATTCCATGCAAAATGCTCCGGTTTCAt        >  1:2579791/1‑36 (MQ=255)
                             aaaaCGAAAATTCCATGCAAAATGCTCCGGTTTCAt        <  1:5108744/36‑1 (MQ=255)
                             aaaaCGAAAATTCCATGCAAAATGCTCCGGTTTCAt        >  1:3153451/1‑36 (MQ=255)
                             aaaaCGAAAATTCCATGCAAAATGCTCCGGTTTCAt        >  1:3253187/1‑36 (MQ=255)
                              aaaCGAAAATTCCATGCAAAATGCTCCGGTTTCATg       <  1:4445335/36‑1 (MQ=255)
                              aaaCGAAAATTCCATGCAAAATGCTCCGGTTTCATg       <  1:5454124/36‑1 (MQ=255)
                               aaCGAAAATTCCATGCAAAATGCTCCGGTTTCATgt      >  1:2598278/1‑36 (MQ=255)
                               aaCGAAAATTCCATGCAAAATGCTCCGGTTTCATgt      >  1:3807703/1‑36 (MQ=255)
                               aaCGAAAATTCCATGCAAAATGCTCCGGTTTCATgt      >  1:3693603/1‑36 (MQ=255)
                                aCGAAAATTCCATGCAAAATGCTCCGGTTTCATgtc     <  1:2968643/36‑1 (MQ=255)
                                   aaaaTTCCATGCAAAATGCTCGGGTTTCATgtcgtc  <  1:1811992/36‑1 (MQ=255)
                                   aaaaTTCCATGCAAAATGCTCCGGTTTCATgtcgtc  <  1:2049824/36‑1 (MQ=255)
                                   aaaaTTCCATGCAAAATGCTCCGGTTTCATgtcgtc  >  1:3867700/1‑36 (MQ=255)
                                   aaaaTTCCATGCAAAATGCTCCGGTTTCATgtcgtc  >  1:1268200/1‑36 (MQ=255)
                                   aaaaTTCCACGCAAAATGCTCCGGTTTCATgtcgtc  <  1:806183/36‑1 (MQ=255)
                                   |                                   
GGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTC  >  REL606/3947762‑3947832

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG