Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,495,4370A→C 4.1% 0.2 0.0 ‑0.072395eutQhypothetical protein
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (2/1):  ref base (34/35):  total (36/36)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.07e-01

CACGCAATGCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTGCTACCATCAATCACTTTGATGCCG  >  REL606/2495403‑2495472
                                  |                                   
cACGCAATGCGGCTCTGCGCCGTCGAAACGACCaaa                                    <  1:802054/36‑1 (MQ=255)
cACGCAATGCGGCTCTGCGCCGTCGAAACGACCaaa                                    >  1:2842135/1‑36 (MQ=255)
 aCGCAATGCGGCTCTGCGCCGTCGAAACGACCAAAc                                   >  1:7872959/1‑36 (MQ=255)
 aCGCAATGCGGCTCTGCGCCGTCGAAACGACCAAAc                                   >  1:1521717/1‑36 (MQ=255)
 aCGCAATGCGGCTCTGCGCCGTCGAAACGACCAAAc                                   >  1:2623409/1‑36 (MQ=255)
  cGCAATGCGGCTCTGCGCCGTCGAAACGACCAAACt                                  >  1:4376455/1‑36 (MQ=255)
  cGCAATGCGGCTCTGCGCCGTCGAAACGACCAAACt                                  <  1:2476770/36‑1 (MQ=255)
   gCAATGCGGCTCTGCGCCGTCGAAACGACCAAACtt                                 >  1:2095674/1‑36 (MQ=255)
    cAATGCGGCTCTGCGCCGTCGAAACGACCAAACTTg                                <  1:6273128/36‑1 (MQ=255)
     actGCGGCTCTGCGCCGTCGAAACGACCAAACTTGa                               <  1:3073647/34‑1 (MQ=255)
     aaTGCGGCTCTGCGCCGTCGAAACGACCACACTTGa                               >  1:4580177/1‑36 (MQ=255)
     aaTGCGGCTCTGCGCCGTCGAAACGACCAAACTTGa                               >  1:155486/1‑36 (MQ=255)
      aTGCGGCTCTGCGCCGTCGAAACGACCAAACTTGAc                              <  1:404746/36‑1 (MQ=255)
      aTGCGGCTCTGCGCCGTCGAAACGACCAAACTTGAc                              <  1:5596869/36‑1 (MQ=255)
       tGCGGCTCTGCGGCGTCGAAACGACCAAACTTGACg                             <  1:880080/36‑1 (MQ=255)
       tGCGGCTCTGCGCCGTCGAAACGACCAAACTTGACg                             <  1:2125600/36‑1 (MQ=255)
        gCGGCTCTGCGCCGTCGAAACGACCAAACTTGACgc                            <  1:2252454/36‑1 (MQ=255)
         cGGCTCTGCGCCGTCGAAACGACCAAACTTGACgct                           >  1:68779/1‑36 (MQ=255)
         cGGCTCTGCGCCGTCGAAACGACCAAACTTGACgct                           <  1:7363326/36‑1 (MQ=255)
          ggCTCTGCGCCGTCGAAACGACCAAACTTGACgctg                          >  1:2731887/1‑36 (MQ=255)
          ggCTCTGCGCCGTCGAAACGACCAAACTTGACgctg                          >  1:406293/1‑36 (MQ=255)
          ggCTCTGCGCCGTCGAAACGAACAAACTTGACgctg                          >  1:4774059/1‑36 (MQ=255)
           gCTCTGCGCCGTCGAAACTACCAAACTTGACgctgc                         <  1:6184791/36‑1 (MQ=255)
           gCTCTGCGCCGTCGAAACGACCAAACTTGACgctgc                         >  1:2882051/1‑36 (MQ=255)
            ctctGCGCCGTCGAAACGACCAAACTTGACgctgct                        >  1:6608160/1‑36 (MQ=255)
             tctGCGCCGTCGAAACGACCAAACTTGACGCTGCTa                       <  1:5870432/36‑1 (MQ=255)
             tctGCGCCGTCGAAACGACCAAACTTGACGCTGCTa                       <  1:3285697/36‑1 (MQ=255)
             tctGCGCCGTCGAAACGACCAAACTTGACGCTGCTa                       >  1:4033334/1‑36 (MQ=255)
              ctGCGCCGTCGAAACGACCAAACTTGACGCTGCTAc                      <  1:2258799/36‑1 (MQ=255)
               tGCGCCGTCGAAACGACCAAACTTGACGCTGCTAcc                     <  1:3128646/36‑1 (MQ=255)
               tGCGCCGTCGAAACGACCAAACTTGACGCTGCTAcc                     >  1:7754064/1‑36 (MQ=255)
                gcgcCGTCGAAACGACCAAACTTGACGCTGCTACCa                    <  1:7681902/36‑1 (MQ=255)
                 cgcCGTCGAAACGACCAAACTTGACGCTGCTACCat                   <  1:4578865/36‑1 (MQ=255)
                 cgcCGTCGAAACGACCAAACTTGACGCTGCTACCat                   <  1:3999579/36‑1 (MQ=255)
                 cgcCGTCGAAACGACCAAACTTGACGCTGCTACCat                   <  1:1051143/36‑1 (MQ=255)
                  gcCGTCGAAACGACCAAACTTGACGCTGCTACCatc                  <  1:991576/36‑1 (MQ=255)
                  gcCGTCGAAACGACCAAACTTGACGCTGCTACCatc                  <  1:1071321/36‑1 (MQ=255)
                   ccGTCGAAACGACCAAACTTGACGCTGCTACCatca                 <  1:1091045/36‑1 (MQ=255)
                   ccGTCGAAACGACCAAACTTGACGCTGCTACCatca                 <  1:7838922/36‑1 (MQ=255)
                    cGTCGAAACGACCAAACTTGACGCTGCTACCatcaa                >  1:7840395/1‑36 (MQ=255)
                    cGTCGAAACGACAAAACTTGACGCTGCTACCatcaa                >  1:3323408/1‑36 (MQ=255)
                     gTCGAAACGACCACACTTGACGCTGCTACCatcaat               >  1:2484984/1‑36 (MQ=255)
                     gTCGAAACGACCAAACTTGACGCTGCTACCatcaat               >  1:5097817/1‑36 (MQ=255)
                     gTCGAAACGACCAAACTTGACGCTGCTACCatcaat               >  1:317192/1‑36 (MQ=255)
                     gTCGAAACGACCAAACTTGACGCTGCTACCatcaat               <  1:4087388/36‑1 (MQ=255)
                     gTCGAAACGACCAAACTTGACGCGGCTACCatcaag               >  1:3521081/1‑35 (MQ=255)
                      tCGAAACGACCAAACTTGACGCTGCTACCatcaatc              <  1:1840985/36‑1 (MQ=255)
                      tCGAAACGACCAAACTTGACGCTGCTACCatcaatc              >  1:494925/1‑36 (MQ=255)
                      tCGAAACGACCAAACTTGACGCGGCGACCatcaatc              >  1:1562367/1‑36 (MQ=38)
                       cGAAACGACCACACTTGACGCTGCTACCatcaatca             <  1:5527728/36‑1 (MQ=255)
                        gAAACGACCAAACTTGACGCTGCTACCATCAATCAc            <  1:1940237/36‑1 (MQ=255)
                        gAAACGACCAAACTTGACGCTGCTACCATCAATCAc            <  1:7202880/36‑1 (MQ=255)
                         aaaCGACCAAACTTGACGCTGCTACCATCAATCACt           <  1:7798915/36‑1 (MQ=255)
                         aaaCGACCAAACTTGACGCTGCTACCATCAATCACt           >  1:1828336/1‑36 (MQ=255)
                         aaaCGACCAAACTTGACGCTGCTACCATCAATCACt           >  1:2617934/1‑36 (MQ=255)
                         aaaCGACCAAACTTGACGCTGCTACCATCAATCACt           >  1:1595607/1‑36 (MQ=255)
                         aaaCGACCAAACTTGACGCTGCTACCATCAATCACt           >  1:2935583/1‑36 (MQ=255)
                          aaCGACCAAACTTGACGCTGCTACCATCAATCACtt          <  1:6389307/36‑1 (MQ=255)
                          aaCGACCAAACTTGACGCTGCTACCATCAATCACtt          >  1:91010/1‑36 (MQ=255)
                           aCGACCAAACTTGACGCTGCTACCATGAATCACttt         >  1:367082/1‑36 (MQ=255)
                           aCGACCAAACTTGACGCTGCTACCATCAATCACttt         >  1:7714040/1‑36 (MQ=255)
                            cGACCAAACTTGACGCTGCTACCATCAATCACTTTg        >  1:2129430/1‑36 (MQ=255)
                            cGACCAAACTTGACGCTGCTACCATCAATCACTTTg        >  1:4897888/1‑36 (MQ=255)
                            cGACCAAACTTGACGCTGCTACCATCAATCACTTTg        <  1:7853005/36‑1 (MQ=255)
                             gACCAAACTTGACGCTGCTACCATCAATCACTTTGa       >  1:991664/1‑36 (MQ=255)
                             gACCAAACTTGACGCTGCTACCATCAATCACTTTGa       >  1:6742577/1‑36 (MQ=255)
                             caCCAAACTTGACGCTGCTACCGTCAATCACTTTGa       <  1:7773191/35‑1 (MQ=255)
                              aCCAAACTTGACGCTGCTACCATCAATCACTTTGAt      <  1:3657861/36‑1 (MQ=255)
                              aCCAAACTTGACGCTGCTACCATCAATCACTTTGAt      >  1:3337093/1‑36 (MQ=255)
                              aCCAAACTTGACGCTGCTACCATCAATCACTTTGAt      <  1:7482739/36‑1 (MQ=255)
                              aCCAAACTTGACGCTGCTACCATCAATCACTTTGAt      >  1:586649/1‑36 (MQ=255)
                               ccAAACTTGACGCTGCTACCATCAATCACTTTGATg     <  1:7475505/36‑1 (MQ=255)
                                cAAACTTGACGCTGCTACCATCAATCACTTTGATgc    <  1:3441133/36‑1 (MQ=255)
                                cAAACTTGACGCTGCTACCATCAATCACTTTGATgc    <  1:2082561/36‑1 (MQ=255)
                                 aaaCTTGACGCTGCTACCATCAATCATTTTGATgcc   >  1:4988570/1‑36 (MQ=255)
                                 aaaCTTGACGCTGCTACCATCAATCACTTTGATgcc   >  1:152010/1‑36 (MQ=255)
                                 aaaCTTGACGCTGCTACCATCAATCACTTTGATgcc   <  1:7448691/36‑1 (MQ=255)
                                 aaaCTTGACGCTGCTACCATCAATCACTTTGATgcc   >  1:7031126/1‑36 (MQ=255)
                                  aaCTTGACGCTGCTACCATCAATCACTTTGATgccg  <  1:1499489/36‑1 (MQ=255)
                                  aaCTTGACGCTGCTACCATCAATCACTTTGATgccg  <  1:280534/36‑1 (MQ=255)
                                  aaCTTGACGCTGCTACCATCAATCACTTTGATgccg  >  1:2940403/1‑36 (MQ=255)
                                  aaCTTGACGCTGCTACCATCAATCACTTTGATgccg  <  1:3888125/36‑1 (MQ=255)
                                  aaCTTGACGCTGCTACCATCAATCACTTTGATgccg  >  1:6799054/1‑36 (MQ=255)
                                  aaCTTGACGCTGCTACCATCAATCACTTTGATgccg  <  1:342015/36‑1 (MQ=255)
                                  |                                   
CACGCAATGCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTGCTACCATCAATCACTTTGATGCCG  >  REL606/2495403‑2495472

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG