Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,627,3840C→A 8.4% 1.0 ‑10.5 0.0741591ydfEpredicted protein
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/9):  ref base (48/17):  total (48/26)
Fisher's exact test strand distribution p-value = 2.83e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ACCAAAAAGCGGAACAAAACAGCCCAAATGCGTTACAAAACGAGCCAGAAACGAAACAGCCTGAACCAGTG  >  REL606/1627349‑1627419
                                   |                                   
aCCAAAAAGCGGAACAAAACAGCCCAAATGCGTTAc                                     <  1:1056012/36‑1 (MQ=255)
aCCAAAAAGCGGAACAAAACAGCCCAAATGCGTTAc                                     >  1:5568719/1‑36 (MQ=255)
aCCAAAAAGCGGAACAAAACAGCCCAAATGCGTTAc                                     >  1:2195754/1‑36 (MQ=255)
   aaaaaGCGGAACAAAACAGCCCAAATGCGTTACaaa                                  >  1:1490092/1‑36 (MQ=255)
   aaaaaGCGGAACAAAACAGCCCAAATGCGTTACaaa                                  >  1:3433812/1‑36 (MQ=255)
   aaaaaGCGGAAC.AAACAGCCCAAATGCGTTACaaaa                                 >  1:3724796/1‑36 (MQ=38)
    aaaaGCGGAACAAAACAGCCCAAATGCGTTACaaaa                                 >  1:5576651/1‑36 (MQ=255)
    aaaaGCGGAACAAAACAGCCCAAATGCGTTACaaaa                                 >  1:5895765/1‑36 (MQ=255)
    aaaaGCGGAACAAAACAGCCCAAATGCGTTACaaaa                                 <  1:1508180/36‑1 (MQ=255)
    aaaaGCGGAACAAAACAGCCCAAATGCGTTACaaaa                                 <  1:1229817/36‑1 (MQ=255)
    aaaaGCGAAACAAAACAGCCCAAATGCGTTACaaaa                                 <  1:4988161/36‑1 (MQ=255)
     aaaGCGGAACAAAACAGCCCAAATGCGTTACAAAAc                                >  1:1198314/1‑36 (MQ=255)
      aaGCGGAACAAAACAGCCCAAATGCGTTACAAAACg                               >  1:279643/1‑36 (MQ=255)
      aaGCGGAACAAAACAGCCCAAATGCGTTACAAAACg                               >  1:2806156/1‑36 (MQ=255)
       aGCGGAACAAAACAGCCCAAATGCGTTACAAAACGa                              >  1:1253645/1‑36 (MQ=255)
        gCGGAACAAAACAGCCCAAATGCGTTACAAAACGAg                             <  1:2750790/36‑1 (MQ=255)
        gCGGAACAAAACAGCCCAAATGCGTTACAAAACGAg                             >  1:709268/1‑36 (MQ=255)
        gCGGAACAAAACAGCCCAAATGCGTTACAAAACGAg                             <  1:6587399/36‑1 (MQ=255)
        gCGGAACAAAACAGCCCAAATGCGTTACAAAACGAg                             >  1:3147971/1‑36 (MQ=255)
        gCGGAACAAAACAGCCCAAATGCGTTACAAAACGAg                             <  1:1453068/36‑1 (MQ=255)
         cGGAACAAAACAGCCCAAATGCGTTACAAAACGAGc                            <  1:4712680/36‑1 (MQ=255)
         cGGAACAAAACAGCCCAAATGCGTTACAAAACGAGc                            >  1:1090609/1‑36 (MQ=255)
         cGGAACAAAACAGCCCAAATGCGTTACAAAACGAGc                            <  1:2484871/36‑1 (MQ=255)
          ggAACAAAACAGCCCAAATGCGTTACAAAACGAGcc                           >  1:3076256/1‑36 (MQ=255)
          ggAACAAAACAGCCCAAATGCGTTACAAAACGAGcc                           >  1:3735931/1‑36 (MQ=255)
          ggAACAAAACAGCCCAAATGCGTTACAAAACGAGcc                           >  1:3351307/1‑36 (MQ=255)
          ggAACAAAACAGCCCAAATGCGTTACAAAACGAGcc                           <  1:2191421/36‑1 (MQ=255)
          ggAACAAAACAGCCCAAATGCGTTACAAAACGAGcc                           >  1:1777330/1‑36 (MQ=255)
           gAACAAAACAGCCCAAATGCGTTACAAAACGAGCCa                          >  1:2867653/1‑36 (MQ=255)
           gAACAAAACAGCCCAAATGCGTTACAAAACGAGCCa                          >  1:1068864/1‑36 (MQ=255)
            aaCAAAACAGCCCAAATGCGTTACAAAACGAGCCAg                         >  1:4453564/1‑36 (MQ=255)
            aaCAAAACAGCCCAAATACCTTACAAAACGATCCAg                         >  1:1699031/1‑36 (MQ=255)
             aCAAAACAGCCCAAATGCGTTACAAAACGAGCCAga                        <  1:4184537/36‑1 (MQ=255)
              cAAAACAGCCCAAATGCGTTACAAAACGAGCCAgaa                       >  1:6268620/1‑36 (MQ=255)
              cAAAACAGCCCAAATGCGTTACAAAACGAGCCAgaa                       >  1:7193339/1‑36 (MQ=255)
              aaaaaCAGCCCAAATGCGTTACAAAACGAGCCAgaa                       <  1:3842454/35‑1 (MQ=255)
               aaaaCAGCCCAAATGCGTTACAAAACGAGCCAgaaa                      >  1:559635/1‑36 (MQ=255)
               aaaaCAGCCCAAATGCGTTACAAAACGAGCCAgaaa                      >  1:631986/1‑36 (MQ=255)
               aaaaCAGCCCAAATGCGTTACAAAACGAGCCAgaaa                      >  1:7217759/1‑36 (MQ=255)
               aaaaCAGCCCAAATGCGTTACAAAACGAGCCAgaaa                      <  1:3807535/36‑1 (MQ=255)
               aaaaCAGCCCAAATGCGTTACAAAACGAGCCAgaaa                      >  1:1117645/1‑36 (MQ=255)
                aaaCAGCCCAAATGCGTTACAAAACGAGCCAgaaac                     >  1:829096/1‑36 (MQ=255)
                aaaCAGCCCAAATGCGTTACAAAACGAGCCAgaaac                     >  1:3953848/1‑36 (MQ=255)
                aaaCAGCCCAAATGCGTTACAAAACGAGCCAgaaac                     <  1:3044544/36‑1 (MQ=255)
                aaaCAGCCCAAATGCGGTACAAAACGAGCCAgatac                     >  1:5741277/1‑36 (MQ=37)
                 aaCAGCCCAAATGCGTTACAAAACGAGCCAgaaacg                    >  1:5981753/1‑36 (MQ=255)
                  aCAGCCCAAATGCGTTACAAAACGAGCCAgaaacga                   >  1:1831086/1‑36 (MQ=255)
                   cAGCCCAAATGCGTTACAAAACGAGCCAgaaacgaa                  >  1:4984138/1‑36 (MQ=255)
                   cAGCCCAAATGCGTTACAAAACGAGCCAgaaacgaa                  >  1:4996950/1‑36 (MQ=255)
                   cAGCCCAAATGCGTTACAAAACGAGCCAgaaacgaa                  <  1:1087277/36‑1 (MQ=255)
                   cAGCCCAAATGCGTTACAAAACGAGCCAgaaacgaa                  >  1:680298/1‑36 (MQ=255)
                    aGCCCAAATGCGTTACAAAACGAGCCAgaaacgaaa                 >  1:775105/1‑36 (MQ=255)
                    aGCCCAAATGCGTTACAAAACGAGCCAgaaacgaaa                 >  1:1467987/1‑36 (MQ=255)
                    aGCCCAAATGCGTTACAAAACGAGCCAgaaacgaaa                 <  1:1232270/36‑1 (MQ=255)
                      cccAAATGCGTTACAAAACGAGCCAGAAACGAAACa               >  1:2622535/1‑36 (MQ=255)
                       ccAAATGCGTTACAAAACGAGCCAGAAACGAAACAg              >  1:3687689/1‑36 (MQ=255)
                       ccAAATGCGTTACAAAACGAGCCAGAAACGAAACAg              >  1:679819/1‑36 (MQ=255)
                        cAAATGCGTTACAAAACGAGCCAGAAACGAAACAGc             >  1:6099724/1‑36 (MQ=255)
                        cAAATGCGTTAAAAAACGAGCCAGAAACGAAACAGc             <  1:4216090/36‑1 (MQ=255)
                          aaTGCGTTACAAAACGAGCCAGAAACGAAACAGCCt           >  1:1689554/1‑36 (MQ=255)
                          aaTGCGTTACAAAACGAGCCAGAAACGAAACAGCCt           <  1:5808254/36‑1 (MQ=255)
                          aaTGCGTTAAAAAACGAGCCAGAAACGAAACAGCCt           <  1:484190/36‑1 (MQ=255)
                           aTGCGTTACAAAACGAGCCAGAAACGAAACAGCCTg          >  1:4043050/1‑36 (MQ=255)
                           aTGCGTTACAAAACGAGCCAGAAACGAAACAGCCTg          >  1:4206615/1‑36 (MQ=255)
                           aTGCGTTACAAAACGAGCCAGAAACGAAACAGCCTg          >  1:5095368/1‑36 (MQ=255)
                           aTGCGTTAAAAAACGAGCCAGAAACGAAACAGCCTg          <  1:2401665/36‑1 (MQ=255)
                            tGCGTTAAAAAACGAGCCAGAAACGAAACAGCCTGa         <  1:1429166/36‑1 (MQ=255)
                              cattAAAAAACGAGCCAGAAACGAAACAGCCTTAAc       <  1:5601470/34‑1 (MQ=37)
                              cGTTAAAAAACGAGCCAGAAACGAAACAGCCTGAAc       <  1:2586256/36‑1 (MQ=255)
                                ttACAAAACGAGCGAGAAACGAAACAGCCTGAACCa     >  1:3242728/1‑36 (MQ=255)
                                ttACAAAACGAGCCAGAAACGAAACAGCCTGAACCa     <  1:6014627/36‑1 (MQ=255)
                                ttAAAAAACGAGCCAGAAACGAAACAGCCTGAACCa     <  1:113692/36‑1 (MQ=255)
                                ttAAAAAACGAGCCAGAAACGAAACAGCCTGAACCa     <  1:954092/36‑1 (MQ=255)
                                ttAAAAAACGAGACAGAAACGAAACAGCCTGAACCa     <  1:5055398/36‑1 (MQ=38)
                                 tACAAAACGAGCCAGAAACGAAACAGCCTGAACCAg    >  1:3845831/1‑36 (MQ=255)
                                  aCAAAACGAGCCAGAAACGAAACAGCCTGAACCAGt   >  1:5735935/1‑36 (MQ=255)
                                  aCAAAACGAGCCAGAAACGAAACAGCCTGAACCAGt   >  1:5398173/1‑36 (MQ=255)
                                   cAAAACGAGCCAGAAACGAAACAGCCTGAACCAGTg  >  1:3470350/1‑36 (MQ=255)
                                   |                                   
ACCAAAAAGCGGAACAAAACAGCCCAAATGCGTTACAAAACGAGCCAGAAACGAAACAGCCTGAACCAGTG  >  REL606/1627349‑1627419

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG