Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,820,0050G→T 8.8% 0.3 ‑10.1 0.071112xniexonuclease IX
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (9/0):  ref base (17/45):  total (26/45)
Fisher's exact test strand distribution p-value = 4.20e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TAAGGGGGTTATCGTGGCTGTTCATTTGCTTATTGTCGATGCACTGAATCTTATTCGTCGCATTCATGCC  >  REL606/2819971‑2820040
                                  |                                   
tAAGGGGGTTATCGTGGCTGTTCATTTGCTTATTGt                                    <  1:5018069/36‑1 (MQ=255)
 aaGGGGGTTATCGTGGCTGTTCATTTTCTTATTTTc                                   >  1:7454605/1‑36 (MQ=37)
  aGGGGGTTATCGTGGCTGTTCATTTTCTTATTTTCg                                  >  1:5809521/1‑36 (MQ=37)
  aGGGGGTTATCGTGGCTGTTCATTTGCTTATTGTCg                                  <  1:2547422/36‑1 (MQ=255)
  aGGGGGTTATCGTGGCTGGTCATTTGCTTTTTTTCg                                  >  1:2373467/1‑36 (MQ=25)
   gggggTTATCGTGGCTGTTTATTTTCTTATTTTCGa                                 >  1:2227026/1‑36 (MQ=25)
   gggggTTATCGTGGCTGTTCATTTTCTTATTTTCGa                                 >  1:5109467/1‑36 (MQ=37)
   gggggTTATCGTGGCTGTTCATTTGCTTATTGTCGa                                 <  1:3200776/36‑1 (MQ=255)
   gggggTTATCGTGGCTGTTCATTTGCTTATTGTCGa                                 <  1:7042485/36‑1 (MQ=255)
   gggggTTATCGTGGCTGTTCATTTGCTTATTGTCGa                                 <  1:3890510/36‑1 (MQ=255)
    ggggTTATCGTGGCTGTTCATTTGCTTATTGTCGAt                                <  1:693113/36‑1 (MQ=255)
    ggggTTATCGTGGCTGTTCATTTGCTTATTGTCGAt                                <  1:6203791/36‑1 (MQ=255)
    ggggTTATCGTGGCTGTTCATTTGCTTATTGTCGAt                                <  1:3050047/36‑1 (MQ=255)
    ggggTTATCGTGGCTGTTCATTTGCTTATTGTCGAt                                <  1:2492966/36‑1 (MQ=255)
     gggTTATCGTGGCTGTTCATTTTCTTTTTTTCTATg                               >  1:3819801/1‑36 (MQ=16)
      ggTTATCGTGGCTGTTCATTTGCTTATTGTCGATGc                              <  1:2340031/36‑1 (MQ=255)
      ggTTATCGTGGCTGTTCATTTGCTTATTGTCGATGc                              <  1:986799/36‑1 (MQ=255)
        ttATCGTGGCTGTTCCTTTTCTTATTGTCTATGCAc                            >  1:2676757/1‑36 (MQ=25)
        ttATCGTGGCTGTTCATTTGCTTATTTTCGATGCAc                            >  1:6954818/1‑36 (MQ=255)
        ttATCGTGGCTGTTCATTTGCTTATTGTCGATGCAc                            <  1:4876149/36‑1 (MQ=255)
        ttATCGTGGCTGTACATTTGCTTATTGTCGATGCAc                            <  1:5752922/36‑1 (MQ=255)
         tATCGTGGCTGTTGATTTGCTTATTGTCGATGCACt                           <  1:1182130/36‑1 (MQ=255)
         tATCGTGGCTGTTCATTTGCTTATTGTCGATGCACt                           <  1:486129/36‑1 (MQ=255)
            cGTGGCTGTTCATTTGCTTATTGTCGATGCACTGaa                        <  1:7186077/36‑1 (MQ=255)
            cGTGGCTGTTCATTTGCTTATTGTCGATGCACTGaa                        <  1:5290463/36‑1 (MQ=255)
             gTGGCTGTTCATTTGCTTATTGTCGATGCACTGAAt                       <  1:2419563/36‑1 (MQ=255)
              ttgCTGTTCATTTGCTTATTGTCGATGCACTGAATc                      <  1:5360055/34‑1 (MQ=255)
              tGGCTGTTCATTTGCTTATTGTCGATGCACTGAATc                      >  1:6119075/1‑36 (MQ=255)
              tGGCTGTT.ATTTGCTTATTTTCGATGCACTGAAttt                     >  1:3544009/1‑34 (MQ=25)
                gCTGTTCATTTGCTTATTTTCGATGCACTGAATCtt                    >  1:1783204/1‑36 (MQ=255)
                gCTGTTCATTTGCTTATTGTCGATGCACTGAATCtt                    >  1:5112134/1‑36 (MQ=255)
                 cTGTTCATTTGCTTATTGTCGATGCACTGAATCTTa                   <  1:5930513/36‑1 (MQ=255)
                  tGTTCATTTGCTTATTGTCGATGCACTGAATCTTAt                  >  1:7485378/1‑36 (MQ=255)
                  tGTTCATTTGCTTATTGTCGATGCACTGAATCTTAt                  >  1:7496774/1‑36 (MQ=255)
                  tGTTCATTTGCTTATTGTCGATGCACTGAATCTTAt                  >  1:2519378/1‑36 (MQ=255)
                  tGTTCATTTGCTTATTGTCGATGCACTGAATCTTAt                  <  1:5494208/36‑1 (MQ=255)
                   gTTCATTTGCTTATTGTCGATGCACTGAATCTTAtt                 <  1:3070555/36‑1 (MQ=255)
                    ttCATTTGCTTATTGTCGATGCACTGAATCTTATtc                >  1:6177967/1‑36 (MQ=255)
                     tCATTTGCTTATTGTCGATGCACTGAATCTTATtcg               >  1:7111843/1‑36 (MQ=255)
                     tCATTTGCTTATTGTCGATGCACTGAATCTTATtcg               >  1:5144738/1‑36 (MQ=255)
                      cATTTGCTTATTGTCGATGCACTGAATCTTATtcgt              >  1:6104340/1‑36 (MQ=255)
                      cATTTGCTTATTGTCGATGCACTGAATCTTATtcgt              >  1:4997538/1‑36 (MQ=255)
                      cATTTGCTTATTGTCGATGCACTGAATCTTATtcgt              >  1:6965929/1‑36 (MQ=255)
                      cATTTGCTTATTGTCGATGCACTGAATCTTATtcgt              <  1:3117863/36‑1 (MQ=255)
                       aTTTGCTTATTGTCGATGCACTGAATCTTATtcgtc             <  1:4962846/36‑1 (MQ=255)
                       aTTTGCTTATTGTCGATGCACTGAATCTTATtcgtc             <  1:170507/36‑1 (MQ=255)
                       aTTTGCTTATTGTCGATGCACTGAATCTTATtcgtc             <  1:7086423/36‑1 (MQ=255)
                       aTTTGCTTATTGTCGATGCACTGAATCTTATtcgtc             >  1:7455184/1‑36 (MQ=255)
                        tttGCTTATTGTCGATGCACTGAATCTTATtcgtcg            <  1:1590627/36‑1 (MQ=255)
                         ttGCTTATTGTCGATGCACTGAATCTTATTCGTCGc           <  1:3559815/36‑1 (MQ=255)
                         ttGCTTATTGTCGATGCACTGAATCTTATTCGTCGc           <  1:4302050/36‑1 (MQ=255)
                         ttGCTTATTGTCGATGCACTGAATCTTATTCGTCGc           <  1:1853419/36‑1 (MQ=255)
                           gCTTATTGTCGATGCACTGAATCTTATTCGTCGCAt         <  1:622845/36‑1 (MQ=255)
                           gCTTATTGTCGATGCACTGAATCTTATTCGTCGCAt         <  1:4423535/36‑1 (MQ=255)
                            cTTATTGTCGATGCACTGAATCTTATTCGTCGCAtt        <  1:7495999/36‑1 (MQ=255)
                            cTTATTGTCGATGCACTGAATCTTATTCGTCGCAtt        >  1:6746951/1‑36 (MQ=255)
                            cTTATTGTCGATGCACTGAATCTTATTCGTCGCAtt        <  1:6987242/36‑1 (MQ=255)
                             ttATTGTCGATGCACTGAATCTTATTCGTCGCATTc       <  1:6547328/36‑1 (MQ=255)
                             ttATTGTCGATGCACTGAATCTTATTCGTCGCATTc       <  1:1634214/36‑1 (MQ=255)
                             ttATTGTCGATGCACTGAATCTTATTCGTCGCATTc       <  1:790700/36‑1 (MQ=255)
                             ttATTGTCGATGCACTGAATCTTATTCGTCGCATTc       <  1:3982757/36‑1 (MQ=255)
                             ttATTGTCGATGCACTGAATCTTATTCGTCGCATTc       <  1:3435253/36‑1 (MQ=255)
                              tATTGTCGATGCACTGAATCTTATTCGTCGCATTCa      <  1:2533932/36‑1 (MQ=255)
                              tATTGTCGATGCACTGAATCTTATTCGTCGCATTCa      <  1:1153520/36‑1 (MQ=255)
                              tATTGTCGATGCACTGAATCTTATTCGTCGCATTCa      >  1:3389327/1‑36 (MQ=255)
                               aTTGTCGATGTTCTGAATCTTATTCGTCGCATTCAt     <  1:4050406/36‑1 (MQ=37)
                               aTTGTCGATGCACTGAATCTTATTCGTCGCATTCAt     >  1:1051239/1‑36 (MQ=255)
                                ttGTCGATGCACTGAATCTTATTCGTCGCATTCATg    >  1:3369093/1‑36 (MQ=255)
                                ttGTCGATGCACTGAATCTTATTCGTCGCATTCATg    <  1:1234438/36‑1 (MQ=255)
                                 tGTCGATGCACTGAATCTTATTCGTCGCATTCATGc   <  1:6393353/36‑1 (MQ=255)
                                 tGTCGATGCACTGAATCTTATTCGTCGCATTCATGc   <  1:999028/36‑1 (MQ=255)
                                  gTCGATGCACTGAATCTTATTCGTCGCATTCATGcc  >  1:7167725/1‑36 (MQ=255)
                                  gTCGATGCACTGAATCTTATTCGTCGCATTCATGcc  <  1:859790/36‑1 (MQ=255)
                                  |                                   
TAAGGGGGTTATCGTGGCTGTTCATTTGCTTATTGTCGATGCACTGAATCTTATTCGTCGCATTCATGCC  >  REL606/2819971‑2820040

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG