Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,315,8480C→A12.1% 0.8 ‑14.0 0.075238yhcPp‑hydroxybenzoic acid efflux system component
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/14):  ref base (43/18):  total (43/32)
Fisher's exact test strand distribution p-value = 8.41e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAA  >  REL606/3315813‑3315883
                                   |                                   
tGCGCGGATCATCGCAATGATGATCACCAGTCCGGc                                     >  1:267995/1‑36 (MQ=255)
tGCGCGGATCATCGCAATGATGATCACCAGTCCGGc                                     <  1:3341041/36‑1 (MQ=255)
 gcgcGGATCATCGCAATGATGATCACCAGTCCGGCa                                    >  1:2232706/1‑36 (MQ=255)
  cgcgGATCATCGCAATGATGATCACCAGTCCGGCaa                                   <  1:7921081/36‑1 (MQ=255)
   gcgGATCATCGCAATGATGATCACCAGTCCGGCAAt                                  <  1:6444203/36‑1 (MQ=255)
   gcgGATCATCGCAATGATGATCACCAGTCCGGCAAt                                  >  1:5891784/1‑36 (MQ=255)
    cgGATCATCGCAATGATGATCACCAGTCCGGCAATa                                 <  1:4405367/36‑1 (MQ=255)
    cgGATCATCGCAATGATGATCACCAGTCCGGCAATa                                 >  1:1493913/1‑36 (MQ=255)
    cgGATCATCGCAATGATGATCACCAGTCCGGCAATa                                 <  1:185946/36‑1 (MQ=255)
    cgGATCATCGCAATGATGATCACCAGGCCGGCAATa                                 >  1:1971838/1‑36 (MQ=255)
    aggATCATCGAAATGATGATCAAAAGTCCGGCaaga                                 >  1:2823283/2‑34 (MQ=19)
     ggATCATCGCAATGATGATCACCAGTCCGGCAATAc                                >  1:1784143/1‑36 (MQ=255)
      gATCATCGCAATGATGATCACCATTCCGGCAATACa                               <  1:4995737/36‑1 (MQ=255)
      gATCATCGCAATGATGATCACCAGTCCGGCAATACa                               >  1:1137482/1‑36 (MQ=255)
      gATCATCGCAATGATGATCACCAGTCCGGCAATACa                               >  1:5219913/1‑36 (MQ=255)
       atcatcGCAATGATGATCACCAGTCCGGCAATACAg                              <  1:5021944/36‑1 (MQ=255)
       atcatcGCAATGATGATCACCAGTCCGGCAATACAg                              <  1:1816981/36‑1 (MQ=255)
       atcatcGCAATGATGATCACCAGTCCGGCAATACAg                              >  1:7641088/1‑36 (MQ=255)
        tcatcGCAATGATGATCACCAGTCCGGCAATACAGc                             >  1:6244395/1‑36 (MQ=255)
        tcatcGCAATGATGATCACCAGTCCGGCAATACAGc                             <  1:4647818/36‑1 (MQ=255)
        tcatcGCAATGATGATCACCAGTCCGGCAATACAGc                             <  1:1447917/36‑1 (MQ=255)
        tcatcGCAATGATGATCACCAGTCCGGCAATACAGc                             <  1:7126663/36‑1 (MQ=255)
        tcatcGCAATGATGATCACCAGTCCGGCAATACAGc                             >  1:933321/1‑36 (MQ=255)
        tcatcGCAATGATGATAACCAGTTCGGCAATACAGc                             <  1:2800113/36‑1 (MQ=37)
         catcGCAATGATGATCACCAGTCCGGCAATACAGcc                            <  1:4526357/36‑1 (MQ=255)
         catcGCAATGATGATCACCAGTCCGGCAATACAGcc                            <  1:3743918/36‑1 (MQ=255)
         catcGCAATGATGATCACCAGTCCGGCAATACAGcc                            >  1:7479218/1‑36 (MQ=255)
         catcGCAATGATGATCACCAGTCCGGCAATACAGcc                            <  1:1232964/36‑1 (MQ=255)
          atcGCAATGATGATCACCAGTCCGGCAATACAGCCa                           <  1:2818950/36‑1 (MQ=255)
          atcGCAATGATGATCACCAGTCCGGCAATACAGCCa                           <  1:7922119/36‑1 (MQ=255)
          atcGCAATGATGATCACCAGTCCGGCAATACAGCCa                           >  1:235605/1‑36 (MQ=255)
           tcGCAATGATGATCACCAGTCCGGCAATACAGTCaa                          <  1:1630667/36‑1 (MQ=255)
           tcGCAATGATGATCACCAGTCCGGCAATACAGCCaa                          >  1:6988615/1‑36 (MQ=255)
           tcGCAATGATGATCACCAGTCCGGCAATACAGCCaa                          >  1:2278989/1‑36 (MQ=255)
            cGCAATGATGATCACCAGTCCGGCAATACAGCCAAt                         >  1:5173169/1‑36 (MQ=255)
            cGCAATGATGATCACCAGTCCGGCAATACAGCCAAt                         >  1:2528108/1‑36 (MQ=255)
             gCAATGATGATCACCAGTCCGGCAATACAGCCaaca                        >  1:5231054/1‑34 (MQ=255)
             gCAATGATGATCACCAGTCCGGCAATACAGCCAATa                        >  1:4972105/1‑36 (MQ=255)
             gCAATGATGATCACCAGTCCGGCAATACAGCCAATa                        >  1:285799/1‑36 (MQ=255)
              cAATGATGATCACCAGTCCGGCAATACAGCCAATaa                       >  1:1026963/1‑36 (MQ=255)
              cAATGATGATCACCAGTCCGGCAATACAGCCAATaa                       >  1:5751827/1‑36 (MQ=255)
              cAATGATAACCACCAGTCCGGCAATACAGCCAATaa                       <  1:4690102/36‑1 (MQ=37)
               aaTGATGATCACCAGTCCGGCAATACAGCCAATaaa                      >  1:890270/1‑36 (MQ=255)
               aaTGATGATCACCAGTCCGGCAATACAGCCAATaaa                      >  1:2539179/1‑36 (MQ=255)
                atgatgATCACCAGTCCGGCAATACAGCCAATAAAt                     >  1:7841/1‑36 (MQ=255)
                atgatgATCACCAGTCCGGAAATACAGCCAATAAAt                     <  1:3327844/36‑1 (MQ=255)
                aaaatgatCACAAGTCCGGAAATACAGCCAATAAAt                     <  1:1049617/33‑1 (MQ=21)
                 tgatgaTCACCAGTCCGGCAATACAGCCAATAAAtg                    >  1:5247317/1‑36 (MQ=255)
                  gatgatCACCAGTCCGGCAATACAGCCAATAAAtgt                   >  1:7294415/1‑36 (MQ=255)
                   atgatCACCAGTCCGGCAATACAGCCAATAAAtgtg                  >  1:1829941/1‑36 (MQ=255)
                   atgatCACAAGTCCGCAAATACAGCCAATAAAtgtg                  <  1:6082329/36‑1 (MQ=37)
                    tgatCACCAGTCCGGCAATACAGCCAATAAATGTGc                 >  1:3197902/1‑36 (MQ=255)
                    tgatCACCAGTCCGGCAATACAGCCAATAAATGTGc                 >  1:5172152/1‑36 (MQ=255)
                     gatCACCAGTCCGGCAATACAGCCAATAAATGTGcc                >  1:2640493/1‑36 (MQ=255)
                      atCACCAGTCCGGCAATACAGCCAATAAATGTGCCg               >  1:6198453/1‑36 (MQ=255)
                       tCACCAGTCCGGCAATACAGCCAATAAATGTGCCga              >  1:490945/1‑36 (MQ=255)
                       tCACCAGTCAGGAAATACAGCCAATAAATGTGCCga              <  1:6666042/36‑1 (MQ=37)
                       acACAAGTCCGTAAATCCAGCCAATAAATGTGCCga              <  1:4009852/35‑1 (MQ=16)
                        cACCAGTCCGGCAATACAGCCAATAAATGTGCCgat             >  1:2041743/1‑36 (MQ=255)
                        cACCAGTCCGGAAATACAGCCAATAAATGTGCCgat             <  1:4700231/36‑1 (MQ=255)
                        cACAAGTCCGGAAATACAGCCAATAAATGTGCCgat             <  1:7013904/36‑1 (MQ=37)
                         aCCAGTCCGGCAATACAGCCAATAAATGTGCCgatg            >  1:5628968/1‑36 (MQ=255)
                         aCCAGTCCGGCAATACAGCCAATAAATGTGCCgatg            >  1:4653298/1‑36 (MQ=255)
                         aCCAGTCCGGAAATACAGCCAATAAATGTGCCgatg            <  1:5082163/36‑1 (MQ=255)
                         aCCAGTCCGCAAATCCAGCCAATAAATGTGCCgatg            <  1:1266630/36‑1 (MQ=25)
                          caaGTCCGGAAATACAGCCAATAAATGTGCCgatga           <  1:4911555/34‑1 (MQ=38)
                           gaGTCCGGAAATACAGCCAATAAATGTGCCgatgat          <  1:6865781/35‑1 (MQ=255)
                           cAGTCCGGCAATACAGCCAATAAATGTGACgatgat          >  1:5931456/1‑36 (MQ=255)
                           aaGTCCGGAAATACAGCCAATAAATGTGCCgatgat          <  1:2552366/35‑1 (MQ=255)
                             gTCCGGCAATACAGCCAATAAATGTGCCGATGATgc        >  1:1774434/1‑36 (MQ=255)
                             gTCCGCAAATACAGCCAATAAATGTGCCGATGATgc        <  1:5477138/36‑1 (MQ=37)
                              tCCGGCAATACAGCCAATAAATGTGCCGATGATgcg       >  1:4661940/1‑36 (MQ=255)
                              tCCGGCAATACAGCCAATAAATGTGCCGATGATgcg       >  1:2111921/1‑36 (MQ=255)
                               ccGGAAACACAGCCAATAAATGTGCCGATGATgcgc      <  1:2941903/36‑1 (MQ=37)
                                cGGCAATACAGCCAATAAATGTGCCGATGATGCGCa     >  1:352431/1‑36 (MQ=255)
                                 ggCAATACAGCCAATAAATGTGCCGATGATGCGGaa    >  1:1312994/1‑36 (MQ=255)
                                  gCAATACAGCCAATAAATGTGCCGATGATGCGCaaa   >  1:6737264/1‑36 (MQ=255)
                                   cAATACAGCCAATAAATGTGCCGATGATGCGCaaaa  >  1:1413124/1‑36 (MQ=255)
                                   cAATACAGCCAATAAATGTGCCGATGATGCGCaaaa  >  1:1805474/1‑36 (MQ=255)
                                   |                                   
TGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAA  >  REL606/3315813‑3315883

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