Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,383,4600G→T14.0% 5.4 ‑10.3 0.079619purAadenylosuccinate synthetase
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (13/1):  ref base (21/44):  total (34/45)
Fisher's exact test strand distribution p-value = 3.47e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CGACATCCTGACTTCTATGGTGGTTGACGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGTGGCGATTTC  >  REL606/4383426‑4383495
                                  |                                   
cGACATCCTGACTTCTATGGTGGTTGACGTTTCTGa                                    >  1:3457030/1‑36 (MQ=255)
 gccATCCTTACTTCTATGGTGGTGGACGTTTCTGAc                                   <  1:5729076/34‑1 (MQ=25)
 gACATCCTGACTTCTATGGTGGTTTACGTTTCTTAc                                   >  1:1698343/1‑36 (MQ=37)
 gACATCCTGACTTCTATGGTGGTTGACGTTTCTGAc                                   <  1:6146197/36‑1 (MQ=255)
 gACATCCTGACTTCTATGGTGGTTGACGTTTCTGAc                                   <  1:675494/36‑1 (MQ=255)
  aCATCCTGACTTCTATGGTGGTTGACGTTTCTGAcc                                  <  1:7045377/36‑1 (MQ=255)
  aCATCCTGACTTCTATGGTGGTTGACGTTTCTGAcc                                  <  1:1609153/36‑1 (MQ=255)
   cATCCTGACTTCTATGGTGGTTGACGTTTCTTACCt                                 >  1:2751130/1‑36 (MQ=255)
   cATCCTGACTTCTATGGTGGTTGACGTTTCTGACCt                                 <  1:3620034/36‑1 (MQ=255)
   cATCCTGACTTCTATGGTGGTTGACGTTTCTGACCt                                 <  1:7181994/36‑1 (MQ=255)
    aTCCTGACTTCTATGTTGGTTAACGTTTCTTACCTg                                >  1:7099339/1‑36 (MQ=25)
    aTCCTGACTTCTATGGTGGTTGACGTTTCTTACCtt                                >  1:5201690/1‑35 (MQ=255)
    aTCCTGACTTCTATGGTGGTTGACGTTTCTGACCTg                                <  1:2585380/36‑1 (MQ=255)
    aTCCTGACTTCTATGGTGGTTGACGTTTCTGACCTg                                <  1:2353408/36‑1 (MQ=255)
     tCCTGTCTTCTATGGCGATTGACGTTTCTGACCTGc                               <  1:4141843/36‑1 (MQ=37)
     tCCTGACTTCTATGGTGTTTGACGTTTCTGACCTGc                               <  1:6813200/36‑1 (MQ=255)
     tCCTGACTTCTATGGTGGTTGACGTTTCTGACCTGc                               <  1:4560266/36‑1 (MQ=255)
      ccTGACTTCTATGGTGGTTTACGTTTTTGACCTTCt                              >  1:4769278/1‑36 (MQ=25)
      ccTGACTTCTATGGTGGTTGACGTTTCTTACCTGCt                              >  1:1795988/1‑36 (MQ=255)
      ccTGACTTCTATGGTGGTTGACGTTTCTGACCTGCt                              <  1:364229/36‑1 (MQ=255)
      ccTGACTTCTATGGTGGTTGACGTTTCTGACCTGCt                              >  1:1486418/1‑36 (MQ=255)
        tGACTTCTATGGTGGTTTACGTTTCTTACCTTCTCg                            >  1:3379785/1‑36 (MQ=25)
        tGACTTCTATGGTGGTTGACGTTTCTGACCTGCTCg                            <  1:5355322/36‑1 (MQ=255)
        tGACTTCTATGGCGGTTCGCATTTCTTTCCAGCTCg                            >  1:6761408/1‑36 (MQ=2)
         gACTTCTATGGTGGTTGACGTTTCTTACCTGCTCGa                           >  1:2088116/1‑36 (MQ=255)
          aCTTCTATGGTGGTTGATGTTTCTTACCTGCTCGAc                          >  1:3244126/1‑36 (MQ=25)
          aCTTCTATGGTGGTTGACGTTTCTTACCTGCTCGAc                          >  1:101316/1‑36 (MQ=255)
          aCTTCTATGGTGGTTGACGTTTCTGACCTGCTCGAc                          <  1:7673674/36‑1 (MQ=255)
          aCTTCTATGGTGGTTGACGTTTCTGACCTGCTCGAc                          >  1:387483/1‑36 (MQ=255)
           cTTCTATGGTTGTTGACGTTTCTGACCTGCTCGAcc                         <  1:3573596/36‑1 (MQ=255)
           cTTCTATGGTGTTTTACGTTTCTGCCCTGCTCGAcc                         <  1:3267544/36‑1 (MQ=25)
           cTTCTATGGTGGTTGACGTTTCTGACCTGCTCGAcc                         <  1:5713877/36‑1 (MQ=255)
           cTTCTATGGTGGTTGACGTTTCTGACCTGCTCGAcc                         <  1:7665601/36‑1 (MQ=255)
           cTTCTATAGTGGTTGACGTTTCTGACCTGCTCGAcc                         <  1:1387872/36‑1 (MQ=255)
            ttCTATGGTGGTTGACGTTTCTGACCTGCTCGACCa                        <  1:4346709/36‑1 (MQ=255)
             tCTATGGTTGTTGACGTTTCTGACCTGCTCGACCAg                       <  1:7634777/36‑1 (MQ=255)
             tCTATGGTGGTTGGCGTTTCTTAGCTGCTCTAGCAg                       <  1:7169921/36‑1 (MQ=16)
             tCTATGGTGGTTGACGTTTCTTACCTTCTCGACCAg                       >  1:1786464/1‑36 (MQ=25)
             tCTATGGTGGTTGACGTTTCTTACCTGCTCGACCAg                       >  1:3172156/1‑36 (MQ=255)
             tCTATGGTGGTTGACGTTTCTGACCTGCTCGACCAg                       <  1:3280531/36‑1 (MQ=255)
               tATGGTGGTTGACGTTTCTGACCTGCTCGACCAGgc                     <  1:1545259/36‑1 (MQ=255)
                aTGGTGGTTGACGTTTCTTACCTGCTCGACCAGgcg                    >  1:3630084/1‑36 (MQ=255)
                aTGGTGGTTGACGTTTCTGACCTGCTCGACCAGgcg                    <  1:7535745/36‑1 (MQ=255)
                 tggtggTTGACGTTTCTGACCTGCTCGACCAGgcgc                   <  1:7310679/36‑1 (MQ=255)
                  ggtggtTGACGTTTCTGACCTGCTCGACCAGGcgcg                  >  1:4072174/1‑36 (MQ=255)
                   gtggtTGACGTTTCTGACCTTCTCGACCAGGCGCGt                 >  1:93401/1‑36 (MQ=255)
                   gtggtTGACGTTTCTGACCTGCTCGACCAGGCGCGt                 >  1:465118/1‑36 (MQ=255)
                   gtggtTGACGTTTCTGACCTGCTCGACCAGGCGCGt                 >  1:1929643/1‑36 (MQ=255)
                    tggtTGACGTTTCTGACCTGCTCGACCAGGCGCGTc                <  1:4737717/36‑1 (MQ=255)
                    tggtTGACGTTTCTGACCTGCTCGACCAGGCGCGTc                >  1:1636983/1‑36 (MQ=255)
                     ggtTGACGTTTCTGACCTGCTCGACCAGGCGCGTCa               >  1:5546710/1‑36 (MQ=255)
                     ggtTGACGTTTCTGACCTGCTCGACCAGGCGCGTCa               >  1:2255932/1‑36 (MQ=255)
                     ggtTGACGTTTCTGACCTGCTCGACCAGGCGCGTCa               <  1:6177420/36‑1 (MQ=255)
                     ggtTGACGTTTCTGACCTGCTCGACCAGGCGCGTCa               >  1:5975698/1‑36 (MQ=255)
                     ggtTGACGGTTCTGACCTGCTCGTCCAGGCGCGTCa               <  1:7403988/36‑1 (MQ=38)
                       ttGACGTTTCTGACCTGCTCGACCAGGCGCGTCAGc             >  1:7385455/1‑36 (MQ=255)
                       ttGACGTTTCTGACCTGCTCGACCAGGCGCGTCAGc             <  1:7052521/36‑1 (MQ=255)
                         ggcGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGt           <  1:2461343/34‑1 (MQ=255)
                         gACGTTTCTGACCTGCTCGACCAGGCGCTTCAGCGt           <  1:1042669/36‑1 (MQ=255)
                         gACGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGt           <  1:5776885/36‑1 (MQ=255)
                         gACGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGt           <  1:581231/36‑1 (MQ=255)
                         gACGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGt           <  1:3907685/36‑1 (MQ=255)
                         gACGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGt           <  1:1659633/36‑1 (MQ=255)
                         gACGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGt           >  1:2087210/1‑36 (MQ=255)
                         aaCGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGt           <  1:1050079/35‑1 (MQ=255)
                          aCGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGTg          >  1:5724270/1‑36 (MQ=255)
                          aCGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGTg          <  1:201497/36‑1 (MQ=255)
                          aCGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGTg          >  1:5618469/1‑36 (MQ=255)
                           cGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGTgg         <  1:4354926/36‑1 (MQ=255)
                           cGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGTgg         <  1:4350553/36‑1 (MQ=255)
                            gTTTCTGACCTGCTCGACCAGGCGCGTCAGCGTGGc        >  1:1291346/1‑36 (MQ=255)
                             tttCTGACCTGCTCGACCAGGCGCGTCAGCGTGGCg       <  1:6344860/36‑1 (MQ=255)
                              ttCTGACCTGCTCGACCAGGCGCGTCAGCGTGGCGa      >  1:346225/1‑36 (MQ=255)
                              ttCTGACCTGCTCGACCAGGCGCGTCAGCGTGGCGa      >  1:2565946/1‑36 (MQ=255)
                               tCTGACCTGCTCGACCAGGCGCGTCAGCGTGGCGAt     <  1:4044485/36‑1 (MQ=255)
                                cTGACCTGCTCGACCAGGCGCGTCAGCGTGGCGAtt    >  1:1721131/1‑36 (MQ=255)
                                 tGACCTGCTCGACCAGGCGCGTCAGCGTGGCGAttt   <  1:2638413/36‑1 (MQ=255)
                                 tGACCTGCTCGACCAGGCGCGTCAGCGTGGCGAttt   <  1:1205839/36‑1 (MQ=255)
                                 tGACCTGCTCGACCAGGCGCGTCAGCGTGGCGAttt   >  1:5834788/1‑36 (MQ=255)
                                  gACCTGCTCGGCCAGGCGCGTCAGCGTGGCGATTTc  <  1:5320052/36‑1 (MQ=255)
                                  gACCTGCTCGACCAGGCGCGTCAGCGTGGCGCTTTc  >  1:2414827/1‑36 (MQ=255)
                                  gACCTGCTCGACCAGGCGCGTCAGCGTGGCGATTTc  <  1:6552515/36‑1 (MQ=255)
                                  gACCTGCTCGACCAGGCGAGTCAGCGTGGCGATTTc  >  1:108733/1‑36 (MQ=255)
                                  |                                   
CGACATCCTGACTTCTATGGTGGTTGACGTTTCTGACCTGCTCGACCAGGCGCGTCAGCGTGGCGATTTC  >  REL606/4383426‑4383495

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG