Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,584,6150T→A 8.8% 3.0 ‑6.0 ‑0.070442yhiFpredicted DNA‑binding ranscriptional regulator
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (42/21):  total (42/28)
Fisher's exact test strand distribution p-value = 2.47e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.89e-01

TTCTGGAGAAACTGAAAATATCGCTAAAAACGTTCT.ACTGCCATAAACACAATATCATGATGATCCTCAAT  >  REL606/3584580‑3584650
                                   |                                    
ttCTGGAGAAACTGAAAATATCGCTAAAAACGTTCt                                      >  1:7755844/1‑36 (MQ=255)
ttCTGGAGAAACTGAAAATATCGCTAAAAACGTTCt                                      <  1:3132517/36‑1 (MQ=255)
  cTGGAGAAACTGAAAATATCGCTAAAAACGTTCT.Ac                                   >  1:7374784/1‑36 (MQ=255)
   tGGAGAAACTGAAAATATCGCTAAAAACGTTCT.ACt                                  <  1:4171625/36‑1 (MQ=255)
   tGGAGAAACTGAAAATATCGCTAAAAACGTTCT.ACt                                  <  1:7784131/36‑1 (MQ=255)
   tGGAGAAACTGAAAATATCGCTAAAAACGTTCT.ACt                                  >  1:3111308/1‑36 (MQ=255)
    ggAGAAACTGAAAATATCGCTAAAAACGTTCT.ACTg                                 >  1:4396345/1‑36 (MQ=255)
    ggAGAAACTGAAAATATCGCTAAAAACGTTCT.ACTg                                 >  1:3756322/1‑36 (MQ=255)
    ggAGAAACTGAAAATATCGCTAAAAACGTTCT.ACTg                                 >  1:2652259/1‑36 (MQ=255)
    ggAGAAACTGAAAATATCGCTAAAAACGTTCT.ACTg                                 <  1:7906155/36‑1 (MQ=255)
     gagaAACTGAAAATATCGCTAAAAACGTTCT.ACTGc                                <  1:7423724/36‑1 (MQ=255)
     gagaAACTGAAAATATCGCTAAAAACGTTCT.ACTGc                                <  1:4398514/36‑1 (MQ=255)
      agaAACTGAAAATATCGCTAAAAACGTTCT.ACTGcc                               <  1:3520421/36‑1 (MQ=255)
      agaAACTGAAAATATCGCTAAAAACGTTCT.ACTGcc                               >  1:5563895/1‑36 (MQ=255)
      agaAACTGAAAATATCGCTAAAAACGTTCT.ACTGcc                               >  1:3584173/1‑36 (MQ=255)
      agaAACTGAAAATATCGCTAAAAACGTTCT.ACTGcc                               <  1:2338725/36‑1 (MQ=255)
       gaAACTGAAAATATCGCTAAAAACGTTCT.ACTGCCa                              <  1:3239153/36‑1 (MQ=255)
       gaAACTGAAAATATCGCTAAAAACGTTCT.ACTGCCa                              <  1:7833780/36‑1 (MQ=255)
       gaAACTGAAAATATCGCTAAAAACGTTCT.ACTACCa                              <  1:1716356/36‑1 (MQ=255)
        aaaCTGAAAATATCGCTAAAAACGTTCT.ACTGCCAt                             >  1:7564904/1‑36 (MQ=255)
        aaaCTGAAAATATCGCAAAAACCGTTCT.ACTGCCAt                             <  1:3995003/36‑1 (MQ=37)
         aaCTGAAAATATCGCTAAAAACGTTCT.ACTGCCATa                            >  1:6224798/1‑36 (MQ=255)
         aaCTGAAAATATCGCTAAAAACGTTCT.ACTGCCATa                            >  1:3749693/1‑36 (MQ=255)
         aaCTGAAAATATCGCTAAAAACGTTCT.ACTGCCATa                            >  1:2314081/1‑36 (MQ=255)
          aCTGAAAATATCGCTAAAAACGTTCT.ACTGCCATaa                           >  1:6700385/1‑36 (MQ=255)
          aCTGAAAATATCGCTAAAAACGTTCT.ACTGCCATaa                           >  1:7491716/1‑36 (MQ=255)
            taaaaaTATCGCAAAAAACGTTCT.ACTGCCATAAac                         <  1:2472889/34‑1 (MQ=37)
            taaaaaTATCGAAAAAAACGTTCT.ACTGCCATAAac                         <  1:5324451/34‑1 (MQ=37)
            taaaaaTATCAAAAAAAACGTTCT.ACTGCCATAAac                         <  1:7571151/34‑1 (MQ=21)
            tGAAAATATCGCTAAAAACGTTCT.ACTGCCATAAac                         >  1:1703184/1‑36 (MQ=255)
            gaaaaaTATCGTAAAAAACTTTCT.ACTGCCATAAac                         <  1:6954521/34‑1 (MQ=18)
             gAAAATATCGCTAAAAACGTTCT.ACTGCCATAAaca                        >  1:3604736/1‑36 (MQ=255)
             gAAAATATCGCTAAAAACGTTCT.ACTGCCATAAaca                        >  1:5088198/1‑36 (MQ=255)
             gAAAATATCGCTAAAAACGTTCT.ACTGCCATAAaca                        >  1:89324/1‑36 (MQ=255)
             gAAAATATCGCAAAAAACGTTCT.ACTGCCATAAaca                        <  1:5933026/36‑1 (MQ=255)
             aaaaaTATCGCAAAAAACGTACT.ACTGCCATAAaca                        <  1:7370852/35‑1 (MQ=37)
              aaaaTATCGCTAAAAACGTTCT.ACTGCCATAAacac                       >  1:7170607/1‑36 (MQ=255)
              aaaaTATCGCTAAAAACGTTCT.ACTGCCATAAacac                       >  1:4068882/1‑36 (MQ=255)
              aaaaTATCGCTAAAAACGTTCT.ACTGCCATAAacac                       >  1:1517050/1‑36 (MQ=255)
              aaaaTATCGAAAAAACCGTTCT.ACTGCCATAAacac                       <  1:583436/36‑1 (MQ=25)
                 atatCGCTAAAAACGTTCT.ACTGCCATAAACACAat                    >  1:7168442/1‑36 (MQ=255)
                 atatCGCTAAAAACGTTCT.ACTGCCATAAACACAat                    >  1:2854618/1‑36 (MQ=255)
                  tatCGCTAAAAACGTTCT.ACTGCCATAAACACAata                   >  1:3378800/1‑36 (MQ=255)
                      gCTAAAAACGTTCT.ACTGCCATAAACACAATATCat               >  1:4313285/1‑36 (MQ=255)
                      gCTAAAAACGTTCT.ACTGCCATAAACACAATATCat               >  1:5021763/1‑36 (MQ=255)
                      gCTAAAAACGTTCT.ACTGCCATAAACACAATATCat               >  1:2807292/1‑36 (MQ=255)
                      aaaaaaaaCTTTCA.ACTGCCATAAACACAATATCat               <  1:5230052/33‑1 (MQ=25)
                       aaaaaaaAGTTCCAACTGCCATAAACACAATATCat               <  1:3993326/34‑1 (MQ=16)
                        tAAAAACGTTCT.ACTGCCATAAACACAATATCatga             >  1:4523575/1‑36 (MQ=255)
                        aaaaaaAGTCCA.ACTGCCATAAACACAATATCatga             <  1:2347215/35‑1 (MQ=21)
                         aaaaaCGTTCT.ACTGCCATAAACACAATATCatgat            >  1:3667611/1‑36 (MQ=255)
                         aaaaaCGTTCT.ACTGCCATAAACACAATATCatgat            >  1:2152453/1‑36 (MQ=255)
                         aaaaaCGTTCT.ACTGCCATAAACACAATATCatgat            <  1:2296317/36‑1 (MQ=255)
                         aaaaaCGTTAT.ACTGCCATAAACACAATATCatgat            <  1:155222/36‑1 (MQ=255)
                          aaaaTGTGCA.ACTGCAATAAACACAATATCatgatg           <  1:5059448/36‑1 (MQ=16)
                           aaaCTTTCA.ACTGCCATAAACACAATATCAtgatga          <  1:6985016/36‑1 (MQ=37)
                           aaaCGTTCT.ACTGCCATAAACACAATATCAtgatga          >  1:3737398/1‑36 (MQ=255)
                           aaaCGTTCT.ACTGCCATAAACACAATATCAtgatga          <  1:339581/36‑1 (MQ=255)
                            aaCGTTAA.ACTGCCATAAACACAATATCATgatgat         <  1:5722397/36‑1 (MQ=37)
                              cGTTCT.ACTGCCATAAACACAATATCATGATGATcc       >  1:2034543/1‑36 (MQ=255)
                              cGTTCT.ACTGCCATAAACACAATATCATGATGATcc       >  1:2202422/1‑36 (MQ=255)
                               gTTCT.ACTGCCATAAACACAATATCATGATGATCCt      >  1:2014331/1‑36 (MQ=255)
                               gTTCT.ACTGCCATAAACACAATATCATGATGATCCt      >  1:5416420/1‑36 (MQ=255)
                               gTTCT.ACTGCCATAAACACAATATCATGATGATCCt      >  1:1220714/1‑36 (MQ=255)
                               gTTCT.ACTGCCATAAACACAATATCATGATGATCCt      >  1:1068112/1‑36 (MQ=255)
                               gTTCA.ACTGCCATAAACACAATATCATGATGATCCt      <  1:1197236/36‑1 (MQ=255)
                                ttCT.ACTGCCATAAACACAATATCATGATGATCCTc     >  1:1474309/1‑36 (MQ=255)
                                ttCT.ACTGCCATAAACACAATATCATGATGATCCTc     >  1:852313/1‑36 (MQ=255)
                                 tCT.ACTGCCATAAACACAATATCATGATGATCCTCa    >  1:7139187/1‑36 (MQ=255)
                                 tCT.ACTGCCATAAACACAATATCATGATGATCCTCa    >  1:2976721/1‑36 (MQ=255)
                                  cT.ACTGCCATAAACACAATATCATGATGATCCTCaa   >  1:7315283/1‑36 (MQ=255)
                                  cT.ACTGCCATAAACACAATATCATGATGATCCTCaa   >  1:6758426/1‑36 (MQ=255)
                                   t.ACTGCCATAAACACAATATCATGATGATCCTCAAt  >  1:1805330/1‑36 (MQ=255)
                                   |                                    
TTCTGGAGAAACTGAAAATATCGCTAAAAACGTTCT.ACTGCCATAAACACAATATCATGATGATCCTCAAT  >  REL606/3584580‑3584650

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG