Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,135,984 A→G 31.4% S26P (TCC→CCC)  dinI ← DNA damage‑inducible protein I

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,135,9840A→G31.4% 3.1 ‑1.9 0.01876dinIDNA damage‑inducible protein I
Reads supporting (aligned to +/- strand):  new base (5/1):  ref base (5/7):  total (10/8)
Fisher's exact test strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGC  >  REL606/1135949‑1136017
                                   |                                 
ttCATTATCAGGAAACGCATACTGAATACGGCGGGa                                   >  1:4318198/1‑36 (MQ=255)
ttCATTATCAGGAAACGCATACTGAATACGGCGGGa                                   <  1:4189701/36‑1 (MQ=255)
ttCATTAACAGGAAACGCATACTGAATACGGCGGGa                                   <  1:4017648/36‑1 (MQ=255)
   aTTATCAGGAAACGCATACTGAATACGGCGGGgaaa                                >  1:2756948/1‑35 (MQ=255)
     tATCAGGAAACGCATGCTGAATACGGCGGGGAAGtt                              <  1:1332604/36‑1 (MQ=37)
      aTCAGGAAACGCATACTGAATACGGCGGGGAAggcg                             >  1:1156278/1‑33 (MQ=37)
      aTCAGGAAACGCATACTGAATACGGAGGGGAAggtt                             >  1:4014784/1‑33 (MQ=25)
       tCAGGAAACGCATACTGAATACGGCGGGAAAGTTCg                            <  1:2281957/36‑1 (MQ=255)
          ggAAACGCATACTGAATACGGCGGGAAAGTTCgcca                         <  1:2155271/36‑1 (MQ=255)
           gAAACGCATACTGAATACGGCGGGGAAGTTTgcccg                        >  1:5176359/1‑34 (MQ=37)
           gAAACGCATACTGAATACGGCGGGGAAGTTCgcccg                        >  1:124421/1‑34 (MQ=38)
            aaaCGCATACTGAATACGGCGGGGAAGTTCtccccg                       >  1:3930479/1‑33 (MQ=37)
              aCGCATACTGAATACGGCGGTAAAGTTCtccagcca                     <  1:3203205/36‑1 (MQ=37)
                     cTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTc              <  1:972768/36‑1 (MQ=255)
                     cTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTc              <  1:5325399/36‑1 (MQ=255)
                       gAATACGGCGGAAAAGTTCGCCAGCCAGGGCGTCaa            <  1:2092367/36‑1 (MQ=255)
                         aTACGGCGGGAAAGTTCGCCAGCCCGGGCGGCAATa          >  1:2731479/1‑36 (MQ=37)
                             ggcggGAAAGTTCGCCAGCCAGGGCGTCAATAGccc      >  1:5427779/1‑36 (MQ=255)
                              gcggGAAAGTTCGCCAGCCAGGGCGTCAATAGcccc     >  1:2719231/1‑36 (MQ=255)
                               cggAAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCa    <  1:3076804/36‑1 (MQ=255)
                                gggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAg   >  1:5667860/1‑36 (MQ=255)
                                 ggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCCGc  >  1:1011319/1‑36 (MQ=255)
                                   |                                 
TTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGC  >  REL606/1135949‑1136017

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 15 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG