Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 340,545 C→T 26.4% L182L (CTG→CTA mhpR ← DNA‑binding transcriptional activator, 3HPP‑binding

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606340,5450C→T26.4% 2.7 ‑3.7 0.015546mhpRDNA‑binding transcriptional activator, 3HPP‑binding
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (3/8):  total (7/8)
Fisher's exact test strand distribution p-value = 2.56e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCAAACGTCGCCCGACCATTGCCCGGTG  >  REL606/340514‑340580
                               |                                   
aGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGaa                                 <  1:2481414/36‑1 (MQ=255)
  gccagccaGGTCAGGCCCGAGGCTGTTTTCAGAAGc                               <  1:78257/36‑1 (MQ=255)
  gccagccaGGTCAGGCCCGAGGCGGTTTTTAGGAgg                               >  1:4538872/1‑35 (MQ=37)
     agcctgGTCAGGCCCGAGGCGGTTTTCAGAAGCGGc                            <  1:428650/36‑1 (MQ=255)
     agccagGTCAGGCCCGAGGCGGTTTTTTGAAGAggg                            >  1:1831883/1‑35 (MQ=25)
     agccagGTCAGGCCCGAGGCGGTTTTCAGAAGCGGc                            <  1:571791/36‑1 (MQ=255)
      gccagGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCa                           <  1:3337266/36‑1 (MQ=255)
         agGTCAGGCCCGAGGCGGTTTTTAGAAGCGGTaaaa                        >  1:1603267/1‑35 (MQ=25)
             cAGGCCCGAGGCGGGTTTTAAAAGCGGCTAACGGCg                    >  1:5087412/1‑36 (MQ=16)
              aGGCCCGAGGCGGTTTTCAGAAGCGGCAAACGTCGc                   <  1:3655526/36‑1 (MQ=255)
                     aGGCGGTTTTCAGAAGCGGCAAACGTCGCCCGACCa            >  1:5707178/1‑36 (MQ=255)
                       gCGGTTTTCAGAAGAGGCAAACGTCCCCCGACCAtt          >  1:3990773/1‑36 (MQ=38)
                           ttttCAGAAGCGGCAAACGTCGCCCGACCATTGccc      <  1:5018939/36‑1 (MQ=255)
                             ttCAGAAGCGGCAAACGTCGCCCGACCATTGCCCgg    >  1:1527191/1‑36 (MQ=255)
                             ttCAGAAGCGGCAAACGTCGCCCGACCATTGCCCgg    <  1:4325736/36‑1 (MQ=255)
                               cAGAAGCGGCAAACGTCGCCCGACCATTGCCCgggg  >  1:713647/1‑34 (MQ=255)
                               |                                   
AGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCAAACGTCGCCCGACCATTGCCCGGTG  >  REL606/340514‑340580

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