diff --git a/workflow/rna-seq.nf b/workflow/rna-seq.nf index 454f6e7d1de1bf96bbb8fb548cae8966e0c63eb2..6dde60106d5c085917e61c56f131e6370f227842 100755 --- a/workflow/rna-seq.nf +++ b/workflow/rna-seq.nf @@ -64,11 +64,12 @@ process getData { input: val repRID_getData + executor 'local' path cookies, stageAs: 'deriva-cookies.txt' from bdbag path bagit output: - path ("**/*.R{1,2}.fastq.gz") into fastqs + path ("*.R{1,2}.fastq.gz") into fastqs script: """ @@ -81,7 +82,7 @@ process getData { echo "LOG: \${replicate}" >>${repRID_getData}.getData.err unzip ${bagit} 2>>${repRID_getData}.getData.err echo "LOG: replicate bdbag unzipped" >>${repRID_getData}.getData.err - sh ${baseDir}/scripts/bdbagFetch.sh \${replicate} 2>>${repRID_getData}.getData.err + sh ${baseDir}/scripts/bdbagFetch.sh \${replicate} ${repRID_getData} 2>>${repRID_getData}.getData.err echo "LOG: replicate bdbag fetched" >>${repRID_getData}.getData.err """ } @@ -91,7 +92,6 @@ process getData { */ process trimData { tag "${repRID_trimData}" - publishDir "${outDir}/tempOut/trimmed", mode: "symlink", pattern: "*_val_{1,2}.fq.gz" publishDir "${logsDir}/trimData", mode: 'symlink', pattern: "\${repRID_trimData}.trimData.*" input: @@ -99,7 +99,8 @@ process trimData { file(fastq) from fastqs output: - path ("*_val_{1,2}.fq.gz", type: 'file', maxDepth: '0') + path ("*.fq.gz") into fastqs_trimmed + val ends script: """ @@ -109,10 +110,12 @@ process trimData { else ncore=`nproc` fi - if [ -z ${fastq[1]} ] + if [ '${fastq[1]}' == 'null' ] then + ends='se' trim_galore --gzip -q 25 --illumina --length 35 --basename ${repRID_trimData} -j \${ncore} ${fastq[0]} 1>>${repRID_trimData}.trimData.log 2>>${repRID_trimData}.trimData.err; else + ends='pe' trim_galore --gzip -q 25 --illumina --length 35 --paired --basename ${repRID_trimData} -j \${ncore} ${fastq[0]} ${fastq[1]} 1>>${repRID_trimData}.trimData.log 2>>${repRID_trimData}.trimData.err; fi """ diff --git a/workflow/scripts/bdbagFetch.sh b/workflow/scripts/bdbagFetch.sh index 65d309d64e28cd2767b02e8080b8d7fb320f5f39..902222a2ebb6aa7e978f0a820ad3c04472395848 100644 --- a/workflow/scripts/bdbagFetch.sh +++ b/workflow/scripts/bdbagFetch.sh @@ -1,3 +1,7 @@ #!/bin/bash -bdbag --resolve-fetch all --fetch-filter filename\$*fastq.gz $1 \ No newline at end of file +bdbag --resolve-fetch all --fetch-filter filename\$*fastq.gz $1 +for i in $(find */ -name "*.R*.fastq.gz"); do + path=${2}$(echo ${i##*/} | grep -o "\.R.\.fastq\.gz"); + mv ${i} ./${path} +done; \ No newline at end of file