Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 874,932 +G 7.7% coding (355/759 nt) deoR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606874,9321.→G 7.7% 3.7 ‑0.7 ‑0.426coding (355/759 nt)deoRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (14/10):  total (16/10)
Fisher's exact test strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.97e-01

CAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGC.GGTAAAAGGGATTTCATTATCAATCGCTTCAATAA  >  REL606/874897‑874967
                                    |                                   
cAGCGCCAGAAAGGTATTTAGCGAATAACAAAcggc                                      <  1:2307710/36‑1 (MQ=255)
 aGCGCCAGAAAGGTATTTAGCGAATAACAAAcggc.g                                    >  1:5550710/1‑36 (MQ=255)
 aGCGCCAGAAAGGTATTTAGCGAATAACAAAcggc.g                                    <  1:1612225/36‑1 (MQ=255)
   cgcCAGAAAGGTATTTAGCGAATAACAAACGGC.GGt                                  >  1:5405841/1‑36 (MQ=255)
    gcCAGAAAGGTATTTAGCGAATAACAAAcggc.ggga                                 >  1:1808727/1‑34 (MQ=255)
    gcCAGAAAGGTATTTAGCGAATAACAAACGGC.GGTa                                 <  1:5673198/36‑1 (MQ=255)
    gcCAGAAAGGTATTTAGCGAATAACAAACGGC.GGTa                                 >  1:2298649/1‑36 (MQ=255)
         aaaGGTATTTAGCGAATAACAAACGGC.GGTgaaaag                            >  1:1713794/1‑34 (MQ=37)
          aaGGTATTTAGCGAATAACAAACGGCGGGTAAAAgg                            >  1:23310/1‑36 (MQ=37)
            ggTATTTAGCGAATAACAAACGGC.GGGAAAAAGGag                         >  1:2968639/1‑35 (MQ=25)
             gTATTTAGCGAATAACAAACGGC.GGAAAAAaggggt                        >  1:2305770/1‑33 (MQ=25)
               aTTTAGCGAATAACAAACGGCGGGCAAAAGGGTttt                       >  1:2512062/1‑36 (MQ=19)
               aTTTAGCGAATAACAAACGGC.GGTAAAAGGGATTTc                      <  1:3922690/36‑1 (MQ=255)
               aTTTAGCGAATAACAAACGGC.GGTAAAAAGGGTTTc                      >  1:2509085/1‑36 (MQ=37)
               aTTTAGCGAATAACAAACGGC.GGGAAAAAGGGttta                      >  1:862602/1‑35 (MQ=21)
                 ttAGCGAATAACAAACGGC.GGTAAAAAGGATTTTAt                    >  1:4732601/1‑36 (MQ=37)
                   aGCGAATAACAAACGGC.GGGAAAAAGGGTTTCCTTa                  >  1:4455960/1‑36 (MQ=16)
                      gAATAACAAACGGC.GGCAAAAAGGGTTTTACTAttt               >  1:3284689/1‑34 (MQ=12)
                       aaTAACAAACGGC.GGGAAAAAGGGTTTTACTATCaa              >  1:789944/1‑36 (MQ=12)
                       aaTAACAAACGGC.GGAAAAAAGGATTTCGTTATCaa              >  1:656975/1‑36 (MQ=25)
                        aTAACAAACGGC.GGTAAAAGGGATTTCATTATCAAt             <  1:1032470/36‑1 (MQ=255)
                        aTAACAAACGGC.GGTAAAAGGGATTTCATTATCAAt             <  1:5888580/36‑1 (MQ=255)
                        aTAACAAACGGC.GGTAAAAGGGATTTCATTATCAAt             <  1:4545974/36‑1 (MQ=255)
                         tAACAAACGGC.GGTAAAAGGGATTTCATTATCAATc            <  1:4982687/36‑1 (MQ=255)
                         tAACAAACGGC.GGTAAAAGGGATTTCATTATCAATc            <  1:1862772/36‑1 (MQ=255)
                         tAACAAACGGC.GGTAAAAGGGATTTCATTATCAATc            <  1:3319823/36‑1 (MQ=255)
                          aaCAAACGGC.GGTAAAAAGGGTTTTATTATCaaaac           >  1:2699963/1‑33 (MQ=16)
                          aaCAAACGGC.GGCAAAAGGGATCTCATTATCaaacc           >  1:346929/1‑33 (MQ=25)
                               aCGGC.GGTAAAAGGGATTTCTTTATCAATCGCTTCa      <  1:3273037/36‑1 (MQ=255)
                               aCGGC.GGTAAAAGGGATTTCATTATCAATCGCTTCa      <  1:2725810/36‑1 (MQ=255)
                                 ggc.ggTAAAAGGGATTTCATTATCAATCGCTTCaat    <  1:1492507/36‑1 (MQ=255)
                                 ggc.ggTAAAAGGGATTTCATTATCAATCGCTTCaat    <  1:5856816/36‑1 (MQ=255)
                                  gc.ggTAAAAGGTATTTCATTATCAATCGCTTCaata   <  1:4186121/36‑1 (MQ=255)
                                  gc.ggTAAAAGGGGTTTCATTTTCAATCGCTTCaata   >  1:12978/1‑36 (MQ=25)
                                   c.ggTAAAAGGGATTTCATTATCAATCGCTTCaataa  <  1:5822279/36‑1 (MQ=255)
                                    |                                   
CAGCGCCAGAAAGGTATTTAGCGAATAACAAACGGC.GGTAAAAGGGATTTCATTATCAATCGCTTCAATAA  >  REL606/874897‑874967

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG