Skip to content
Snippets Groups Projects
Commit 22a3597f authored by Laurent Modolo's avatar Laurent Modolo
Browse files

RNAseq_sen1D.nf: fix bowtie2 index problem

parent b39575a3
No related branches found
No related tags found
No related merge requests found
...@@ -221,8 +221,14 @@ process mapping_fastq { ...@@ -221,8 +221,14 @@ process mapping_fastq {
file "*_bowtie2_report.txt" into mapping_fastq_report file "*_bowtie2_report.txt" into mapping_fastq_report
script: script:
index_id = index[0]
for (index_file in index) {
if (index_file =~ /.*\.1\.bt2/ && !(index_file =~ /.*\.rev\.1\.bt2/) ) {
index_id = ( index_file =~ /(.*)\.1\.bt2/)[0][1]
}
}
""" """
bowtie2 --very-sensitive -p ${task.cpus} -x ${index[0].baseName} \ bowtie2 --very-sensitive -p ${task.cpus} -x ${index_id} \
-1 ${reads[0]} -2 ${reads[1]} 2> \ -1 ${reads[0]} -2 ${reads[1]} 2> \
${pair_id}_bowtie2_report.txt | \ ${pair_id}_bowtie2_report.txt | \
samtools view -Sb - > ${pair_id}.bam samtools view -Sb - > ${pair_id}.bam
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment