Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,139,263 +T 6.1% coding (913/918 nt) oxyR → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,139,2631.→T 6.1% 3.9 ‑1.6 ‑0.233coding (913/918 nt)oxyRDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (13/18):  total (15/18)
Fisher's exact test strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.76e-01

TGGATGGCCATTTCGATAAAGTTTTAAAACAGGC.GG.TTTAAACCGTTTAACGCAGCTACCCGATAGGCTTC  >  REL606/4139228‑4139298
                                     |                                   
tggatggCCATTTCGATAAAGTTTTAAAACAGGC.gg                                      >  1:2926123/1‑36 (MQ=255)
  gatggCCATTTCGATAAAGTTTTAAAACAGGC.gggt                                    >  1:3924901/1‑34 (MQ=255)
  gagggCCATTTCCATAAAGTGTTAAAACAGGC.GG.tt                                   >  1:3187345/4‑36 (MQ=25)
    tggCCATTTCGATAAAGTTTTAAAACAGGC.GG.TTTa                                 <  1:692279/36‑1 (MQ=255)
     ggCCATTTCGATAAAGTTTTAAAACAGGC.GG.gtttt                                >  1:835079/1‑34 (MQ=25)
     ggCCATTTCGATAAAGTTTTAAAACAGGC.GG.GTTat                                >  1:4103213/1‑35 (MQ=255)
        cATTTCGATAAAGTTTTAAAACAGGC.GG.TTTAAAcc                             <  1:996112/36‑1 (MQ=255)
         aTTTCGATAAAGTTTTAAAACAGGC.GG.TTTTAAACg                            >  1:1942063/1‑36 (MQ=37)
         aTTTCGATAAAGTTTTAAAACAGGC.GG.TTTAAAACg                            >  1:610073/1‑36 (MQ=255)
          tttCGATAAAGTTTTAAAACATGC.GGTTTTAAAcct                            >  1:4946915/1‑35 (MQ=21)
          tttCGATAAAGTTTTAAAACAGGCGGG.TTTaaaact                            >  1:3066049/1‑33 (MQ=25)
          tttCGATAAAGTTTTAAAACAGGC.GG.TTTAAACCGt                           <  1:39085/36‑1 (MQ=255)
           ttCGATAAAGTTTTAAAACAGGC.GG.TTTAAACCGtt                          <  1:3690037/36‑1 (MQ=255)
            tCGATAAAGTTTTAAAACAGGC.GG.TTTAAACCGttt                         <  1:3244005/36‑1 (MQ=255)
             cGATAAAGTTTTAAAACAGGC.GG.TTTAAACCGTTTa                        <  1:4445382/36‑1 (MQ=255)
              gATAAAGTTTTAAAACAGGCGGG.TTTAAAGCTtttt                        >  1:2706428/1‑35 (MQ=16)
              gATAAAGTTTTAAAACAGGC.GG.TTTAAACCGTTTaa                       <  1:2792730/36‑1 (MQ=255)
                tAAAGTTTTAAAACAGGC.GG.TTTAAACCCTTTAACg                     >  1:5789111/1‑36 (MQ=255)
                 aaaGTTTTAAAACAGGC.GGTTTTAAACCCTTAaatc                     >  1:3626338/1‑34 (MQ=12)
                  aaGTTTTAAAACAGGC.GG.TTTAAACCGTTTAACGCa                   <  1:2809624/36‑1 (MQ=255)
                  aaGTTTTAAAACAGGC.GG.TTTAAACCGTTTAACGCa                   <  1:2738895/36‑1 (MQ=255)
                   aGTTTTAAAACAGGC.GG.TTTAAACCGTTTAACGCAg                  <  1:4205159/36‑1 (MQ=255)
                     ttttAAAACAGGC.GG.TTTAAACCGTTTAACGCAGCt                <  1:2361602/36‑1 (MQ=255)
                     ttttAAAACAGGC.GG.TTTAAACCGTTTAACGCAGCt                <  1:1269026/36‑1 (MQ=255)
                     ttttAAAACAGGC.GG.TTTAAACCGTTTAACGCAGCt                <  1:846716/36‑1 (MQ=255)
                       ttAAAACAGGC.GG.TTTAAACCGTTTAACGCAGCTAc              <  1:4419617/36‑1 (MQ=255)
                       ttAAAACAGGC.GG.GTTTAAACGTTTTACCCAGCttc              >  1:4141570/1‑34 (MQ=9)
                         aaaaCAGGC.GG.TTTAAACCGTTTAACGCAGCTAccc            <  1:3564178/36‑1 (MQ=255)
                         aaaaCAGGC.GG.TTTAAAAAGCTTAAATCAGCTAccc            >  1:2974480/1‑36 (MQ=16)
                           aaCAGGC.GG.TTTTAACAGTTTAAAGCAGGTAcccca          >  1:1574980/1‑34 (MQ=16)
                           aaCAGGC.GG.TTTAAACCGTTTTTGGCAGCTTACCGa          >  1:4366107/1‑36 (MQ=12)
                            aaaGGC.GG.TTTAAACCGTTTAACGCAGCTACCCGAt         <  1:2989899/34‑1 (MQ=255)
                            aCAGGC.GG.TTTAAACCGTTTAACGCAGCTACCCGAt         <  1:123230/36‑1 (MQ=255)
                              aGGC.GG.TTTAAACCGTTTAACGCAGCTACCCGATAg       <  1:5302663/36‑1 (MQ=255)
                              aGGC.GG.TTTAAACCGTTTAACGCAGCTACCCCATAg       >  1:4744011/1‑36 (MQ=255)
                               ggC.GG.TTTAAACCGTTTAACGAAGCTACCCGATAgg      >  1:3919636/1‑36 (MQ=255)
                                gC.GG.TTTAAACCGTTTAACGCAGCTACCCGATAGGc     >  1:5550607/1‑36 (MQ=255)
                                   gg.TTTAAACCGTTTAACGCAGCTACCCGATAGGCtt   >  1:3428675/1‑36 (MQ=255)
                                   gg.TTTAAAACGTTTAACGCAGCTACCCGATAGGCtt   <  1:1347608/36‑1 (MQ=255)
                                    g.TTTAAACCGTTTAACGCAGCTACCCGATAGGCTTc  >  1:3867944/1‑36 (MQ=255)
                                    g.TTTAAACCGTTTAACGCAGCTACCCGATAGGCTTc  >  1:3090950/1‑36 (MQ=255)
                                     |                                   
TGGATGGCCATTTCGATAAAGTTTTAAAACAGGC.GG.TTTAAACCGTTTAACGCAGCTACCCGATAGGCTTC  >  REL606/4139228‑4139298

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 13 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG