Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,236,623 +T 5.7% coding (577/609 nt) lexA → LexA repressor

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,236,6231.→T 5.7% 3.4 ‑2.4 ‑0.235coding (577/609 nt)lexALexA repressor
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (9/24):  total (11/24)
Fisher's exact test strand distribution p-value = 9.24e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.07e-01

TTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGC.GG.TTG.GGGTTATTCGCAACGGCGACTGGCTGTAACAT  >  REL606/4236588‑4236658
                                     |                                    
ttCGTCAGCAGAGCTTCACCATTGAAGGGCTGGC.gg                                       <  1:130328/36‑1 (MQ=255)
ttCGTCAGCAGAGCTTCACCATTGAAGGGCTGGC.gg                                       >  1:776979/1‑36 (MQ=255)
ttCGTCAGCAGAGCTTCACCATTGAAGGGCTGGC.gg                                       >  1:5741155/1‑36 (MQ=255)
 tCGTCAGCAGAGCTTCACCATTGAAGGGCTGGC.GG.t                                     <  1:5628202/36‑1 (MQ=255)
 tCGTCAGCAGAGCTTCACCATTGAAGGGCTGGC.GG.t                                     <  1:5694808/36‑1 (MQ=255)
  cGTCAGCAGAGCTTCACCATTGAAGGGCTGGC.GG.tt                                    <  1:3633535/36‑1 (MQ=255)
  cGTCAGCAGAGCTTCACCATTGAAGGGCTGGC.GG.tt                                    <  1:7410856/36‑1 (MQ=255)
   gccagcagAGCTTCACCATTGAAGGGCTGGC.GG.TTg                                   <  1:7930364/34‑1 (MQ=255)
   gTCAGCAGAGCTTCACCATTGAAGGGCTGGC.GG.TTg                                   <  1:2811325/36‑1 (MQ=255)
     cagcagAGCTTCACCATTGAAGGGCTGGC.GG.GTg.gg                                >  1:1894523/1‑36 (MQ=255)
      agcagAGCTTCACCATTGAAGGGCTGGC.GG.TTg.ggg                               <  1:1264730/36‑1 (MQ=255)
      agcagAGCTTCACCATTGAAGGGCTGGC.GG.GTt.ggg                               >  1:1940621/1‑36 (MQ=37)
      agcagAGCTTCACCATTGAAGGGCTGGC.GG.ATt.ggg                               >  1:5957940/1‑36 (MQ=37)
      agcagAGCTTCACCATTGAAGGGCTGGC.GG.ATg.ggg                               >  1:233490/1‑36 (MQ=255)
      agaagaaCTTCACCATTGAAGGGCTGGC.GG.TTt.ggg                               <  1:5774753/33‑1 (MQ=21)
       gcagAGCTTCACCATTGAAGGGCTGGC.GG.TTG.GGGt                              <  1:3321618/36‑1 (MQ=255)
         agagCTTCACCATTGAAGGGCTGGC.GG.TTG.GGGTTa                            <  1:5661476/36‑1 (MQ=255)
          gagCTTCACCATTGAAGGGCTGGG.GG.TTG.GGGTTAt                           <  1:7506368/36‑1 (MQ=255)
          gagCTTCACCATTGAAGGGCTGGC.GG.TTGTGGGttt                            >  1:5139953/1‑35 (MQ=37)
          gagCTTCACCATTGAAGGGCTGGC.GG.GTG.GGGttgt                           >  1:5892819/1‑34 (MQ=37)
           agCTTCACCATTGAAGGGCTGGCGGG.TTC.GGGtttt                           >  1:5211897/1‑34 (MQ=18)
           agCTTCACCATTGAAGGGCTGGC.GGTTTG.GGGtttt                           >  1:3880690/1‑34 (MQ=25)
           agCTTCACCATAGAAGGGCTGGC.GG.TTG.GGGTTAtt                          <  1:1755491/36‑1 (MQ=255)
            gCTTCACCATTGAAGGGCTGGC.GGTTTG.GGGTTTtt                          >  1:4390369/1‑36 (MQ=21)
            gCTTCACCATTGAAGGGCTGGC.GG.TTG.GGGTTATTc                         <  1:6306049/36‑1 (MQ=255)
             cTTCACCATTGAAGGGCTGGC.GG.TTG.GGGTTATTCg                        <  1:2745634/36‑1 (MQ=255)
               tCACCATTGAAGGGCTGGC.GG.TTG.GGGTTATTCGCa                      <  1:2586377/36‑1 (MQ=255)
                cACCATTGAAGGGCTGGC.GT.TTG.GGGTTATTCGCaa                     <  1:5427494/36‑1 (MQ=255)
                    aTTGAAGGGCTGGC.GG.TTG.GGGTTATTCGCAACGGc                 <  1:7125413/36‑1 (MQ=255)
                     ttGAAGGGCTGGC.GG.TTG.GGGTTATTCGCAACGGCg                <  1:3137008/36‑1 (MQ=255)
                     ttGAAGGGCTGGC.GG.TTG.GGGTTATTCGCAACGGCg                <  1:1717161/36‑1 (MQ=255)
                       gAAGGGCTGGC.GG.TTG.GGGTTATTCGCAACGGCGAc              <  1:1306606/36‑1 (MQ=255)
                        aaGGGCTGGC.GG.TTT.GGGTTTTTTGCAACGGCGagg             >  1:7032813/1‑34 (MQ=18)
                        aaGGGCTGGC.GG.TTG.GGGTTATTCGCAACGGCGACt             <  1:3960355/36‑1 (MQ=255)
                         aGGGCTGGC.GG.TTG.GGGTTATTCGCAACGGCGACTg            <  1:7769533/36‑1 (MQ=255)
                          gggCTGGC.GG.TTT.GGGTTATTTTCCAAGGGGACTgg           >  1:7207363/1‑36 (MQ=9)
                           ggCTGGC.GG.TTG.GGGTTATTCGCAACGGCGACTGGc          <  1:5913683/36‑1 (MQ=255)
                            gCTGGC.GG.TTG.GGGTTATTCGCAACGGCGACTGGCt         <  1:7017898/36‑1 (MQ=255)
                              tGGC.GG.TTG.GGGTTATTCGCAACGGCGACTGGCTGt       <  1:1123886/36‑1 (MQ=255)
                              tGGC.GG.TTG.GGGTTATTCGCAACGGCGACTGGCTGt       <  1:6610209/36‑1 (MQ=255)
                                gC.GG.TTG.GGGTTATTCGCAACGGCGACTGGCTGTaa     <  1:4771282/36‑1 (MQ=255)
                                 c.GG.TTG.GGGTTATTCGCAACGGCGACTGGCTGTAAc    <  1:6382775/36‑1 (MQ=255)
                                   gg.TTG.GGGTTATTCGCAACGGCGACTGGCTGTAACa   <  1:1397283/36‑1 (MQ=255)
                                   gg.TTG.GGGTTATTCGCAACGGCGACTGGCTGTAACa   <  1:7693415/36‑1 (MQ=255)
                                   gg.TTG.GGGTTATTCGAAACGGCGACTGGATGTaaga   >  1:4069596/1‑34 (MQ=25)
                                    g.TTG.GGGTTATTCGCAACGGCGACTGGCTGTAACat  <  1:2772350/36‑1 (MQ=255)
                                     |                                    
TTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGC.GG.TTG.GGGTTATTCGCAACGGCGACTGGCTGTAACAT  >  REL606/4236588‑4236658

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 18 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG