Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,735,9940C→G13.2% 5.8 ‑13.4 0.073145ascBcryptic 6‑phospho‑beta‑glucosidase
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (11/0):  ref base (14/47):  total (26/47)
Fisher's exact test strand distribution p-value = 1.47e-06
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CACATGGCGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGTTT  >  REL606/2735959‑2736029
                                     |                                   
cacaTGGCGAGCATCGAATGGC.GGGGA.AACCAGGGc                                     >  1:1045539/1‑36 (MQ=25)
 acaTGGCGAGCATCGAATGGC.GGTGA.AACTGGGGCt                                    <  1:3809306/36‑1 (MQ=255)
   aTGGCGAGCATCGAATGGC.GGTGA.AACTGGGGCTgg                                  <  1:5190139/36‑1 (MQ=255)
   aTGGCGAGCATCGAATGGC.GGGGG.AAATGGGGGTgg                                  >  1:2663627/1‑36 (MQ=16)
   aTGGCGAGCATCGAATGGC.GGGGA.AAATGGGGGTgg                                  >  1:3153511/1‑36 (MQ=25)
    tGGCGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGa                                 <  1:4615166/36‑1 (MQ=255)
    tGGCGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGa                                 <  1:309059/36‑1 (MQ=255)
    tGGCGAGCATCGAATGGC.GGGGA.AAAAGGGGGTggg                                 >  1:2156053/1‑35 (MQ=16)
     ggggAGCATCGAATGAC.GGTGA.AACTGGGGCTGGaa                                <  1:96130/33‑1 (MQ=37)
     ggCGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGaa                                <  1:2218393/36‑1 (MQ=255)
     ggCGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGaa                                <  1:2414174/36‑1 (MQ=255)
      gCGAGCATCGAATGGG.GGTGA.AACTGGGGCTGGaaa                               <  1:5840954/36‑1 (MQ=255)
      gCGAGCATCGAATGGC.GGTGG.AAACGGGGGTgggga                               >  1:2574793/1‑33 (MQ=12)
       cGAGCATCGAATGGG.GGTGA.AACTGGGGCTGGaaaa                              <  1:5186746/36‑1 (MQ=255)
       cGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGaaaa                              <  1:1186745/36‑1 (MQ=255)
        gAGCATCGAATGGC.GGTGG.AACTGGGGGTGGgaaaa                             >  1:1869216/1‑36 (MQ=25)
        gAGCATCGAATGGC.GGTGG.AAATGGGGGTGGgaaaa                             >  1:2974262/1‑36 (MQ=19)
        gAGCATCGAATGGC.GGTGG.AAAAGGGGGTGGgaaaa                             >  1:894135/1‑36 (MQ=16)
        gAGCATCGAATGGC.GGTAA.AACTGGGGGGGGgaaaa                             >  1:626920/1‑36 (MQ=19)
         aGCATCGAATGGC.GGTGAGAATGGGGGGTGGaaaaa                             >  1:5269231/1‑36 (MQ=16)
         aGCATCGAATGGC.GGTGA.AACTGGGGCTGGAAAAAc                            <  1:5097304/36‑1 (MQ=255)
          gCATCGAATGGC.GGGGC.AAAAGGGGCTGGaaaaaaa                           >  1:1503270/1‑34 (MQ=16)
           cATCGAATGGC.GGCGA.AACTGGGGGGGGGAAAACGt                          >  1:4858330/1‑36 (MQ=16)
            aTCGAATGGCGGGGGA.AACCGGGGTGGGAAAAAGGt                          >  1:3321385/1‑36 (MQ=2)
            aTCGAATGGC.GGTGA.AACTGGGGCTGGAAAAACGtt                         <  1:1583146/36‑1 (MQ=255)
            aTCGAATGGC.GGTGA.AACTGGGGCTGGAAAAACGtt                         <  1:4426043/36‑1 (MQ=255)
             tCGAATGGC.GGTGG.AACTGGGGGTGGGAAAACGgtt                        >  1:1841974/1‑36 (MQ=16)
              cGAATGGC.GGTGA.AACTGGGGCTGGAAAAACGtttt                       <  1:5331580/36‑1 (MQ=255)
              cGAATGGC.GGTGA.AACTGGGGCTGGAAAAACGtttt                       <  1:3743109/36‑1 (MQ=255)
                aaTGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCa                     <  1:1477905/36‑1 (MQ=255)
                aaTGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCa                     <  1:3897914/36‑1 (MQ=255)
                  tGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGt                   <  1:5598876/36‑1 (MQ=255)
                  tGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGt                   <  1:1769221/36‑1 (MQ=255)
                  tGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGt                   <  1:501513/36‑1 (MQ=255)
                   ggC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGtt                  <  1:281118/36‑1 (MQ=255)
                   ggC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGtt                  <  1:5273897/36‑1 (MQ=255)
                    gC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTg                 <  1:4441975/36‑1 (MQ=255)
                    gC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTg                 <  1:1790599/36‑1 (MQ=255)
                    gC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGGTg                 >  1:2843954/1‑36 (MQ=255)
                     c.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTgc                <  1:4810524/36‑1 (MQ=255)
                     c.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTgc                <  1:5906418/36‑1 (MQ=255)
                      ccgTGA.AACTGGGGCTGGAAAAGCGTTTTCAGTTgc                <  1:79645/34‑1 (MQ=37)
                       ggTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTgcg               >  1:5043584/1‑36 (MQ=255)
                       ggTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTgcg               >  1:1693302/1‑36 (MQ=255)
                        gTGA.AACTTGGGCTGGAAAAACGTTTTCAGTTgcgc              <  1:1142138/36‑1 (MQ=255)
                        gTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTgcgc              >  1:1271737/1‑36 (MQ=255)
                        gTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTgcgc              >  1:2888334/1‑36 (MQ=255)
                         tGA.AACTGTGGCTGGAAAAGCGTTTTCAGTTGcgcg             <  1:1660651/36‑1 (MQ=25)
                         tGA.AACTGGGGCTGGAAAAACGTTTTCAGTTGcgcg             <  1:4266001/36‑1 (MQ=255)
                           a.aaCTGGGTCTGGAGAAACGTTTTCAGTTGCGCGAt           <  1:37462/36‑1 (MQ=37)
                           a.aaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGAt           >  1:638271/1‑36 (MQ=255)
                           a.aaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGAt           <  1:2384924/36‑1 (MQ=255)
                           a.aaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGAt           >  1:2532993/1‑36 (MQ=255)
                           a.aaCTGGGGCTGGAAAAACGTTTTCAGNNGCGCGAt           >  1:2437183/1‑36 (MQ=255)
                             aaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATg          >  1:3569296/1‑36 (MQ=255)
                             aaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATg          >  1:287523/1‑36 (MQ=255)
                              aCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGa         >  1:2122137/1‑36 (MQ=255)
                              aCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGa         <  1:4791490/36‑1 (MQ=255)
                              aCTGGGGCTGGAAAAACGTTTTCAGATGCGCGATGa         <  1:3646449/36‑1 (MQ=255)
                               cTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGAc        <  1:4224646/36‑1 (MQ=255)
                               cTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGAc        >  1:3207572/1‑36 (MQ=255)
                                tGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGACg       <  1:3780398/36‑1 (MQ=255)
                                gggggCTGGAAAAACGTTTTCAGTTGCGCGATGACg       <  1:3344639/35‑1 (MQ=255)
                                  gggCTGGAAAAACGTTTTCAGTTGCGCGATGACGAg     <  1:3269157/36‑1 (MQ=255)
                                   ggCTGGAAAACCGTTTTCATTTGCGCGATGACGAGt    <  1:2728480/36‑1 (MQ=37)
                                   ggCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGt    <  1:5227466/36‑1 (MQ=255)
                                   ggCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGt    <  1:3182888/36‑1 (MQ=255)
                                   ggCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGt    >  1:1207338/1‑36 (MQ=255)
                                   ggCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGt    <  1:2263815/36‑1 (MQ=255)
                                    gCTTGAAAAACGTTTTCAGTTGCGCGATGACGAGtt   <  1:2987946/36‑1 (MQ=255)
                                    gCTGGAAAACCGTTTTCAGTGGCGCGATGACGAGtt   <  1:2156379/36‑1 (MQ=38)
                                    gCTGGAAAAACGTTTTCAGTTGCGCGATGCCGAGtt   <  1:5567448/36‑1 (MQ=255)
                                    gCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGtt   <  1:5567741/36‑1 (MQ=255)
                                    gCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGtt   <  1:853519/36‑1 (MQ=255)
                                     cTGGAAAAACGTTTTCAGTTGCGCGATGACGAGttt  >  1:4607416/1‑36 (MQ=255)
                                     |                                   
CACATGGCGAGCATCGAATGGC.GGTGA.AACTGGGGCTGGAAAAACGTTTTCAGTTGCGCGATGACGAGTTT  >  REL606/2735959‑2736029

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG