Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,476,2330A→T 7.7% 0.6 ‑7.1 ‑0.067347idnOgluconate 5‑dehydrogenase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (16/44):  total (23/44)
Fisher's exact test strand distribution p-value = 8.21e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.53e-01

CCGCTTGCGATACCAGGAACACGGCGGTCTGGTTTACTGCGATCACATCATTCCACTCTTGTTCAGGGAAC  >  REL606/4476198‑4476268
                                   |                                   
ccGCTTGCGATACCAGGAACACGGCGGTCTGGTTTa                                     <  1:6174989/36‑1 (MQ=255)
 cGCTTGCGATACCAGGAACACGGCGGTCTGGTTTAc                                    <  1:4596123/36‑1 (MQ=255)
  gCTTGCGATACCAGGAACACGGCGGTCTGGTTTACt                                   >  1:5435042/1‑36 (MQ=255)
   cTTGCGATACCAGGAACACGGCGGTCTGGTTTACTg                                  <  1:3484378/36‑1 (MQ=255)
   cTTGCGATACCAGGAACACGGCGGTCTGGTTTACTg                                  <  1:7080035/36‑1 (MQ=255)
    ttGCGATACCAGGAACACGGCGGTCTGGTTTACTGc                                 <  1:2715740/36‑1 (MQ=255)
    ttGCGATACCAGGAACACGGCGGTCTGGTTTACTGc                                 <  1:1766959/36‑1 (MQ=255)
    gtGCGATACCAGGAACACGGCGGTCTGGTTTACTGc                                 <  1:3848713/35‑1 (MQ=255)
     tGCGATACCAGGAACACGGCGGTCTGGTTTACTGCg                                <  1:1337075/36‑1 (MQ=255)
      gCGATACCAGGAACACGGCGGTCTGGTTTACTGCGa                               <  1:5524820/36‑1 (MQ=255)
       cGATACCAGGAACACGGCGGTCTGGTTTACTGCGAt                              <  1:2059501/36‑1 (MQ=255)
        gATACCAGGAACACGGCGGTCTTGTTTGCTGCGATc                             >  1:995267/1‑36 (MQ=37)
         aTACCAGGAACACGGCGGTCTGGTTTACTGCGATca                            <  1:4371408/36‑1 (MQ=255)
          tACCAGGAACACGGCGGTCTGGTTTACTGCGATcac                           <  1:1216132/36‑1 (MQ=255)
           aCCAGGAACACGGCGGTCTGGTTTACTGCGATtaca                          >  1:5642074/1‑36 (MQ=255)
           aCCAGGAACACGGCGGTCTGGTTTACTGCGATcacc                          <  1:5050983/36‑2 (MQ=255)
           aCCAGGAACACGGCGGTCTGGTTTACTGCGATcaca                          <  1:557261/36‑1 (MQ=255)
           aCCAGGAACACGGCGGTCTGGGTTTCTTCGATcaca                          >  1:4860055/1‑36 (MQ=25)
            ccAGGAACACGGCGGTCTGGTTTTCTGCGATcacac                         >  1:1962430/1‑35 (MQ=255)
            ccAGGAACACGGCGGTCTGGTTTACTGCGATCAcat                         <  1:7635459/36‑1 (MQ=255)
            ccAGGAACACGGCGGTCTGGTTTACTGCGATCAcat                         <  1:6146694/36‑1 (MQ=255)
            ccAGGAACACGGCGGTCTGGTTTACTGCGATCAcat                         <  1:3748067/36‑1 (MQ=255)
             cAGGAACACGGCGGTCTGGTTTTCTTCGATCCcatc                        >  1:4887561/1‑36 (MQ=25)
             cAGGAACACGGCGGTCTGGTTTACTGCGATCAcatc                        <  1:5256539/36‑1 (MQ=255)
             cAGGAACACGGCGGTCTGGTTTACTGCGATCAcatc                        <  1:989260/36‑1 (MQ=255)
             cAGGAACACGGCGGTCTGGTTTACTGCGATCAcatc                        <  1:2174609/36‑1 (MQ=255)
              aGGAACACGGCGGTCTGGTTTTTTGCGATCAcatca                       >  1:2948396/1‑36 (MQ=37)
              aGGAACACGGCGGTCTGGTTTACTGCGATCAcatca                       <  1:5357906/36‑1 (MQ=255)
              aGGAACACGGCGGTCTGGGTTTCTTTTATCAcctca                       >  1:4499729/1‑36 (MQ=12)
              aGGAACACGGCGGGCTGGTTTACTGCGATTAcatca                       >  1:1521823/1‑36 (MQ=37)
               ggAACACGGCGGTCTGGTTTACTGCGATCAcatcat                      <  1:7262903/36‑1 (MQ=255)
               ggAACACGGCGGTCTGGTTTACTGCGATCAcatcat                      <  1:5559471/36‑1 (MQ=255)
                gAACACGGCGGTCTGGTTTACTGCGATCACATCAtt                     <  1:1981341/36‑1 (MQ=255)
                 aaCACGGCGGTCTGGTTTACTGCGATCACATCATTc                    <  1:3363531/36‑1 (MQ=255)
                 aaCACGGCGGTCTGGTTTACTGCGATCACATCATTc                    <  1:3033919/36‑1 (MQ=255)
                 aaCACGGCGGTCCGGTGTATTGCGCTCATAGCATTc                    >  1:1518988/1‑36 (MQ=12)
                  acacGGCGGTCTGGTTTACTGCGATCACATCATTcc                   <  1:5337330/36‑1 (MQ=255)
                   cacGGCGGTCTGGTTTACTGCGATCACATCATTCCa                  <  1:5216043/36‑1 (MQ=255)
                    acGGCGGTCTGGTTTACTGCGATCACATCATTCCAc                 <  1:7268518/36‑1 (MQ=255)
                     cggcggTCTGGTTTACTGCGATCACATCATTCCAct                <  1:1842368/36‑1 (MQ=255)
                      ggcggTCTGGTTTTCTGCGATCACATCATTTCActc               >  1:4148231/1‑36 (MQ=37)
                      ggcggTCTGGTTTACTGCGATCACATCATTCCActc               <  1:7657546/36‑1 (MQ=255)
                       gcggTCTGGTTTACTGCGATCACATCATTCCActct              <  1:7248561/36‑1 (MQ=255)
                       gcggTCTGGTTTACTGCGATCACATCATTCCActct              >  1:6327009/1‑36 (MQ=255)
                       gcggTCTGGTTTACTGCGATCACATCATTCCActct              >  1:5507597/1‑36 (MQ=255)
                       gcggTCTGGTTTACTGCGATCACATCATTCCActct              <  1:2541513/36‑1 (MQ=255)
                       gcggTCTGGTTTACTGCGATCACATCATTCCActct              <  1:304787/36‑1 (MQ=255)
                       gcggTCTGGTTTACTGCGATCACATCATTCCActct              <  1:4914653/36‑1 (MQ=255)
                       gcggTCTGGTTTACTGCGATCACATCATTCCActct              <  1:3666200/36‑1 (MQ=255)
                        cggTCTGGTTTACTGCGATCACATCATTCCACTCtt             >  1:6110251/1‑36 (MQ=255)
                        cggTCTGGTTTACTGCGATCACATCATTCCACTCtt             <  1:1593157/36‑1 (MQ=255)
                        cggTCTGGTTTACTGCGATCACATCATTCCACTCtt             <  1:7872659/36‑1 (MQ=255)
                         ggTCTTATTTACTTCGATCACATCATTCCACTCTTg            <  1:4051116/36‑1 (MQ=25)
                         ggTCTGGTTTACTGCGATCACATCATTCCACTCTTg            >  1:6245780/1‑36 (MQ=255)
                         ggTCTGGTTTACTGCGATCACATCATTCCACTCTTg            >  1:3934801/1‑36 (MQ=255)
                          gTCTGGTTTACTGCGATCACATCATTCCACTCTTGt           >  1:1544719/1‑36 (MQ=255)
                           tCTGGTTTACTGCGATCACATCATTCCACTCTTGtt          <  1:886486/36‑1 (MQ=255)
                            cTGGTTTACTGCGATCACATCATTCCACTCTTGTTc         <  1:237146/36‑1 (MQ=255)
                            cTGGTTTACTGCGATCACATCATTCCACTCTTGTTc         >  1:3164888/1‑36 (MQ=255)
                              ggTTTACTGCGATCACATCATTCCACTCTTGTTCAg       >  1:7425278/1‑36 (MQ=255)
                               gTTTACTGCGATCACATCATTCCACTCTTGTTCAgg      >  1:1213892/1‑36 (MQ=255)
                               gTTTACTGCGATCACATCATTCCACTCTTGTTCAgg      >  1:3437504/1‑36 (MQ=255)
                               gTTTACTGCGATCACATCATTCCACTCTTGTTCAgg      <  1:2209368/36‑1 (MQ=255)
                                tttACTGCGATCACATCATTCCACTCTTGTTCAggg     <  1:3456282/36‑1 (MQ=255)
                                tttACTGCGATCACATCATTCCACTCTTGTTCAggg     >  1:3555848/1‑36 (MQ=255)
                                tttACTGCGATCACATCATTCCACTCTTGTTCAggg     <  1:5220988/36‑1 (MQ=255)
                                 ttACTGCGATCACATCATTCCACTCTTGTTCAGGGa    >  1:7759352/1‑36 (MQ=255)
                                  tACTGCGATCACATCATTCCACTCTTGTTCAGGGaa   <  1:6772287/36‑1 (MQ=255)
                                   aCTGCGATCACATCATTCCTCTCTTGTTCAGGGAAc  <  1:1391124/36‑1 (MQ=255)
                                   aCTGCGATCACATCATTCCACTCTTGTTCAGGGAAc  >  1:3328620/1‑36 (MQ=255)
                                   |                                   
CCGCTTGCGATACCAGGAACACGGCGGTCTGGTTTACTGCGATCACATCATTCCACTCTTGTTCAGGGAAC  >  REL606/4476198‑4476268

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG