Skip to content
Snippets Groups Projects
Unverified Commit 99a62a30 authored by Laurent Modolo's avatar Laurent Modolo
Browse files

create RNASeq.nf

parent e71dcd50
No related branches found
No related tags found
No related merge requests found
log.info "fastq files : ${params.fastq}"
Channel
.fromFilePairs( params.fastq )
.ifEmpty { error "Cannot find any fastq files matching: ${params.fastq}" }
.set { fastq_files }
process adaptor_removal {
tag "$pair_id"
publishDir "results/fastq/adaptor_removal/", mode: 'copy'
input:
set pair_id, file(reads) from fastq_files
output:
set pair_id, "*_cut_R{1,2}.fastq.gz" into fastq_files_cut
script:
"""
cutadapt -a AGATCGGAAGAG -g CTCTTCCGATCT -A AGATCGGAAGAG -G CTCTTCCGATCT \
-o ${pair_id}_cut_R1.fastq.gz -p ${pair_id}_cut_R2.fastq.gz \
${reads[0]} ${reads[1]} > ${pair_id}_report.txt
"""
}
......@@ -9,7 +9,7 @@ process sample_fasta {
file fasta from fasta_file
output:
file "sample.fasta" into fasta_sample
file "*_sample.fasta" into fasta_sample
script:
"""
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment