fix chipseeker mm10

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")
 { 

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]])