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

bwa: add single_end mapping

parent c64d2dd0
No related branches found
No related tags found
No related merge requests found
...@@ -32,7 +32,7 @@ process mapping_fastq { ...@@ -32,7 +32,7 @@ process mapping_fastq {
output: output:
tuple val(file_id), path("*.bam"), emit: bam tuple val(file_id), path("*.bam"), emit: bam
tuple val(file_id), path("${file_id}_bwa_report.txt"), emit: report tuple val(file_id), path("${id}_bwa_report.txt"), emit: report
script: script:
if (file_id.containsKey('library')) { if (file_id.containsKey('library')) {
...@@ -43,6 +43,7 @@ if (file_id.containsKey('library')) { ...@@ -43,6 +43,7 @@ if (file_id.containsKey('library')) {
id = file_id id = file_id
} }
bwa_mem_R = "@RG\\tID:${library}\\tSM:${library}\\tLB:lib_${library}\\tPL:illumina" bwa_mem_R = "@RG\\tID:${library}\\tSM:${library}\\tLB:lib_${library}\\tPL:illumina"
if (reads instanceof List)
""" """
bwa mem -t ${task.cpus} \ bwa mem -t ${task.cpus} \
-R '${bwa_mem_R}' \ -R '${bwa_mem_R}' \
...@@ -50,5 +51,14 @@ ${index_id} ${reads[0]} ${reads[1]} 2> \ ...@@ -50,5 +51,14 @@ ${index_id} ${reads[0]} ${reads[1]} 2> \
${id}_bwa_report.txt | \ ${id}_bwa_report.txt | \
samtools view -@ ${task.cpus} -Sb - > ${id}.bam samtools view -@ ${task.cpus} -Sb - > ${id}.bam
""" """
else
"""
bwa mem -t ${task.cpus} \
-R '${bwa_mem_R}' \
${index_id} ${reads} 2> \
${id}_bwa_report.txt | \
samtools view -@ ${task.cpus} -Sb - > ${id}.bam
"""
} }
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment