Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,806,8100A→G 7.4% 0.5 ‑7.5 0.087290uhpTsugar phosphate antiporter
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (28/50):  total (36/51)
Fisher's exact test strand distribution p-value = 5.70e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TACAAATAGCAGAGAGGATCAGCATGAACGGCAGGAATTGTTTGGTGTTTTTGCCGTCGGCGTAGTAGGAA  >  REL606/3806775‑3806845
                                   |                                   
taaaaaAAGCAGAGAGGATCAGCATGACCGGCAGGa                                     <  1:5697566/33‑1 (MQ=25)
tACAAATAGCAGAGAGGATCAGCATGAACGGCAGGa                                     >  1:1018306/1‑36 (MQ=255)
 aCAAATAGCAGAGAGGATCAGCATGAACGGCAGGaa                                    >  1:1839792/1‑36 (MQ=255)
  cAAATAGCAGAGAGGATCAGCATGAACGGCAGGAAt                                   <  1:6050263/36‑1 (MQ=255)
   aaaTAGCAGAGAGGATCAGCATGAACGGCAGGGAtt                                  >  1:5656411/1‑36 (MQ=255)
   aaaTAGCAGAGAGGATCAGCATGAACGGCAGGGAtt                                  >  1:6532676/1‑36 (MQ=255)
   aaaTAGCAGAGAGGATCAGCATGAACGGCAGGAAtt                                  >  1:6916200/1‑36 (MQ=255)
   aaaTAGCAGAGAGGATCAGCATGAACGGCAGGAAtt                                  <  1:6531384/36‑1 (MQ=255)
   aaaTAGCAGAGAGGATCAGCATGAACGGCAGGAAtt                                  <  1:6123806/36‑1 (MQ=255)
   aaaTAGCAGAGAGGATCAGAATGAACGGCAGGAAtt                                  <  1:5333255/36‑1 (MQ=255)
      tAGCAGAGAGGATCAGCATGAACGGCAGGGATTTtt                               >  1:3918428/1‑36 (MQ=37)
      tAGCAGAGAGGATCAGCATGAACGGCAGGAATTGtt                               <  1:3331613/36‑1 (MQ=255)
       aGCAGAGAGGATCAGCATGAACGGCAGGCATTGttt                              <  1:3577408/36‑1 (MQ=255)
       aGCAGAGAGGATCAGCATGAACGGCAGGAATTGttt                              <  1:6216709/36‑1 (MQ=255)
       aGCAGAGAGGATCAGCATGAACGGCAGGAATTGttt                              >  1:6298644/1‑36 (MQ=255)
        gCAGAGAGGATCAGCATGAACGGCAGGGATTTtttt                             >  1:3355736/1‑35 (MQ=37)
         cAGAGAGGATCAGCATGAACGGCAGGGATTGTTTgg                            >  1:1893109/1‑36 (MQ=255)
         cAGAGAGGATCAGCATGAACGGCAGGAATTGTTTgg                            >  1:1944336/1‑36 (MQ=255)
         cAGAGAGGATCAGAATGAACGGCAGGAATTGTTTgg                            <  1:5636980/36‑1 (MQ=255)
          agagagGATCAGCATGAACGGCAGGGATTGttttgg                           >  1:4376510/1‑33 (MQ=37)
          agagagGATCAGCATGAACGGCAGGAATTGTTTGgt                           <  1:4227694/36‑1 (MQ=255)
          agagagGATCAGCATGAACGGCAGGAATTGTTTGgt                           >  1:4775952/1‑36 (MQ=255)
          agagagGATCAGCATGAACGGCAGGAATTGTTTGgt                           <  1:3510515/36‑1 (MQ=255)
          agagagGATCAGCATGAACGGCAGGAATTGTTTGgt                           <  1:341657/36‑1 (MQ=255)
           gagagGATCAGCATGAACGGCAGGAATTGTTTGgtg                          <  1:4004550/36‑1 (MQ=255)
           gagagGATCAGCATGAACGGCAGGAATTGTTTGgtg                          <  1:4906075/36‑1 (MQ=255)
           gagagGATCAGCATGAACGGCAGGAATTGTTTGgta                          <  1:2899885/36‑2 (MQ=255)
            agagGATCAGCATGAACGGCAGGAATTGTTTGgtgt                         <  1:4939760/36‑1 (MQ=255)
            agagGATCAGCATGAACGGCAGGAATTGTTTGgtgt                         >  1:1037445/1‑36 (MQ=255)
            agagGATCAGCATGAACGGCAGGAATTGTTTGgtgt                         <  1:2198523/36‑1 (MQ=255)
            agagGATCAGCATGAACGGCAGGAATTGTTTGgtgt                         >  1:3957218/1‑36 (MQ=255)
             gagGATCAGCATGAACGGCAGGGATTGTTTTGTGtt                        >  1:4872317/1‑36 (MQ=37)
             gagGATCAGCATGAACGGCAGGAATTGTTTGGTGtt                        >  1:3063818/1‑36 (MQ=255)
              gggATCAGCATGAACGGGAGGAATTGTTTGGTGttt                       <  1:565603/35‑1 (MQ=255)
              agGATCAGCATGAACGGCAGGAATTGTTTGGTGttt                       <  1:4932333/36‑1 (MQ=255)
              agGATCAGCATGAACGGCAGGAATTGTTTGGTGttt                       <  1:5008412/36‑1 (MQ=255)
              agGATCAGCATGAACGGCAGGAATTGTTTGGTGttt                       >  1:1203640/1‑36 (MQ=255)
               ggATCAGCATGAACGGCAGGAATTGTTTGGTGtttt                      <  1:897220/36‑1 (MQ=255)
               ggATCAGCATGAACGGCAGGAATTGTTTGGTGtttt                      <  1:3953805/36‑1 (MQ=255)
                gATCAGCATGAACGGCAGGAATTGTTTGGTGttttt                     >  1:1712965/1‑36 (MQ=255)
                gATCAGCATGAACGGCAGGAATTGTTTGGTGttttt                     <  1:5892826/36‑1 (MQ=255)
                gATCAGCATGAACGGCAGGAATTGTTTGGTGttttt                     <  1:221715/36‑1 (MQ=255)
                 ataaGCATGAACGGCAGGAATTGTTTGGTGTTTTTg                    <  1:2799087/33‑1 (MQ=255)
                 aTCAGCATGAACGGCAGGGATTGTTTTGTGTTTTTg                    >  1:6230394/1‑36 (MQ=37)
                 aTCAGCATGAACGGCAGGAATTGTTTGGTGtttttt                    >  1:6766183/1‑35 (MQ=255)
                 aTCAGCATGAACGGCAGGAATTGTTTGGTGtttttt                    >  1:5330597/1‑35 (MQ=255)
                 aTCAGCATGAACGGCAGGAATTGTTTGGTGTTTTTg                    <  1:1516090/36‑1 (MQ=255)
                 aTCAGCATGAACGGCAGGAATTGTTTGGTGTTTTTg                    <  1:756713/36‑1 (MQ=255)
                  tCAGCATGAACGGCAGGAATTGTTTGGTGttttttc                   >  1:3438670/1‑34 (MQ=255)
                  tCAGCATGAACGGCAGGAATTGTTTGGTGttttttc                   >  1:1599704/1‑34 (MQ=255)
                  tCAGCATGAACGGCAGGAATTGTTTGGTGTTTTTGc                   >  1:6664403/1‑36 (MQ=255)
                  tCAGCATGAACGGCAGGAATTGTTTGGTGTTTTTGc                   <  1:1494639/36‑1 (MQ=255)
                  tCAGCATGAACGGCAGGAATTGTTTGGTGTTTTTGc                   <  1:4707325/36‑1 (MQ=255)
                  tCAGCATGAACGGCAGGAATTGTTTGGTGTTTTTGc                   <  1:650701/36‑1 (MQ=255)
                  tCAGCATGAACGGCAGGAATTGTTTGGTGTTTTCGc                   <  1:7478627/36‑1 (MQ=255)
                   cAGCATGAACGGCAGGAATTGTTTGGTGTTTTTTcc                  >  1:1528434/1‑36 (MQ=255)
                   cAGCATGAACGGCAGGAATTGTTTGGTGTTTTTGcc                  <  1:1641138/36‑1 (MQ=255)
                   cAGCATGAACGGCAGGAATTGTTTGGTGTTTTTGcc                  <  1:1141348/36‑1 (MQ=255)
                     gCATGAACGGCAGGAATTGTTTGGTGTTTTTGCCGt                <  1:4033915/36‑1 (MQ=255)
                     gCATGAACGGCAGGAATTGTTTGGTGTTTTTGCCGt                >  1:2726859/1‑36 (MQ=255)
                     gCATGAACGGCAGGAATTGTTTGGTGTTTTTGCCGt                <  1:6917736/36‑1 (MQ=255)
                     gCATGAACGGCAGGAATTGTTTGGTGTTTTTGCCGt                >  1:7243768/1‑36 (MQ=255)
                      cATGAACGGCAGGAATTGTTTGGTGTTTTTGCCGTc               <  1:4796228/36‑1 (MQ=255)
                       aTGAACGGCAGGAATTGTTTGGTGTTTTTTCCGTCg              >  1:4281237/1‑36 (MQ=255)
                        tGAACGGCAGGAATTGTTTGGTGTTTTTTCCGTCgg             >  1:6798707/1‑36 (MQ=255)
                        tGAACGGCAGGAATTGTTTGGTGTTTTTGCCGTCgg             >  1:202102/1‑36 (MQ=255)
                        ggAACGGCAGGAATTGTTTGGTGTTTTTGCCGTCgg             <  1:3389075/35‑1 (MQ=255)
                         gAACGGCAGGAATTGTTTGGTGTTTTTGCCGTCGGc            <  1:4221363/36‑1 (MQ=255)
                          aaCGGCAGGAATTGTTTGGTGTTTTTGCCGTCGGCg           <  1:7052155/36‑1 (MQ=255)
                           aCGGCAGGAATTGTTTGGTGTTTTTTCCGTCGGCgt          >  1:6441067/1‑36 (MQ=255)
                           aCGGCAGGAATTGTTTGGTGTTTTTGCCGTCGGCgt          <  1:5170224/36‑1 (MQ=255)
                           aCGGCAGGAATTGTTTGGTGTTATTGCCGTCGGCgt          <  1:4314462/36‑1 (MQ=255)
                            cGGCAGGAATTGTTTGGTGTTTTTGCCGTCGGCgta         <  1:7417059/36‑1 (MQ=255)
                            cGGCAGGAATTGTTTGGTGTTTTTGCCGTCGGCgta         <  1:7501898/36‑1 (MQ=255)
                            cGGCAGGAATTGTTTGGTGTTTTTGCCGTCGGCgta         >  1:882413/1‑36 (MQ=255)
                             ggCAGGAATTGTTTGGTGTTTTTGCCGTCGGCgtag        >  1:2235747/1‑36 (MQ=255)
                             ggCAGGAATTGTTCGGTGTTTTTGCCGTCGGCgtag        >  1:6718807/1‑36 (MQ=255)
                              gCAGGAATTGTTTGGTGTTTTTGCCGTCGGCgtagt       <  1:3877969/36‑1 (MQ=255)
                              gCAGGAATTGTTTGGTGTTTTTGCCGTCGGCgtagt       <  1:1406457/36‑1 (MQ=255)
                              gCAGGAATTGTTTGGTGTTTTTGCCGTCGGCgtagt       >  1:4475312/1‑36 (MQ=255)
                              gCAGGAATTGTTTGGTGTTTTTGCCGTCGGCgtagt       <  1:3494907/36‑1 (MQ=255)
                               cAGGAATTGTTTGGTGTTTTTGCCGTCGGCgtagta      <  1:3589431/36‑1 (MQ=255)
                                aGGAATTGTTTGGTGTTTTTGCCGTCGGCGtagtag     >  1:7151265/1‑36 (MQ=255)
                                aGGAATTGTTTGGTGTTTTTGCCGTCGGCGtagtag     <  1:1899332/36‑1 (MQ=255)
                                aGGAATTGTTTGGCTTTTTTGCCTTCGGCGtggtag     >  1:5333190/1‑36 (MQ=16)
                                 ggAATTGTTTGGTGTTTTTGCCGTCGGCGTAGTAgg    <  1:682540/36‑1 (MQ=255)
                                  gAATTGTTTGGTGTTTTTGCCGTCGGCGTAGTAGGa   >  1:6542938/1‑36 (MQ=255)
                                  gAATTGTTTGGTGTTTTTGCCGTCGGCGTAGTAGGa   <  1:6342706/36‑1 (MQ=255)
                                  gAATTGTTTGGTGTTTTTGCCGTCGGCGTAGTAGGa   <  1:3233800/36‑1 (MQ=255)
                                  gAATTGTTTGGTGTTTTTGCCGTCGGCGTAGTAGGa   <  1:5153618/36‑1 (MQ=255)
                                   aaTTGTTTGGTGTTTTTGCCGTCGGCGTAGTAGGaa  <  1:7356388/36‑1 (MQ=255)
                                   |                                   
TACAAATAGCAGAGAGGATCAGCATGAACGGCAGGAATTGTTTGGTGTTTTTGCCGTCGGCGTAGTAGGAA  >  REL606/3806775‑3806845

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG