Skip to content

Commit

Permalink
some minor changes to the code
Browse files Browse the repository at this point in the history
  • Loading branch information
lucacozzuto committed Aug 10, 2018
1 parent d38c0cf commit 7ea70bf
Show file tree
Hide file tree
Showing 4 changed files with 32 additions and 30 deletions.
2 changes: 1 addition & 1 deletion lib/FeatureCounter.groovy
Original file line number Diff line number Diff line change
Expand Up @@ -26,7 +26,7 @@
def public htseqCount() {

"""
htseq-count ${this.extrapars} -s ${this.strand} -f bam ${this.input} ${this.annotation} > ${this.output}
htseq-count ${this.extrapars} -s ${this.strand} -f bam ${this.input} ${this.annotation} ${this.extrapars} > ${this.output}
"""
}

Expand Down
26 changes: 14 additions & 12 deletions lib/Misc.groovy
Original file line number Diff line number Diff line change
Expand Up @@ -100,11 +100,12 @@
def public makeAlnProfilesWithDeepTools() {

"""
bamCoverage --bam ${this.input} -o ${this.output} \\
--binSize 10 \\
--normalizeUsing RPGC \\
--effectiveGenomeSize ${this.genome_size} \\
--extendReads ${this.read_size} \\
bamCoverage --bam ${this.input} -o ${this.output} \
--binSize 10 \
--normalizeUsing RPGC \
--effectiveGenomeSize ${this.genome_size} \
--extendReads ${this.read_size} \
${this.extrapars} \
-p ${this.cpus}
"""
}
Expand Down Expand Up @@ -136,13 +137,14 @@
def public markDuplicateWithPicard() {

"""
java -jar ${this.java_path} MarkDuplicates \\
INPUT=${this.input} \\
OUTPUT=${this.output} \\
ASSUME_SORTED=true \\
REMOVE_DUPLICATES=true \\
METRICS_FILE=${this.id}.metrics.txt \\
VALIDATION_STRINGENCY=LENIENT
java -jar ${this.java_path} MarkDuplicates \
INPUT=${this.input} \
OUTPUT=${this.output} \
ASSUME_SORTED=true \
REMOVE_DUPLICATES=true \
METRICS_FILE=${this.id}.metrics.txt \
VALIDATION_STRINGENCY=LENIENT \
${this.extrapars}
"""
}

Expand Down
26 changes: 13 additions & 13 deletions lib/PeakCaller.groovy
Original file line number Diff line number Diff line change
Expand Up @@ -45,12 +45,12 @@

"""
export PYTHONPATH=""
macs2 callpeak \\
-t ${this.sample} \\
-f BAM \\
-n ${this.id} \\
--extsize ${this.frag_len} \\
${extrapars} \\
macs2 callpeak \
-t ${this.sample} \
-f BAM \
-n ${this.id} \
--extsize ${this.frag_len} \
${extrapars} \
-g ${this.genome_size} 2>${this.id}_log.txt
"""
}
Expand All @@ -61,13 +61,13 @@
def private peakCallWithEpic() {

"""
epic --treatment ${this.sample} \\
--control ${this.control} \\
--chromsizes ${this.chr_sizes} \\
--effective-genome-fraction ${this.genome_size} \\
-o ${this.id}.out \\
-b ${this.id}.bed \\
--fragment-size ${this.frag_len} \\
epic --treatment ${this.sample} \
--control ${this.control} \
--chromsizes ${this.chr_sizes} \
--effective-genome-fraction ${this.genome_size} \
-o ${this.id}.out \
-b ${this.id}.bed \
--fragment-size ${this.frag_len} \
${extrapars}
"""
}
Expand Down
8 changes: 4 additions & 4 deletions lib/QualityChecker.groovy
Original file line number Diff line number Diff line change
Expand Up @@ -116,7 +116,7 @@
-T "Fingerprints" \
--plotFile ${this.output}.png \
--outRawCounts ${this.output}_rawcounts.tab \
--numberOfProcessors ${this.cpus} \
--numberOfProcessors ${this.cpus} ${this.extrapars} \
"""
}

Expand All @@ -133,7 +133,7 @@
-b ${this.input} \
--smartLabels \
--outFileName ${this.output} \
--numberOfProcessors ${this.cpus} \
--numberOfProcessors ${this.cpus} ${this.extrapars} \
"""
}

Expand Down Expand Up @@ -186,7 +186,7 @@
--regionBodyLength 5000 \
--afterRegionStartLength 3000 \
--numberOfProcessors ${this.cpus} \
--skipZeros -o ${this.output}
--skipZeros -o ${this.output} ${this.extrapars}
"""
}

Expand All @@ -205,7 +205,7 @@
--beforeRegionStartLength 3000 \
--afterRegionStartLength 3000 \
--numberOfProcessors ${this.cpus} \
--skipZeros -o ${this.output}
--skipZeros -o ${this.output} ${this.extrapars}
"""
}

Expand Down

0 comments on commit 7ea70bf

Please sign in to comment.