Skip to content
Snippets Groups Projects
Commit ec41d5e0 authored by vvanoost's avatar vvanoost
Browse files

modify params.fastq to include the .gz format and modify the output to include...

modify params.fastq to include the .gz format and modify the output to include _001 and include the 5' adapter
parent 742e683a
No related branches found
No related tags found
No related merge requests found
...@@ -10,7 +10,7 @@ ...@@ -10,7 +10,7 @@
* for paired-end data * for paired-end data
*/ */
params.fastq = "$baseDir/data/fastq/*_R{1,2}_001.fastq" params.fastq = "$baseDir/data/fastq/*_R{1,2}_001.fastq.gz"
log.info "fastq files : ${params.fastq}" log.info "fastq files : ${params.fastq}"
...@@ -27,11 +27,11 @@ process adaptor_removal { ...@@ -27,11 +27,11 @@ process adaptor_removal {
set pair_id, file(reads) from fastq_files set pair_id, file(reads) from fastq_files
output: output:
set pair_id, "*_cut_R{1,2}.fastq.gz" into fastq_files_cut set pair_id, "*_cut_R{1,2}_001.fastq.gz" into fastq_files_cut
script: script:
""" """
cutadapt -a TGGAATTCTCGGGTGCCAAGG -g CCTTGGCACCCGAGAATTCCA -A TGGAATTCTCGGGTGCCAAGG -G CCTTGGCACCCGAGAATTCCA \ cutadapt -a TGGAATTCTCGGGTGCCAAGG -A GTTCAGAGTTCTACAGTCCGACGATC \
-o ${pair_id}_cut_R1.fastq.gz -p ${pair_id}_cut_R2.fastq.gz \ -o ${pair_id}_cut_R1.fastq.gz -p ${pair_id}_cut_R2.fastq.gz \
${reads[0]} ${reads[1]} > ${pair_id}_report.txt ${reads[0]} ${reads[1]} > ${pair_id}_report.txt
""" """
......
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