diff --git a/workflow/conf/biohpc.config b/workflow/conf/biohpc.config index 5842b58b2e50dbb43b9a11fb67645391b0e8d557..732ddb895dd3faba09af8190719a340d33bff535 100755 --- a/workflow/conf/biohpc.config +++ b/workflow/conf/biohpc.config @@ -5,13 +5,17 @@ process { // Process specific configuration withName:splitData { - container = 'docker://bicf/bdbag:1.0' + container = 'docker://bicf/gudmaprbkfilexfer:1.1' } withName:getData { - container = 'docker://bicf/bdbag:1.0' + container = 'docker://bicf/gudmaprbkfilexfer:1.1' } withName:trimData { - container = 'docker://bicf/trimgalore:1.0' + container = 'docker://bicf/trimgalore:1.1' + queue = '256GB,256GBv1,384GB' + } + withName:alignReads { + container = 'docker://bicf/hisat2:1.0' queue = '256GB,256GBv1,384GB' } } diff --git a/workflow/rna-seq.nf b/workflow/rna-seq.nf index 035faa8ee8b3f2b95f298e4edc8f074a0b695587..df727d59d752230dedad70650e4c3dcc8de77318 100755 --- a/workflow/rna-seq.nf +++ b/workflow/rna-seq.nf @@ -3,6 +3,8 @@ // Define input variables params.deriva = "/project/BICF/BICF_Core/shared/gudmap/cookies/deriva-cookies.txt" params.bdbag = "${baseDir}/../test_data/Study_Q-Y4H0.zip" +params.spikein = "false" +params.species = "human" params.outDir = "${baseDir}/../output" @@ -92,7 +94,7 @@ process trimData { set repID, reads from trimming output: - path ("*_val_{1,2}.fq.gz", type: 'file', maxDepth: '0') + set repID, path ("*_val_{1,2}.fq.gz", type: 'file', maxDepth: '0') into aligning script: """ @@ -100,3 +102,15 @@ process trimData { trim_galore --gzip --max_n 1 --paired --basename \${rep} -j `nproc` ${reads[0]} ${reads[1]} 1>>\${rep}.trimData.log 2>>\${rep}.trimData.err; """ } + +/* + * alignReads: aligns the reads to a reference database +*/ +process alignReads { + tag "align-${repID}" + publishDir "${outDir}/tempOut/aligned", mode: "symlink" + + input: + set repID, fqs from aligning + +