Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,139,261 +G 5.7% coding (911/918 nt) oxyR → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,139,2611.→G 5.7% 3.1 ‑1.5 ‑0.135coding (911/918 nt)oxyRDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (15/18):  total (17/18)
Fisher's exact test strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.01e-01

AATGGATGGCCATTTCGATAAAGTTTTAAAACAGGC.GG.TTTAAACCGTTTAACGCAGCTACCCGATAGGC  >  REL606/4139226‑4139295
                                    |                                   
aaTGGATGGCCATTTCGATAAAGTTTTAAAACAGGc                                      >  1:1708044/1‑36 (MQ=255)
aaTGGATGGCCATTTCGATAAAGTTTTAAAACAGGc                                      <  1:1947533/36‑1 (MQ=255)
aaTGGATGGCCATTTCGATAAAGTTTTAAAACAGGc                                      >  1:2118862/1‑36 (MQ=255)
aaTGGATGGCCATTTATAGAAAGATTTAAAACACGc                                      >  1:1387781/1‑36 (MQ=16)
  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:846716/36‑1 (MQ=255)
                       ttttAAAACAGGC.GG.TTTAAACCGTTTAACGCAGCt             <  1:1269026/36‑1 (MQ=255)
                       ttttAAAACAGGC.GG.TTTAAACCGTTTAACGCAGCt             <  1:2361602/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)
                                    |                                   
AATGGATGGCCATTTCGATAAAGTTTTAAAACAGGC.GG.TTTAAACCGTTTAACGCAGCTACCCGATAGGC  >  REL606/4139226‑4139295

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 14 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG