diff --git a/workflow/scripts/runChipseeker.R b/workflow/scripts/runChipseeker.R
index e238a4aa4d7af33c9455abded45e0d4fb9c1e10c..b20df7833b2ee7bcd65dec1908712f941bf8479a 100644
--- a/workflow/scripts/runChipseeker.R
+++ b/workflow/scripts/runChipseeker.R
@@ -18,8 +18,8 @@ txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene
 }
 if(genome=="GRCm38")
 { 
-library(TxDb.Hsapiens.UCSC.mm10.knownGene)
-txdb <- TxDb.Hsapiens.UCSC.mm10.knownGene
+library(TxDb.Mmusculus.UCSC.mm10.knownGene)
+txdb <- TxDb.Mmusculus.UCSC.mm10.knownGene
 }
 if(genome=="GRCh38")
 { 
diff --git a/workflow/scripts/runMemechip.py b/workflow/scripts/runMemechip.py
index 36f84c512e79c1f6381bc2066954a849c4fc6b0d..b1f848f0376a449329594d1fe3c67b3b0e454c84 100644
--- a/workflow/scripts/runMemechip.py
+++ b/workflow/scripts/runMemechip.py
@@ -29,7 +29,7 @@ def prepare_argparser():
   argparser.add_argument("-l","--limit",dest = "limit",type=int,default=-1, help="Top limit of peaks")
   return(argparser)
 
-def rc():
+def rc(seq):
   comps = {'A':"T",'C':"G",'G':"C",'T':"A","N":"N"}
   return ''.join([comps[x] for x in seq.upper()[::-1]])