Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,190,6530A→G17.0% 7.9 ‑18.6 0.068328yhaOpredicted transporter
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (13/0):  ref base (8/46):  total (22/46)
Fisher's exact test strand distribution p-value = 8.02e-09
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CAGGTAGGAAGCACTATCGTTGGTGATGGCGGTGGAATAGACG.AACATCCAAATCACCAGCATTACGAA  >  REL606/3190618‑3190686
                                   |                                  
cAGGTAGTAAGCGCTAGCGTTTGTGATGGCGGTGGa                                    <  1:4846399/36‑1 (MQ=16)
cAGGTAGGAAGCACTATCGTTGGTGATGGCGGTGGa                                    <  1:5462381/36‑1 (MQ=255)
 aGGTAGGAAGCACTATCGTTGGTGATGGCGGTGGaa                                   <  1:361399/36‑1 (MQ=255)
  ggTAGGAAGCACTATCGTTGGTGATGGCGGTGGAAt                                  <  1:3442150/36‑1 (MQ=255)
  ggTAGGAAGCACTATCGTTGGTGATGGCGGGGGGAt                                  >  1:2628469/1‑36 (MQ=37)
   gTAGGAAGCCCTATCGTTGGTGATGGCGGTGGAATa                                 <  1:2155792/36‑1 (MQ=255)
   gTAGGAAGCACTATCGTTGGGGATGGAGGTGGGAtt                                 >  1:4867720/1‑35 (MQ=21)
    tAGGAAGCACTATCGTTGGTGATGGCGGTGGGATat                                >  1:2806677/1‑35 (MQ=255)
    tAGGAAGCACTATCGTTGGTGATGGAGGTGGGATAg                                >  1:4742681/1‑36 (MQ=25)
     aGGAAGCACTATCGTTGGTGATGGCGGTGGGAAAGa                               >  1:5174985/1‑36 (MQ=37)
     aGGAAGCACTATCGTTGGTGATGGCGGTGGAATAGa                               <  1:351911/36‑1 (MQ=255)
      ggAAGCACTATCGTTGGTGATGGCGGTGGAATAGAc                              <  1:3795236/36‑1 (MQ=255)
      ggAAGCACTATCGTTGGTGATGGCGGTGGAATAGAc                              <  1:4699123/36‑1 (MQ=255)
      ggAAGCACTATCGTTGGTGATGGCGGGGGGATATAc                              >  1:2176753/1‑36 (MQ=25)
      ggAAGCACTATCGTTGGTGATGGCGGCGGGATATAc                              >  1:853211/1‑36 (MQ=25)
       gAAGCACTATCGTTGGTGATGGCGGTGGAATAGACg                             <  1:1306161/36‑1 (MQ=255)
        aaGCACTATCGTTGGTGATGGCGGTGGAATAGACG.a                           <  1:5142038/36‑1 (MQ=255)
         aGCACTATCGTTGGTGATGGCGGTGGAGTAGCCgga                           >  1:3048375/1‑34 (MQ=25)
         aGCACTATCGTTGGTGATGGCGGTGGAATAGACG.aa                          <  1:1959063/36‑1 (MQ=255)
         aGCACTATCGTTGGTGATGGCGGTGGAATAGACG.aa                          <  1:2861250/36‑1 (MQ=255)
         aGCACTATCGTTGGTGATGGCGGTGGAATAGACG.aa                          <  1:2727598/36‑1 (MQ=255)
          gCACTATCGTTGGTGATGGCGGTGGGATAGACG.AAc                         >  1:5218070/1‑36 (MQ=255)
           cACTATCGTTGATGCTGGCGGTGGAGTCGACG.CACa                        <  1:925770/36‑1 (MQ=12)
            aCTATCGTTGGTGATGGCGGTGGAATAGACG.AACAt                       <  1:1848881/36‑1 (MQ=255)
            aCTATCGTTGGTGATGGCGGTGGAATAGACG.AACAt                       <  1:5125591/36‑1 (MQ=255)
             cTATCGTTGGTGATGGCGGTGGAATAGACG.AACATc                      <  1:2826460/36‑1 (MQ=255)
              tATCGTTGGTGATGGCGGTGGAATAGACG.AACATcc                     <  1:1829696/36‑1 (MQ=255)
              tATCGTTGGTGATGGCGGTGGAATAGACG.AACATcc                     <  1:4233087/36‑1 (MQ=255)
               aTCGTTGGTGATGGCGGTGGGATTGAGG.GACATCCa                    >  1:5705312/1‑36 (MQ=16)
                tCGTTTTTGATGGCGGTTGACTATACG.AACATCCaa                   <  1:4612134/36‑1 (MQ=16)
                tCGTTGGTGATGGCGGTGGAATAGACG.AACATCCaa                   <  1:4736975/36‑1 (MQ=255)
                tCGTTGGTGATGGCGGTGGAATAGACG.AACATCCaa                   <  1:1422840/36‑1 (MQ=255)
                tCGTTGGTGATGGCGGTGGAATAGACG.AACATCCaa                   <  1:713727/36‑1 (MQ=255)
                 cGTTGGTGATGGCGGTGGAATAGACG.AACATCCaaa                  <  1:5936325/36‑1 (MQ=255)
                 cGTTGGTGATGGCGGTGGAATAGACG.AACATCCaaa                  <  1:5133978/36‑1 (MQ=255)
                  gTTGGTGATGGCGGTGGGATATAAG.AACATCCaaaa                 >  1:3973772/1‑35 (MQ=21)
                  gTTGGTGATGGCGGTGGGAGTGACGAAAAATCCaaa                  >  1:2424280/1‑36 (MQ=9)
                  gTTGGTGATGGCGGTGGAATAGACG.AACATCCAAAt                 <  1:5049899/36‑1 (MQ=255)
                   ttGGTTATGGCGGTGGAATAGACG.AACATCCGAATc                <  1:4584779/36‑1 (MQ=38)
                   ttGGTGATGGCGGTGGAATAGACG.AACATCCAAATc                <  1:4165756/36‑1 (MQ=255)
                   ttGGTGATGGCGGTGGAATAGACG.AACATCCAAATc                <  1:110193/36‑1 (MQ=255)
                   ttGGTGATGGCGGTGGAATAGACG.AACATCCAAATc                >  1:5614181/1‑36 (MQ=255)
                   ttGGTGATGGCGGTGGAATAGAAG.AACATCCAAAtt                >  1:1710296/1‑35 (MQ=255)
                    tGGTGATGGCGTTGGAATAGACG.AACATCCAAATCa               <  1:2732626/36‑1 (MQ=255)
                    tGGTGATGGCGGTGGAATAGACG.AACATCCAAATCa               <  1:1046117/36‑1 (MQ=255)
                     ggTGATGGCGGTGGTATATACG.AAAAACCAAATCat              >  1:5644230/1‑35 (MQ=19)
                     ggTGATGGCGGTGGAATAGACG.AACATCCAAATCAc              <  1:2821548/36‑1 (MQ=255)
                      gTGATGGCGGTGGAATAGAAG.AAAATCCAAAACAcc             >  1:953590/1‑36 (MQ=25)
                        gATGGCGGTGGGATAGACG.AACAACCAAAACACCAg           >  1:1334986/1‑36 (MQ=37)
                         aTGGCGGTGGGGTAGAAG.AAAAACCAAAAAACCacc          >  1:5460589/1‑34 (MQ=0)
                         aTGGCGGTGGAATAGACG.AACCTCCAAATCACCAGc          <  1:3054137/36‑1 (MQ=255)
                          tGGCGGTGGAATATACG.AACATCCAAATCACCAGCa         <  1:2470345/36‑1 (MQ=255)
                          tGGCGGTGGAATAGACG.AACATCCAAATCACCAGCa         <  1:3662127/36‑1 (MQ=255)
                           ggCGGTGGAATAGACG.AACATCCAAATCACCAGCAt        <  1:514544/36‑1 (MQ=255)
                           ggCGGTGGAATAGACG.AACATCCAAATCACCAGCAt        <  1:5559235/36‑1 (MQ=255)
                            gCGGTGGAATAGACG.AACATCCAAATCACCAGCAtt       <  1:4418340/36‑1 (MQ=255)
                            gCGGTGGAATAGACG.AACATCCAAATCACCAGCAtt       <  1:4289910/36‑1 (MQ=255)
                            gCGGTGGAATAGACG.AACATCCAAATCACCAGCAtt       <  1:1598618/36‑1 (MQ=255)
                             cGGTGAAATAGACG.AACATCCAAATCACCAGCATTa      <  1:4548978/36‑1 (MQ=255)
                              ggTGGAATAGACG.AACATCCAAATCACCAGCATTAc     >  1:2181531/1‑36 (MQ=255)
                              ggTGGAATAGACG.AACATCCAAATCACCAGCATTAc     <  1:2565335/36‑1 (MQ=255)
                              ggTGGAATAGACG.AACATCCAAATCACCAGCATTAc     <  1:3224841/36‑1 (MQ=255)
                              ggTGGAATAGACG.AACATCCAAATCACCAGCATTAc     <  1:3607431/36‑1 (MQ=255)
                               gTGGAATAGACG.AACATCCAAATCACCAGCATTACg    <  1:587518/36‑1 (MQ=255)
                               gTGGAATAGACG.AACATCCAAATCACCAGCATTACg    <  1:4272271/36‑1 (MQ=255)
                                tGGAATAGACG.AACATCCAAATCACCAGCATTACGa   <  1:467384/36‑1 (MQ=255)
                                 ggAATAGACG.AACATCCAAATCCCCAGCATTACGaa  <  1:1778700/36‑1 (MQ=255)
                                 ggAATAGACG.AACATCCAAATCACCAGCATTACGaa  >  1:4749961/1‑36 (MQ=255)
                                 ggAATAGACG.AACATCCAAATCACCAGCATTACGaa  <  1:863515/36‑1 (MQ=255)
                                 ggAATAGACG.AACATACAAATCACCAGCATTACGaa  >  1:4403566/1‑36 (MQ=255)
                                 ggAATAGACG.AAAATCCAAATCACCAGCATTACGaa  >  1:170787/1‑36 (MQ=255)
                                   |                                  
CAGGTAGGAAGCACTATCGTTGGTGATGGCGGTGGAATAGACG.AACATCCAAATCACCAGCATTACGAA  >  REL606/3190618‑3190686

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 14 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG