Skip to content

Commit

Permalink
Remove legacy rseqc strandRule code.
Browse files Browse the repository at this point in the history
  • Loading branch information
ewels committed Nov 21, 2018
1 parent 90b42ca commit 67be8b6
Showing 1 changed file with 1 addition and 7 deletions.
8 changes: 1 addition & 7 deletions main.nf
Original file line number Diff line number Diff line change
Expand Up @@ -554,7 +554,7 @@ if(params.aligner == 'star'){
script:
prefix = reads[0].toString() - ~/(_R1)?(_trimmed)?(_val_1)?(\.fq)?(\.fastq)?(\.gz)?$/
def star_mem = task.memory ?: params.star_memory ?: false
def avail_mem = star_mem ? "--limitBAMsortRAM ${star_mem.toBytes() - 100000000}" : ''
def avail_mem = star_mem ? "--limitBAMsortRAM ${star_mem.toBytes() - 100000000}" : ''
seqCenter = params.seqCenter ? "--outSAMattrRGline ID:$prefix 'CN:$params.seqCenter'" : ''
"""
STAR --genomeDir $index \\
Expand Down Expand Up @@ -712,12 +712,6 @@ process rseqc {
file "*.{txt,pdf,r,xls}" into rseqc_results

script:
def strandRule = ''
if (forward_stranded && !unstranded){
strandRule = params.singleEnd ? '-d ++,--' : '-d 1++,1--,2+-,2-+'
} else if (reverse_stranded && !unstranded){
strandRule = params.singleEnd ? '-d +-,-+' : '-d 1+-,1-+,2++,2--'
}
"""
samtools index $bam_rseqc
infer_experiment.py -i $bam_rseqc -r $bed12 > ${bam_rseqc.baseName}.infer_experiment.txt
Expand Down

0 comments on commit 67be8b6

Please sign in to comment.