|
| 1 | +process SortMeRNA { |
| 2 | + tag {"SortMeRNA ${sample_id}"} |
| 3 | + label 'SortMeRNA_4_3_3' |
| 4 | + container = 'quay.io/biocontainers/sortmerna:4.3.3--h9ee0642_0' |
| 5 | + shell = ['/bin/bash', '-euo', 'pipefail'] |
| 6 | + |
| 7 | + input: |
| 8 | + tuple(sample_id, rg_id, path(fastq_files)) |
| 9 | + path(db_fasta) |
| 10 | + |
| 11 | + output: |
| 12 | + tuple(sample_id, rg_id, path("*_non_rRNA.fastq.gz"), emit: non_rRNA_fastqs) |
| 13 | + path("*_filtered_rRNA.fastq.gz", emit: rRNA_fastqs) |
| 14 | + path("*_rRNA_report.txt", emit: qc_report) |
| 15 | + |
| 16 | + script: |
| 17 | + def refs = db_fasta.collect{ "$it" }.join(" -ref ") |
| 18 | + def report_title = fastq_files[0].simpleName.split("_R1_")[0] |
| 19 | + if (params.single_end) { |
| 20 | + """ |
| 21 | + sortmerna -ref ${refs} \ |
| 22 | + -reads ${fastq_files} \ |
| 23 | + --num_alignments 1 \ |
| 24 | + --threads ${task.cpus} \ |
| 25 | + --fastx \ |
| 26 | + -workdir \${PWD} \ |
| 27 | + --aligned rRNA-reads \ |
| 28 | + --other non-rRNA-reads \ |
| 29 | + --zip-out |
| 30 | + |
| 31 | + mv non-rRNA-reads.fq.gz ${fastq_files[0].simpleName}_non_rRNA.fastq.gz |
| 32 | + mv rRNA-reads.fq.gz ${fastq_files[0].simpleName}_filtered_rRNA.fastq.gz |
| 33 | + mv rRNA-reads.log ${report_title}_rRNA_report.txt |
| 34 | + """ |
| 35 | + } else { |
| 36 | + """ |
| 37 | + sortmerna -ref ${refs} \ |
| 38 | + -reads ${fastq_files[0]} -reads ${fastq_files[1]} \ |
| 39 | + --num_alignments 1 \ |
| 40 | + --threads ${task.cpus} \ |
| 41 | + -workdir \${PWD} \ |
| 42 | + --fastx -paired_in \ |
| 43 | + --aligned rRNA-reads \ |
| 44 | + --other non-rRNA-reads \ |
| 45 | + -out2 \ |
| 46 | + --zip-out |
| 47 | +
|
| 48 | + mv non-rRNA-reads_fwd.fq.gz ${fastq_files[0].simpleName}_non_rRNA.fastq.gz |
| 49 | + mv non-rRNA-reads_rev.fq.gz ${fastq_files[1].simpleName}_non_rRNA.fastq.gz |
| 50 | + mv rRNA-reads_fwd.fq.gz ${fastq_files[0].simpleName}_filtered_rRNA.fastq.gz |
| 51 | + mv rRNA-reads_rev.fq.gz ${fastq_files[1].simpleName}_filtered_rRNA.fastq.gz |
| 52 | + mv rRNA-reads.log ${report_title}_rRNA_report.txt |
| 53 | + """ |
| 54 | + } |
| 55 | +} |
0 commit comments