Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,028,8580C→A11.2% 6.3 ‑9.9 0.067758ugdUDP‑glucose 6‑dehydrogenase
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (0/9):  ref base (42/16):  total (42/25)
Fisher's exact test strand distribution p-value = 4.78e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AGTTCGCCAGTAACTGCTTGGTATCTTTCGGCAGACAATAACCACCATAACCAAACGACGGATTGTTGTAA  >  REL606/2028823‑2028893
                                   |                                   
aGTTCGCCAGTAACTGCTTGGTATCTTTCGGCAGAc                                     <  1:6218477/36‑1 (MQ=255)
aGTTCGCCAGTAACTGCTTGGTATCTTTCGGCAGAc                                     >  1:3705534/1‑36 (MQ=255)
aGTTCGCCAGTAACTGCTTGGTATCTTTCGGCAGAc                                     <  1:3699754/36‑1 (MQ=255)
aGTTCGCAAGTAACTGCTTGGAATCTTTCGGCAGAc                                     >  1:292743/1‑36 (MQ=37)
 gTTCGCCAGTAACTGCTTGGTATCTTTCGGCAGACa                                    <  1:5180820/36‑1 (MQ=255)
 gTTCGCCAGTAACTGCTTGGTATCTTTCGGCAGACa                                    <  1:4499938/36‑1 (MQ=255)
 gTTCGCCAGTAACTGCTTGGTATCTTTCGGCAGACa                                    <  1:6266671/36‑1 (MQ=255)
 gTTCGCCAGTAACTGCTTGGTATCTTTCGGCAGACa                                    <  1:3500014/36‑1 (MQ=255)
  ttCGCCAGTAACTGCTTGGTATCTTTCGGCAGACaa                                   >  1:7071076/1‑36 (MQ=255)
   tCGCCAGTAACTGCTTGGTATCTTTCGGCAGACAAt                                  <  1:6784299/36‑1 (MQ=255)
   tCGCCAGTAACTGCTTGGTATCTTTCGGCAGACAAt                                  >  1:2815905/1‑36 (MQ=255)
   tCGCCAGTAACTGCTTGGTATCTTTCGGCAGACAAt                                  >  1:7621456/1‑36 (MQ=255)
   tCGCCAGTAACTGCTTGGTATCTTTCGGCAGACAAt                                  <  1:7796147/36‑1 (MQ=255)
   tCGCCAGTAACTGCTTGGTATCTTTCGGCAGACAAt                                  >  1:1472228/1‑36 (MQ=255)
    cGCCAGTAACTGCTTGGTATCTTTCGGCAGACAATa                                 <  1:3827766/36‑1 (MQ=255)
    cGCCAGTAACTGCTTGGTATCTTTCGGCAGACAATa                                 >  1:1320800/1‑36 (MQ=255)
    cGCCAGTAACTGCTTGGTATCTTTCGGCAGACAATa                                 >  1:571168/1‑36 (MQ=255)
    cGCCAGTAACTGCTTGGTATCTTTCGGCAGACAATa                                 <  1:1156510/36‑1 (MQ=255)
    cGCCAGTAACTGCTTGGTATCTTTCGGCAGACAATa                                 >  1:6987881/1‑36 (MQ=255)
     gCCAGTAACTGCTTGGTATCTTTCGGCAGACAATaa                                <  1:4563608/36‑1 (MQ=255)
     gCCAGTAACTGCTTGGTATCTTTCGGCAGACAATaa                                <  1:7933606/36‑1 (MQ=255)
     gCCAGTAACTGCTTGGTATCTTTCGGCAGACAATaa                                >  1:2100507/1‑36 (MQ=255)
      ccAGTAACTGCTTGGTATCTTTCGGCAGACAATAac                               >  1:6822120/1‑36 (MQ=255)
      ccAGTAACTGCTTGGTATCTTTCGGCAGACAATAac                               <  1:5314840/36‑1 (MQ=255)
       cAGTAACTGCTTGGTATCTTTCGGCAGACAATAacc                              >  1:7626819/1‑36 (MQ=255)
       cAGTAACTGCTTGGTATCTTTCGGCAGACAATAacc                              >  1:7153543/1‑36 (MQ=255)
       cAGTAACTGCTTGGTATCTTTCGGCAGACAATAacc                              <  1:7514428/36‑1 (MQ=255)
        aGTAACTGCTTGGTATCTTTCGGCAGACAATAacca                             >  1:2295782/1‑36 (MQ=255)
         gTAACTGCTTGGTATCTTTCGGCAGACAATAaccac                            <  1:6559923/36‑1 (MQ=255)
           aaCTGCTTGGTATCTTTCGGCAGACAATAAccacca                          >  1:1118770/1‑36 (MQ=255)
           aaCTGCTTGGTATCTTTCGGCAGACAATAAccacca                          <  1:1715688/36‑1 (MQ=255)
            aCTGCTTGGTATCTTTCGGCAGACAATAACCACCAt                         >  1:5874062/1‑36 (MQ=255)
            aCTGCTTGGTATCTTTCGGCAGACAATAACCACCAt                         >  1:6526352/1‑36 (MQ=255)
            aCTGCTTGGTATCTTTCGGCAGACAATAACCACCAt                         >  1:1362939/1‑36 (MQ=255)
            aCTGCTTGGTATCTTTCGGCAGACAATAACCACCAt                         >  1:664583/1‑36 (MQ=255)
            aCTGCTTGGTATCTTTCGGCAGACAATAACCACCAt                         >  1:1743454/1‑36 (MQ=255)
             cTGCTTGGTATCTTTCGGCAGACAATAACCACCATa                        >  1:954579/1‑36 (MQ=255)
             cTGCTTGGTATCTTTCGGCAGACAATAACCACCATa                        >  1:5679778/1‑36 (MQ=255)
             cTGCTTGGTATCTTTCGGAAAAAAATAACCACCATa                        <  1:406105/36‑1 (MQ=21)
              tGCTTGGTATCTTTCGGCAGACAATAACCACCATaa                       >  1:7607724/1‑36 (MQ=255)
              tGCTTGGTATCTTTCGGCAGACAATAACCACCATaa                       >  1:2732520/1‑36 (MQ=255)
               gCTTGGTATCTTTCGGCAGACAATAACCACCATAAc                      >  1:5315263/1‑36 (MQ=255)
               gCTTGGTATCTTCCGGAAGAAAATAACCACCATAAc                      <  1:525173/36‑1 (MQ=25)
               gCTGGGTATCTTTCGGCAGAAAATAACCACCATAAc                      <  1:2861269/36‑1 (MQ=37)
                 ttGGTATCTTTCGGCAGACAATAACCACCATAACCa                    >  1:119998/1‑36 (MQ=255)
                 ttGGTATCTTTCGGCAGACAATAACCACCATAACCa                    >  1:7068065/1‑36 (MQ=255)
                  tGGTATCTTTCGTAAAAAAATAACCACCATAACCaa                   <  1:3746170/36‑1 (MQ=16)
                  tGGTATCTTTCGGCAGACAATAACCACCATAACCaa                   >  1:3583636/1‑36 (MQ=255)
                  tGGTATCTTTCGGAAAACAATAACCACCATAACCaa                   <  1:7687958/36‑1 (MQ=37)
                   ggTATTTTTCGGAAAAAAAAAACCACCATAACCaaa                  <  1:4345158/36‑1 (MQ=16)
                   ggTATCTTTCGGCAGTCAATAACCACCATAACCaaa                  >  1:5477669/1‑36 (MQ=255)
                   ggTATCTTTCGGCAGACAATAACCACCATAACCaaa                  >  1:2395817/1‑36 (MQ=255)
                   ggTATCTTTCGGCAGACAATAACCACCATAACCaaa                  >  1:528279/1‑36 (MQ=255)
                   ggTATCTTTCGGCAGACAATAACCACCATAACCaaa                  <  1:481340/36‑1 (MQ=255)
                    gTATCTTTCGGCAGACAATAACCACCATAACCAAac                 >  1:3022725/1‑36 (MQ=255)
                     tATCTTTCGGCAGACAATAACCACCATAACCAAacg                >  1:2190349/1‑36 (MQ=255)
                      aTCTTTCGGCAGACAATAACCACCATAACCAAacga               >  1:2305411/1‑36 (MQ=255)
                      aTCTTTCGGCAGACAATAACCACCATAACCAAacga               >  1:6190115/1‑36 (MQ=255)
                       tCTTACGGAAAAAAAAAACCACCATAACCAAacgac              <  1:6133644/36‑1 (MQ=16)
                          ttCGGCAGACAATAACCACCATAACCAAACGACGGa           >  1:6986992/1‑36 (MQ=255)
                           tCGGCAGACAATAACCACCATAACCAAACGACGGAt          >  1:7381669/1‑36 (MQ=255)
                           tCGGAAGAAAATAACCACCATAACCAAACGACGGAt          <  1:1271534/36‑1 (MQ=37)
                            cGGCAGACAATAACCACCATAACCAAACGACGGAtt         >  1:6913723/1‑36 (MQ=255)
                            cGGCAGACAATAACCACCATAACCAAACGACGGAtt         >  1:3529278/1‑36 (MQ=255)
                             ggCAGACAATAACCACCATAACCAAACGACGGAttg        >  1:1311693/1‑36 (MQ=255)
                             gacAGAAAAAAACCACCATAACCAAACGACGGAttg        <  1:4176898/34‑1 (MQ=37)
                              gCAGACAATAACCACCATAACCAAACGACGGAttgt       >  1:2250527/1‑36 (MQ=255)
                              gCAGACAATAACCACCATAACCAAACGACGGAttgt       >  1:3978904/1‑36 (MQ=255)
                               aaGAAAATAACCACCATAACCAAACGACGGAttgtt      <  1:1823018/35‑1 (MQ=255)
                                aGACAATAAACACCATAACCAAACGACGGAttgttg     >  1:2552970/1‑36 (MQ=255)
                                  aCAATAACCACCATAACCAAACGACGGATTGTTGTa   >  1:5708206/1‑36 (MQ=255)
                                   cAATAACCACCATAACCAAACGACGGATTGTTGTaa  >  1:7155790/1‑36 (MQ=255)
                                   |                                   
AGTTCGCCAGTAACTGCTTGGTATCTTTCGGCAGACAATAACCACCATAACCAAACGACGGATTGTTGTAA  >  REL606/2028823‑2028893

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG