Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,735,9850A→G10.7% 4.8 ‑14.8 0.0107136ascBcryptic 6‑phospho‑beta‑glucosidase
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (13/1):  ref base (20/73):  total (33/74)
Fisher's exact test strand distribution p-value = 3.55e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

ATATGATCCCACATGGCGAGCATCGAATGGCGG.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTGCGCGAT  >  REL606/2735950‑2736020
                                     |                                   
atatGATCCCACATGGCGAGCATCGAATGGCGG.GG.a                                     >  1:1488747/1‑36 (MQ=255)
  atGATCCCACATGTCGAGCATCTAATGGCGG.TG.aaa                                   <  1:7459171/36‑1 (MQ=38)
  atGATCCCACATGGCGAGCATCGAATGGCGG.TG.aaa                                   <  1:7529662/36‑1 (MQ=255)
  atGATCCCACATGGCGAGCATCGAATGGCGG.TG.aaa                                   <  1:4342085/36‑1 (MQ=255)
   tgctCCCACATGGCGAGCATCGAATGGCGG.TG.AAAc                                  <  1:5170235/33‑1 (MQ=255)
   tGATCCCACATGGCTAGCATCGAATGGCGG.TG.AAAc                                  <  1:6677890/36‑1 (MQ=255)
   tGATCCCACATGGCGATCATCGAATGGCGG.TG.AAAc                                  <  1:1939918/36‑1 (MQ=255)
   tGATCCCACATGGCGAGCATCGAATGGCGG.TG.AAAc                                  <  1:5545467/36‑1 (MQ=255)
    gATCCCACATGGCGAGCATCGAATGGCGG.TG.GAAcc                                 >  1:7360330/1‑35 (MQ=255)
    gATCCCACATGGCGAGCATCGAATGGCGG.TG.AAACt                                 <  1:5173714/36‑1 (MQ=255)
    gATCCCACATGGCGAGCATCGAATGGCGG.TG.AAACt                                 <  1:6232648/36‑1 (MQ=255)
    gATCCCACATGGCGAGCATCGAATGGCGG.TG.AAACt                                 <  1:2085895/36‑1 (MQ=255)
     aTCCCACATGGCGAGCATCGAATGGCGG.TG.AAACTg                                <  1:1465236/36‑1 (MQ=255)
      tCCCACATGGCGAGCATCGAATGGCGG.TG.AAACTgg                               <  1:5736623/36‑1 (MQ=255)
       cccACATGGCGAGCATCGAATGGCGG.TG.AAACTggg                              <  1:1712621/36‑1 (MQ=255)
        ccACATGGCGAGCATCGAATGGCGG.TG.AAACTgggg                             <  1:6618079/36‑1 (MQ=255)
        ccACATGGCGAGCATCGAATGGCGG.TG.AAACTgggg                             <  1:1486253/36‑1 (MQ=255)
        ccACATGGCGAGCATCGAATGGCGG.TG.AAACTgggg                             <  1:475378/36‑1 (MQ=255)
        ccACATGGCGAGCATCGAATGGCGG.TG.AAACTgggg                             <  1:8385067/36‑1 (MQ=255)
        ccACATGGCGAGCATCGAATGGCGG.GG.GAACTgggg                             >  1:1301547/1‑36 (MQ=37)
         cacaTGGCGAGCATCGAATGGCGG.TG.GAACTggggg                            >  1:5406390/1‑35 (MQ=255)
         cacaTGGCGAGCATCGAATGGCGG.TG.AAACTGTGGc                            <  1:7435561/36‑1 (MQ=255)
         cacaTGGCGAGCATCGAATGGCGG.TG.AAACTGGGGc                            <  1:5561494/36‑1 (MQ=255)
         cacaTGGCGAGCATCGAATGGCGG.TG.AAACTGGGGc                            <  1:3234420/36‑1 (MQ=255)
         cacaTGGCGAGCATCGAATGGCGG.TG.AAAATggggg                            >  1:7551594/1‑35 (MQ=255)
          acaTGGCGAGCATCGAATGGCGG.TG.GAAATgggggg                           >  1:4243182/1‑34 (MQ=37)
          acaTGGCGAGCATCGAATGGCGG.TG.AAACTGGGGCt                           <  1:526716/36‑1 (MQ=255)
          acaTGGCGAGCATCGAATGGCGG.TG.AAACTGGGGCt                           <  1:7478468/36‑1 (MQ=255)
          acaTGGCGAGCATCGAATGGCGG.GG.GAAATgggggt                           >  1:5791412/1‑34 (MQ=21)
           caTGGCGAGCATCGAATGGCGG.TG.GAAATggggggg                          >  1:8230259/1‑33 (MQ=25)
           caTGGCGAGCATCGAATGGCGG.TG.AAACTGGGGCTg                          <  1:7632797/36‑1 (MQ=255)
           caTGGCGAGCATCGAATGGCGG.TG.AAACTGGGGCTg                          <  1:1036269/36‑1 (MQ=255)
           caTGGCGAGCATAGAATGGCGG.TG.AAACGgggggtt                          >  1:4286928/1‑33 (MQ=25)
            aTGGCGAGCATCGAATGGCGG.TG.GAAATGGGGCTgg                         >  1:7650688/1‑36 (MQ=37)
            aTGGCGAGCATCGAATGGCGG.TG.AAACTTGGGCttg                         >  1:2519649/1‑34 (MQ=38)
            aTGGCGAGCATCGAATGGCGG.GG.AAACTGGGGCTgg                         >  1:7066191/1‑36 (MQ=255)
             tGGCGAGCATCGAATGGCGG.TG.AAACTTGGGGTggg                        >  1:6340740/1‑35 (MQ=37)
             tGGCGAGCATCGAATGGCGG.TG.AAACTGGGGGTggg                        >  1:4970005/1‑35 (MQ=255)
             tGGCGAGCATCGAATGGCGG.TG.AAACTGGGGCTGGa                        <  1:1200523/36‑1 (MQ=255)
             tGGCGAGCATCGAATGGCGG.TG.AAACTGGGGCTGGa                        <  1:5992981/36‑1 (MQ=255)
             tGGCGAGCATCGAATGGCGG.TG.AAACTGGGGCTGGa                        <  1:4351067/36‑1 (MQ=255)
             cggCGAGCATCGAATGGCGG.TG.AAACTGGGGCTGGa                        <  1:4614050/35‑1 (MQ=255)
              ggCGAGCATCGAATTGCGG.TG.AAACTGGGGCTGGaa                       <  1:626235/36‑1 (MQ=255)
              ggCGAGCATCGAATGGCGG.TG.GGACTGGGGGTggga                       >  1:4269396/1‑34 (MQ=21)
              ggCGAGCATCGAATGGCGG.TG.GAAATTGGGCGggga                       >  1:3231131/1‑34 (MQ=16)
              ggCGAGCATCGAATGGCGG.TG.AAACTGGGGCTGGaa                       <  1:5148352/36‑1 (MQ=255)
              ggCGAGCATCGAATGGCGG.TG.AAACTGGGGCTGGaa                       <  1:5836506/36‑1 (MQ=255)
              ggCGAACATCGAATGGCGG.TG.AAACTGGGGCTGGaa                       <  1:5216009/36‑1 (MQ=255)
               gCGAGCATCGAATGGCGG.TG.AAACTGGGGCTGGaaa                      <  1:5538425/36‑1 (MQ=255)
               gCGAGCATCGAATGGCGG.TG.AAACTGGGGCTGGaaa                      <  1:231613/36‑1 (MQ=255)
                cGAGCATCGAATGGCGG.TG.AAACTGGGGCTGGaaaa                     <  1:7531947/36‑1 (MQ=255)
                cGAGCATCGAATGGCGG.TG.AAACTGGGGCTGGaaaa                     <  1:1893677/36‑1 (MQ=255)
                 gAGCATCGAATGGCGG.TGAAAACTGGGGGGGGaaaa                     >  1:2118797/1‑36 (MQ=16)
                 gAGCATCGAATGGCGG.TG.AAACTGGGGCTGGaaaaa                    <  1:803035/36‑1 (MQ=255)
                 gAGCATCGAATGGCGG.TG.AAACTGGGGCTGGaaaaa                    <  1:2389999/36‑1 (MQ=255)
                  aGCATCGAATGGCGG.TG.AAACTGGGGCTGGAAAAAc                   <  1:6482206/36‑1 (MQ=255)
                  aGCATCGAATGGCGG.TG.AAACTGGGGCTGGAAAAAc                   <  1:3737854/36‑1 (MQ=255)
                  aGCATCGAATGGCGG.GG.AAACTGGGGCTGGGAAAAc                   >  1:5875267/1‑36 (MQ=37)
                  aGCATCGAATGGCGG.CG.GAACTGGGGGTGGaacaaa                   >  1:924790/1‑35 (MQ=16)
                   gCATCGAATGGCGGCGG.AAATTAGGGGTGGgaaaaa                   >  1:5120213/1‑35 (MQ=0)
                   gCATCGAATGGCGG.TG.AAACTGGGGGTGGgaaaaag                  >  1:2459014/1‑34 (MQ=25)
                   gCATCGAATGGCGG.TG.AAACTGGGGCTGGAAAAACg                  <  1:4555302/36‑1 (MQ=255)
                   gCATCGAATGGCGG.TG.AAACTGGGGCTGGAAAAACg                  <  1:4903341/36‑1 (MQ=255)
                    cATCGAATGGCGG.TG.AAACTGGGGCTGGAAAAACGt                 <  1:4749623/36‑1 (MQ=255)
                      tctaaTGGCGG.TG.AAACTGGGGCTGGAAAAACGttt               <  1:6748780/33‑1 (MQ=255)
                      tCGAATGGCGG.TG.AAACTGGGTCTGGAAAAACGttt               <  1:1075479/36‑1 (MQ=255)
                      tCGAATGGCGG.TG.AAACTGGGGCTGGAAAAACGttt               <  1:1243439/36‑1 (MQ=255)
                      tCGAATGGCGG.TG.AAACTGGGGCTGGAAAAACGttt               <  1:2534714/36‑1 (MQ=255)
                      tCGAATGGCGG.TG.AAACTGGGGCTGGAAAAACGttt               <  1:3120013/36‑1 (MQ=255)
                       cGAATGGCGG.TG.AAACTGGGGCTGGAAAAACGtttt              <  1:8334566/36‑1 (MQ=255)
                        gAATGGCGG.TG.AAACTGGGGCTGGAAAAACGTTTTc             <  1:7097397/36‑1 (MQ=255)
                          aTGGCGG.TGAAAAATGGGGTGGGAAAAAATtttttg            >  1:3807395/1‑34 (MQ=2)
                          aTGGCGG.TG.GAACTGGGGCTGGAAAAACGTTTTCAg           <  1:3362099/36‑1 (MQ=255)
                          aTGGCGG.TG.AAACAGGGGGTGGGAAAAAGTTTTTAg           >  1:20032/1‑36 (MQ=16)
                             gCGG.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTg        <  1:1882446/36‑1 (MQ=255)
                             gCGG.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTg        <  1:3902403/36‑1 (MQ=255)
                             gCGG.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTg        <  1:4251907/36‑1 (MQ=255)
                             gCGG.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTg        <  1:1228703/36‑1 (MQ=255)
                             gCGG.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTg        >  1:4657050/1‑36 (MQ=255)
                             gCGG.TG.AAACTAGGGCTGGAAAAACGTTTTCAGTTg        <  1:2593161/36‑1 (MQ=255)
                              cGG.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTgc       <  1:3238496/36‑1 (MQ=255)
                              cGG.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTgc       <  1:6061069/36‑1 (MQ=255)
                              cGG.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTgc       <  1:4339652/36‑1 (MQ=255)
                               gg.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTgcg      <  1:8451297/36‑1 (MQ=255)
                               gg.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTgcg      >  1:610725/1‑36 (MQ=255)
                               gg.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTgcg      <  1:4470762/36‑1 (MQ=255)
                               gg.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTgcg      <  1:4100827/36‑1 (MQ=255)
                               gg.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTgcg      <  1:2227285/36‑1 (MQ=255)
                               gg.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTgcg      >  1:1025065/1‑36 (MQ=255)
                                g.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTgcgc     <  1:2813153/36‑1 (MQ=255)
                                g.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTgcgc     >  1:2712371/1‑36 (MQ=255)
                                  tG.AAACTGGGGCTGGAAAAACGTTTTCAGTTGcgcg    <  1:965159/36‑1 (MQ=255)
                                   g.AAACTGGGGCTGGAAAAACGTTTTCAGTTGCGCGa   <  1:7127807/36‑1 (MQ=255)
                                   g.AAACTGGGGCTGGAAAAACGTTTTCAGTTGCGCGa   >  1:4340649/1‑36 (MQ=255)
                                   g.AAACTGGGGCTGGAAAAACGTTTTCAGTTGCGCGa   >  1:2897027/1‑36 (MQ=255)
                                   g.AAACTGGGGCTGGAAAAACGTTTTCAGTTGCGCGa   >  1:2281042/1‑36 (MQ=255)
                                     aaaCTGGGGCTGGAAAAACGTTTTCAGTTTCGCGAt  <  1:7651458/36‑1 (MQ=255)
                                     aaaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGAt  >  1:6997752/1‑36 (MQ=255)
                                     aaaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGAt  <  1:5195674/36‑1 (MQ=255)
                                     aaaCTGGGGCTGGAAAAACGTTTTCAGTTGCGCGAt  <  1:3713111/36‑1 (MQ=255)
                                     |                                   
ATATGATCCCACATGGCGAGCATCGAATGGCGG.TG.AAACTGGGGCTGGAAAAACGTTTTCAGTTGCGCGAT  >  REL606/2735950‑2736020

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG