Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,190,6530A→G 7.8% 1.4 ‑9.7 0.076328yhaOpredicted transporter
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (14/54):  total (22/54)
Fisher's exact test strand distribution p-value = 5.86e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CAGGTAGGAAGCACTATCGTTGGTGATGGCGGTGGAATAGACGAACATCCAAATCACCAGCATTACGAAAT  >  REL606/3190618‑3190688
                                   |                                   
cAGGTAGGAAGCTCTATCGTTGGTGATGGCGGTGGa                                     <  1:1313194/36‑1 (MQ=255)
cAGGTAGGAAGCACTATCGTTGGTGATGGCGGTGGa                                     <  1:6257530/36‑1 (MQ=255)
cAGGTAGGAAGCACTATCGTTGGTGATGGCGGTGGa                                     <  1:5147471/36‑1 (MQ=255)
 aGGGAGGAAGCACTAGCGTTGGTGATGGCGGTGGaa                                    <  1:2672542/36‑1 (MQ=37)
  ggTAGGAAGCACTATCGTTGGTGATGGCGGTGGAAt                                   <  1:3075521/36‑1 (MQ=255)
  ggTAGGAAGCACTATCGTTGGTGATGGCGGTGGAAt                                   <  1:7372994/36‑1 (MQ=255)
  ggTAGGAAGAACTATCGTTGGTGATGGAGGGGGGAt                                   >  1:5112431/1‑36 (MQ=16)
   gTAGGAAGCACTATCGTTGGTGATGGCGGTGGAATa                                  <  1:4895267/36‑1 (MQ=255)
   gTAGGAAGCACTATCGTTGGTGATGGCGGTGGAATa                                  <  1:6877208/36‑1 (MQ=255)
   gTAGGAAGCACTATCGTTGGTGATGGCGGAGGAAtt                                  >  1:6715876/1‑35 (MQ=255)
     aGGAAGCACTATCCTTGGTGATGGCGGTGGAATAGa                                >  1:3199020/1‑36 (MQ=255)
      ggAAGCACTATCGTTGGTGATGGCGGTGGAATAGAc                               <  1:6940257/36‑1 (MQ=255)
       gAAGCACTATCGTTGGTGATGGCGGTGGGATAGACg                              >  1:6222748/1‑36 (MQ=255)
       gAAGCACTATCGTTGGTGATGGCGGTGGGATAGACg                              >  1:5586752/1‑36 (MQ=255)
       gAAGCACTATCGTTGGTGATGGCGGTGGAATAGACg                              <  1:6431980/36‑1 (MQ=255)
       gAAGCACTATCGTTGGTGATGGCGGTGGAATAGACg                              <  1:2116897/36‑1 (MQ=255)
        aaGCACTATCGTTGGTGATGGCGGTGGTATTGAAGa                             >  1:1254519/1‑36 (MQ=25)
         atcACTATCGTTGGTGATGGCGGTGGAATAGACGaa                            <  1:7485381/34‑1 (MQ=255)
          gCACTATCGTTGGTGATGGCGGTGGAATAGACGaaa                           >  1:2441819/1‑35 (MQ=255)
          gCACTATCGTTGGTGATGGCGGTGGAATAGACGAAc                           <  1:5193017/36‑1 (MQ=255)
          gCACTATCGTTGGTGATGGCGGTGGAATAGACGAAc                           <  1:4839621/36‑1 (MQ=255)
          gCACTATCGTTGGTGATGGCGGTGGAATAGACGAAc                           <  1:3993178/36‑1 (MQ=255)
            aCTATCGTTGTTGATGGCGGTGGAATAGACGAACAt                         <  1:5181564/36‑1 (MQ=255)
            aCTATCGTTGGTGATGGCGGTGGGATAGACGAAAAt                         >  1:3119522/1‑36 (MQ=37)
            aCTATCGTTGGTGATGGCGGTGGGATAGAAGAACAt                         >  1:4118514/1‑36 (MQ=37)
            aCTATCGTTGGTGATGGCGGTGGAATAGACGAACAt                         <  1:6457155/36‑1 (MQ=255)
            aCTATCGTTGGTGATGGCGGTGGAATAGACGAACAt                         >  1:6847649/1‑36 (MQ=255)
            aCTATCGTTGGCGATGGCGGTGGAATAGACGAACAt                         <  1:3900020/36‑1 (MQ=255)
            aCTATCGGTGGTGATGGCGGTGGAATAGACGAACAt                         <  1:6537541/36‑1 (MQ=255)
             cTATCGTTGGTGATGGCGGTGGAATAGACGAACATc                        <  1:5205264/36‑1 (MQ=255)
             cTATCGTTGGTGATGGCGGTGGAATAGACGAACATc                        <  1:4459255/36‑1 (MQ=255)
              tATCGTTGGTGATTTCGGTTGAATAGACGAACATcc                       >  1:5448095/1‑36 (MQ=37)
              tATCGTTGGTGATGGCGGTGGAATAGACGAAAATcc                       >  1:3535984/1‑36 (MQ=255)
               aTCGTTGGTGATGGCGGTGGGATAGACGAACAGCCa                      >  1:1623040/1‑36 (MQ=37)
               aTCGTTGGTGATGGCGGTGGAATAGACGAACATCCa                      <  1:6361231/36‑1 (MQ=255)
                tCGTTGGTGATGGCGGTGGAATAGACGAACATCCaa                     <  1:687192/36‑1 (MQ=255)
                 cGTTGGTGATGGCGGTGGAATAGACGAACATCCaaa                    <  1:1440858/36‑1 (MQ=255)
                 cGTTGGTGATGGCGGTGGAATAGACGAACATCCaaa                    <  1:4797632/36‑1 (MQ=255)
                  gTTGGTGATGGCGGTGGAATAGACGAACATCCAAAt                   <  1:5020584/36‑1 (MQ=255)
                  gTTGGTGATGGCGGTGGAATAGACGAACATCCAAAt                   <  1:284019/36‑1 (MQ=255)
                  gTTGGTGATGGCGGTGGAATAGACGAACATCCAAAt                   <  1:6796462/36‑1 (MQ=255)
                  gTTGGTGATGGCGGTGGAATAGACGAACATCCAAAt                   <  1:3722875/36‑1 (MQ=255)
                  gTTGGTGATGGCGGTGGAATAGACGAAAATTCAAAt                   >  1:1874030/1‑36 (MQ=37)
                   ttGGTGATGGCGGTGGGAAAGACGAAAATACAAATc                  >  1:4091410/1‑36 (MQ=16)
                   ttGGTGATGGCGGTGGAATAGACGAACATCCAAATc                  <  1:6154331/36‑1 (MQ=255)
                   ttGGTGATGGCGGTGGAATAGACGAACATCCAAATc                  <  1:3833149/36‑1 (MQ=255)
                   ttGGTGATGGCGGTGGAATAGACGAACATCCAAATc                  <  1:89282/36‑1 (MQ=255)
                   ttGGTGATGGCGGTGGAATAGACGAACATCCAAATc                  <  1:1417217/36‑1 (MQ=255)
                    tGGTGATGGCGGTGGAATAGACGAACATCCAAATCa                 <  1:5298733/36‑1 (MQ=255)
                     ggTGATGGCGGTGGAATAGACGAACATCCAAATCAc                <  1:5357117/36‑1 (MQ=255)
                      gTGATGGCGGTGGAATAGACGAACATCCAAATCAcc               <  1:3148979/36‑1 (MQ=255)
                      gTGATGGCGGTGGAATAGACGAACATCCAAATCAcc               <  1:1649122/36‑1 (MQ=255)
                       tGATGGCGGTGGAATAGACGAACATCCAAATCACCa              <  1:5137039/36‑1 (MQ=255)
                       tGATGGCGGTGGAATAGACGAACATCCAAATCACCa              <  1:2758082/36‑1 (MQ=255)
                       tGATGGCGGTGGAATAGACGAACATCCAAATCACCa              <  1:1448181/36‑1 (MQ=255)
                        gATGGCGGTGGAATAGACGAACATTCAAATCACCAg             >  1:7180583/1‑36 (MQ=255)
                        gATGGCGGTGGAATAGACGAACATCCAAATCACCAg             <  1:2328479/36‑1 (MQ=255)
                        gATGGCGGTGGAATAGACGAACATCCAAATCACCAg             <  1:603608/36‑1 (MQ=255)
                        gATGGCGGTGGAATAGACGAACATCCAAATCACCAg             <  1:948650/36‑1 (MQ=255)
                         aTGGCGGTGGAATAGACGAACATCCAAATCACCAGc            <  1:3363803/36‑1 (MQ=255)
                         aTGGCGGTGGAATAGACGAACATCCAAATAACCAGc            >  1:1192/1‑36 (MQ=255)
                          tGGCGGTGGAATAGACGAACATCCAAATCACCAGCa           <  1:7411463/36‑1 (MQ=255)
                           ggaggTGGAATAGACGAACATCCAAATCACCAGCAt          <  1:503745/33‑1 (MQ=255)
                           ggCGGTGGAATAGACGAACATCCAAATCACCCGCAt          >  1:3014371/1‑36 (MQ=255)
                           ggCGGTGGAATAGACGAACATCCAAATCACCAGCAt          >  1:989275/1‑36 (MQ=255)
                           ggCGGTGGAATAGACGAACATCCAAATCACCAGCAt          <  1:5332000/36‑1 (MQ=255)
                           ggCGGTGGAATAGACGAACATCCAAATCACCAGCAt          >  1:4878674/1‑36 (MQ=255)
                           ggCGGTGGAATAGACGAACATCCAAATCACCAGCAt          <  1:340471/36‑1 (MQ=255)
                            gCGGTGGAATAGACGAACATCCAAATCACCAGCAtt         <  1:3510983/36‑1 (MQ=255)
                             cGGTGGAATAGACGAACATCCAAATCACCAGCATTa        <  1:4144455/36‑1 (MQ=255)
                              ggTGGAATAGACGAACATCCAAATCACCAGCATTAc       <  1:637684/36‑1 (MQ=255)
                              ggTGGAATAGACGAACATCCAAATCACCAGCATTAc       <  1:639217/36‑1 (MQ=255)
                              ggTGGAATAGACGAACATCCAAATCACCAGCATTAc       <  1:5035597/36‑1 (MQ=255)
                              ggTGGAATAGACGAACATCCAAATCACCAGCATTAc       <  1:15163/36‑1 (MQ=255)
                               gTGGAATAGACGAACATCCAAATCACCAGCATTACg      >  1:4795968/1‑36 (MQ=255)
                                tGGAATAGACGAACATCCAAATCACCAGCATTACGa     <  1:6030197/36‑1 (MQ=255)
                                 ggAATAGACGAACATCCAAATCACCAGCATTACGaa    <  1:2088667/36‑1 (MQ=255)
                                 ggAATAGACGAACATCCAAATCACCAGCATTACGaa    <  1:4244540/36‑1 (MQ=255)
                                  gAATAGACGAACATCCAAATCACCAGCATTGCGaaa   <  1:381783/36‑1 (MQ=255)
                                  gAATAGACGAACATCCAAATCACCAGCATTACGaaa   >  1:3145266/1‑36 (MQ=255)
                                   aaTAGACGAACATCCAAATCACCAGCATTACGAAAt  <  1:3884681/36‑1 (MQ=255)
                                   |                                   
CAGGTAGGAAGCACTATCGTTGGTGATGGCGGTGGAATAGACGAACATCCAAATCACCAGCATTACGAAAT  >  REL606/3190618‑3190688

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG