Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,584,6040T→A 7.0% 2.1 ‑6.4 ‑0.070431yhiFpredicted DNA‑binding ranscriptional regulator
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (49/16):  total (49/21)
Fisher's exact test strand distribution p-value = 1.68e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.98e-01

AACGGAAGATATTCTGGAGAAACTGAAAATATCGCTAAAAACGTTCTACTGCCATAAACACAATATCATGA  >  REL606/3584569‑3584639
                                   |                                   
aaCGGAAGATATTCTGGAGAAACTGAAAATATCGCt                                     <  1:3472504/36‑1 (MQ=255)
 aCGGAAGATATTCTGGAGAAACTGAAAATATCGCTa                                    >  1:2675998/1‑36 (MQ=255)
 aCGGAAGATATTCTGGAGAAACTGAAAATATCGCTa                                    >  1:6604344/1‑36 (MQ=255)
 aCGGAAGAGATTCTGGAGAAACTGAAAATATCGCTa                                    <  1:6086686/36‑1 (MQ=255)
   gtaaGATATTCTGGAGAAACTGAAAATATCGCTaaa                                  <  1:5509537/34‑1 (MQ=255)
   ggAAGATATTCTGGAGAAACTGAAAATATCGCTaaa                                  >  1:2937580/1‑36 (MQ=255)
   ggAAGATATTCTGGAGAAACTGAAAATATCGCTaaa                                  >  1:3943199/1‑36 (MQ=255)
   ggAAGATATTCTGGAGAAACTGAAAATATCGCTaaa                                  >  1:4551059/1‑36 (MQ=255)
   ggAAGATATTCTGGAGAAACTGAAAATATCGCTaaa                                  >  1:2078957/1‑36 (MQ=255)
   ggAAGATATTCTGGAGAAACTGAAAATATCGCTaaa                                  <  1:1791646/36‑1 (MQ=255)
    gAAGATATTCTGGAGAAACTGAAAATATCGCTaaaa                                 >  1:3886160/1‑36 (MQ=255)
    gAAGATATTCTGGAGAAACTGAAAATATCGCTaaaa                                 >  1:6957131/1‑36 (MQ=255)
    gAAGATATTCTGGAGAAACTGAAAATATCGCTaaaa                                 <  1:2584657/36‑1 (MQ=255)
     aaGATATTCTGGAGAAACTGAAAATATCGCTaaaaa                                >  1:359825/1‑36 (MQ=255)
     aaGATATTCTGGAGAAACTGAAAATATCGCTaaaaa                                >  1:147542/1‑36 (MQ=255)
     aaGATATTCTGGAGAAACTGAAAATATCGCTaaaaa                                >  1:4943866/1‑36 (MQ=255)
     aaGATATTCTGGAGAAACTGAAAATATCGCTaaaaa                                >  1:4287149/1‑36 (MQ=255)
      aGATATTCTGGAGAAACTGAAAATATCGCTAAAAAc                               <  1:981933/36‑1 (MQ=255)
      aGATATTCTGGAGAAACTGAAAATATCGCTAAAAAc                               <  1:4799346/36‑1 (MQ=255)
      aGATATTCTGGAGAAACTGAAAATATCGCTAAAAAc                               >  1:1526290/1‑36 (MQ=255)
       gATATTCTGGAGAAACTGAAAATATCGCTAAAAACg                              >  1:7660349/1‑36 (MQ=255)
       gATATTCTGGAGAAACTGAAAATATCGCTAAAAACg                              >  1:6132742/1‑36 (MQ=255)
        atatTCTGGAGAAACTGAAAATATCGCTAAAAACGt                             >  1:6616577/1‑36 (MQ=255)
        atatTCTGGAGAAACTGAAAATATCGCTAAAAACGt                             <  1:3865640/36‑1 (MQ=255)
          atTCTGGAGAAACTGAAAATATCGCTAAAAACGTTc                           <  1:7587091/36‑1 (MQ=255)
          atTCTGGAGAAACTGAAAATATCGCTAAAAACGTTc                           <  1:3984101/36‑1 (MQ=255)
          atTCTGGAGAAACTGAAAATATCGCTAAAAACGTTc                           <  1:827901/36‑1 (MQ=255)
          atTCTGGAGAAACTGAAAATATCGCTAAAAACGTTc                           >  1:2311943/1‑36 (MQ=255)
           ttCTGGAGAAACTGAAAATATCGCTAAAAACGTTCt                          >  1:7533433/1‑36 (MQ=255)
           ttCTGGAGAAACTGAAAATATCGCTAAAAACGTTCt                          <  1:1428112/36‑1 (MQ=255)
            tCTGGAGAAACTGAAAATATCGCTAAAAACGTTCTa                         >  1:1161809/1‑36 (MQ=255)
            tCTGGAGAAACTGAAAATATCGCTAAAAACGTTCTa                         <  1:4931725/36‑1 (MQ=255)
             cTGGAGAAACTGAAAATATCGCTAAAAACGTTCTAc                        <  1:1683378/36‑1 (MQ=255)
             cTGGAGAAACTGAAAATATCGCTAAAAACGGTCTAc                        >  1:2991814/1‑36 (MQ=255)
              tGGAGAAACTGAAAATATCGCTAAAAACGTTCTACt                       >  1:6322484/1‑36 (MQ=255)
              tGGAGAAACTGAAAATATCGCTAAAAACGTTCTACt                       <  1:595714/36‑1 (MQ=255)
              tGGAGAAACTGAAAATATCGCTAAAAACGTTCTACt                       <  1:2467638/36‑1 (MQ=255)
               ggAGAAACTGAAAATATCGCTAAAAACGTTCTACTg                      <  1:4181751/36‑1 (MQ=255)
               ggAGAAACTGAAAATATCGCTAAAAACGTTCTACTg                      >  1:5306294/1‑36 (MQ=255)
                 agaAACTGAAAATATCGCTAAAAACGTTCTACTGcc                    >  1:5868849/1‑36 (MQ=255)
                  gaAACTGAAAATATCGCTAAAAACGTTCTACTGCCa                   >  1:372292/1‑36 (MQ=255)
                  gaAACTGAAAATATCGCTAAAAACGTTCTACTGCCa                   >  1:6396713/1‑36 (MQ=255)
                   aaaCTGAAAATATCGCTAAAAACGTTCTATTGCCAt                  >  1:1040020/1‑36 (MQ=255)
                   aaaCTGAAAATATCGCTAAAAACGTTCTACTGCCAt                  >  1:3241532/1‑36 (MQ=255)
                   aaaCTGAAAATATCGCTAAAAACGTTCTACTGCCAt                  >  1:7178199/1‑36 (MQ=255)
                   aaaCTGAAAATATCGCTAAAAACGTTCTACTGCCAt                  >  1:1435211/1‑36 (MQ=255)
                     aCTGAAAATATCGCTAAAAACGTTCTACTGCCATaa                >  1:1044063/1‑36 (MQ=255)
                     aCTGAAAATATCGCTAAAAACGTTCTACTGCCATaa                >  1:1897102/1‑36 (MQ=255)
                     aCTGAAAATATCGCTAAAAACGTTCTACTGACATaa                >  1:2800067/1‑36 (MQ=255)
                       tGAAAATATCGCTAAAAACGTTCTACTGCCATAAac              >  1:5822177/1‑36 (MQ=255)
                       tGAAAATATCGCTAAAAACGTTCTACTGCCATAAac              >  1:6740144/1‑36 (MQ=255)
                       aaaaaaCATCCCAAAAACCGTTCAACTGCCATAAac              <  1:2476702/34‑1 (MQ=12)
                        gAAAATATCGCTAAAAACGTTCTACTGCCATaaaaa             >  1:3723184/1‑34 (MQ=255)
                        gAAAATATCGCAAAAAACGTTCTACTGCCATAAaca             <  1:2731919/36‑1 (MQ=255)
                         aaaaTATCGCTAAAAACGTTCTACTGCCATAAacac            >  1:5201676/1‑36 (MQ=255)
                         aaaaTATCGCTAAAAACGTTCTACTGCCATAAacac            >  1:3478030/1‑36 (MQ=255)
                         aaaaTATCGCTAAAAACGTTCTACTGCCATAAacac            >  1:2617281/1‑36 (MQ=255)
                         aaaaTATCCAAAAAAAGATTCCACTGCCATAAacac            <  1:3973101/36‑1 (MQ=9)
                          aaaTATCGCTAAAAACGTTCTACTTCCATAAAcaca           >  1:6780863/1‑36 (MQ=255)
                          aaaTATCGCTAAAAACGTTCTACTGCCATAAAcaca           >  1:7171389/1‑36 (MQ=255)
                          aaaTATCGAAAAAAACGTTCTACTGCCATAAAcaca           <  1:7616051/36‑1 (MQ=38)
                           aaTATCGAAAAAAACGTTCTACTGCCATAAACACaa          <  1:5408286/36‑1 (MQ=38)
                            atatCGCTAAAAACGTTCTACTGCCAGAAACACAat         >  1:7014364/1‑36 (MQ=255)
                              atCGCTAAAAACGTTCTACTGCCATAAACACAatat       >  1:1590894/1‑36 (MQ=255)
                              atCGCTAAAAACGTTCTACTGCCATAAACACAatat       >  1:1042337/1‑36 (MQ=255)
                                cGCTAAAAACGTTCTACTGCCATAAACACAATATCa     >  1:6764257/1‑36 (MQ=255)
                                cGCTAAAAACGTTCTACTGCCATAAACACAATATCa     >  1:2821032/1‑36 (MQ=255)
                                agCTAAAAAAGTTCTACTGAAATAAACACAATATCa     >  1:196330/2‑36 (MQ=25)
                                 gCTAAAAACGTTCTACTGCCATAAACACAATATCat    >  1:6467182/1‑36 (MQ=255)
                                  cTAAAAACGTTCTACTGCCATAAACACAATATCatg   >  1:6090050/1‑36 (MQ=255)
                                  cTAAAAACGTTCTACTGCCATAAACACAATATCatg   >  1:7162272/1‑36 (MQ=255)
                                   tAAAAACGTTCTACTGCCATAAACACAATATCatga  >  1:2722985/1‑36 (MQ=255)
                                   tAAAAACGTTCTACTGCCATAAACACAATATCatga  >  1:5496776/1‑36 (MQ=255)
                                   |                                   
AACGGAAGATATTCTGGAGAAACTGAAAATATCGCTAAAAACGTTCTACTGCCATAAACACAATATCATGA  >  REL606/3584569‑3584639

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