Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,201,9490G→A 4.6% 0.9 ‑0.5 ‑0.165611iclRDNA‑binding transcriptional repressor
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (1/2):  ref base (35/27):  total (36/29)
Fisher's exact test strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.47e-01

TTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGAC  >  REL606/4201915‑4201984
                                  |                                   
ttCCTCATCGTCAAATGAATAACCCCGTTTGCGCGt                                    <  1:5973208/36‑1 (MQ=255)
ttCCTCATCGTCAAATGAATAACCCCGTTTGCGCGt                                    <  1:3970808/36‑1 (MQ=255)
 tCCTCATCGTCAAATGAATAACCCCGTTTGCGCGtt                                   <  1:3963820/36‑1 (MQ=255)
  ccTCATCGTCAAATGAATAACCCCGTTTGCGCGttt                                  <  1:5564362/36‑1 (MQ=255)
  ccTCATCGTCAAATGAATAACCCCGTTTGCGCGttt                                  <  1:2057925/36‑1 (MQ=255)
   cTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTg                                 >  1:6673844/1‑36 (MQ=255)
   cTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTg                                 <  1:6650682/36‑1 (MQ=255)
   cTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTg                                 <  1:6497371/36‑1 (MQ=255)
    tCATCGTCAAATGAATAACCCCGTTTGCGCGTTTgg                                <  1:1092980/36‑1 (MQ=255)
     cATCGTCAAATGAATAACCCCGTTTGCGCGTTTggg                               >  1:5144464/1‑36 (MQ=255)
      aTCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGc                              >  1:6051900/1‑36 (MQ=255)
      aTCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGc                              >  1:635228/1‑36 (MQ=255)
      aTCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGc                              <  1:2631010/36‑1 (MQ=255)
      aTCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGc                              <  1:2528697/36‑1 (MQ=255)
      aTCGTCAAATGAATAACCCCGTTTGCGCATTTGGGc                              <  1:7176346/36‑1 (MQ=255)
        cGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCga                            >  1:7276425/1‑36 (MQ=255)
         gTCAAATGAATAACCCCGTTTGCGCGTTTGGGCgag                           >  1:6850678/1‑36 (MQ=255)
         gTCAAATGAATAACCCCGTTTGCGCGTTTGGGCgag                           >  1:2597838/1‑36 (MQ=255)
          tCAAATGAATAACCCCGTTTGCGCGTTTGGGCgaga                          >  1:6101835/1‑36 (MQ=255)
          tCAAATGAATAACCCCGTTTGCGCGTTTGGGCgaga                          >  1:6969308/1‑36 (MQ=255)
          tCAAATGAATAACCCCGTTTGCGCGTTTGGGCgaga                          <  1:3766034/36‑1 (MQ=255)
           cAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGAt                         >  1:2753704/1‑36 (MQ=255)
           cAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGAt                         <  1:553847/36‑1 (MQ=255)
           cAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGAt                         <  1:5404249/36‑1 (MQ=255)
           cAAATGAATAACCCCGTTTGAGCGTTTGGGCGAGAt                         >  1:6961299/1‑36 (MQ=255)
            aaaTGAATAACCCCGTTTGCGCGTTTTGGCGAGAtc                        <  1:3823521/36‑1 (MQ=255)
            aaaTGAATAACCCCGTTTGCGCGTTTGGGCGAGAtc                        >  1:1468701/1‑36 (MQ=255)
             aaTGAATAACCCCGTTTGCGCGTTTGGGCGAGAtct                       >  1:3276533/1‑36 (MQ=255)
             aaTGAATAACCCCGTTTGCGCGTTTGGGCGAGAtct                       >  1:1471927/1‑36 (MQ=255)
              aTGAATAACCCCGTTTGCGCGTTTGGGCGAGAtctt                      >  1:4295551/1‑36 (MQ=255)
              aTGAATAACCCCGTTTGCGCGTTTGGGCGAGAtctt                      >  1:632328/1‑36 (MQ=255)
              aTGAATAACCCCGTTTGCGCGTTTGGGCGAGAtctt                      >  1:1477813/1‑36 (MQ=255)
              aTGAATAACCCCGCTTGCGAGTTTGGGCGAGAtctt                      >  1:2652046/1‑36 (MQ=37)
               tGAATAACCCCGTTTGCGCGTTTGGGCGAGAtcttc                     >  1:6659288/1‑36 (MQ=255)
                gAATAACCCCGTTTGCGCGTTTGGGCGAGAtcttct                    >  1:225544/1‑36 (MQ=255)
                 aaTAACCCCGTTTGCGCGTTTGGGCGAGATcttctt                   <  1:288061/36‑1 (MQ=255)
                  aTACCCCCGTTTGCGCGTTTGGGCGAGATCTTCttt                  <  1:6064980/36‑1 (MQ=255)
                  aTAACCCCGTTTGCGCGTTTGGGCGGGATCTTCttt                  <  1:5131242/36‑1 (MQ=255)
                  aTAACCCCGTTTGCGCGTTTGGGCGAGATCTTCttt                  <  1:2912403/36‑1 (MQ=255)
                  aTAACCCCGTTTGCGCGTTTGGGCGAGATCTTCttt                  <  1:274271/36‑1 (MQ=255)
                  aTAACCCCGTTTGCGCGATTGGGCGAGATCTTCttt                  >  1:2140367/1‑36 (MQ=255)
                  aTAACCCCGTTTGCGCATTTGGGCGAGATCTTCttt                  <  1:7226467/36‑1 (MQ=255)
                   tAACCCCGTTTGCGCGTTTGGGCGAGATCTTCtttt                 >  1:2783735/1‑36 (MQ=255)
                     aCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTaa               >  1:4034785/1‑36 (MQ=255)
                     aCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTaa               >  1:3571201/1‑36 (MQ=255)
                     aCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTaa               <  1:3303667/36‑1 (MQ=255)
                       cccGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAAt             >  1:5553947/1‑36 (MQ=255)
                        ccGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATg            <  1:324623/36‑1 (MQ=255)
                         cGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGc           >  1:4914211/1‑36 (MQ=255)
                         cGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGc           >  1:4496631/1‑36 (MQ=255)
                         cGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGc           <  1:1572129/36‑1 (MQ=255)
                          gTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGca          <  1:668803/36‑1 (MQ=255)
                           tttGCGCGTTTGGGCGAGATCTTCTTTTAAATGcac         <  1:6866736/36‑1 (MQ=255)
                           tttGCGCGTTTGGGCGAGATCTTCTTTTAAATGcac         >  1:2903892/1‑36 (MQ=255)
                           tttGCGCGTTTGGGCGAGATCTTCTTTTAAATGcac         >  1:4171486/1‑36 (MQ=255)
                           tttGCGCGTTTGGGCGAGATCTTCTTTTAAATGcac         <  1:6625252/36‑1 (MQ=255)
                            ttGCGCGTTTGGGCGAGATCTTCTTTTAAATGcaca        <  1:1729500/36‑1 (MQ=255)
                            ttGCGCGTTTGGGCGAGATCTTCTTTTAAATGcaca        <  1:1665772/36‑1 (MQ=255)
                            ttGCGCGTTTGGGCGAGATCTTCTTTTAAATGcaca        >  1:3389059/1‑36 (MQ=255)
                            ttGCGCGTTTGGGCGAGATCTTCTTTTAAATGcaca        >  1:1506248/1‑36 (MQ=255)
                             tGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAg       >  1:3317562/1‑36 (MQ=255)
                              gcgcGTTTGGGCGAGATCTTCTTTTAAATGCACAgg      >  1:2294869/1‑36 (MQ=255)
                              gcgcGTTTGGGCGAGATCTTCTTTTAAATGCACAgg      >  1:2010491/1‑36 (MQ=255)
                              gcgcGGTTGGGCGAGATCTTCTTTTAAATGCACAgg      >  1:2698344/1‑36 (MQ=255)
                              gcgcATTTGGGCGAGATCTTCTTTTAAATGCACAgg      >  1:5847048/1‑36 (MQ=255)
                               cgcgTTTGGGCGAGATCTTCTTTTAAATGCACAGga     >  1:1527408/1‑36 (MQ=255)
                               cgcgTTTGGGCGAGATCTTCTTTTAAATGCACAGga     <  1:4571004/36‑1 (MQ=255)
                                gcgTTTGGGCGAGATCTTCTTTTAAATGCACAGgag    >  1:3356780/1‑36 (MQ=255)
                                gcgTTTGGGCGAGATCTTCTTTTAAATGCACAGgag    <  1:6810453/36‑1 (MQ=255)
                                  gTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGac  >  1:5426999/1‑36 (MQ=255)
                                  |                                   
TTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGAC  >  REL606/4201915‑4201984

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