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

update picard samtools and bwa for DSL2

parent 0f8b9d13
No related branches found
No related tags found
No related merge requests found
version = "0.7.17"
container_url = "lbmc/bwa:${version}"
process index_fasta {
container = "${container_url}"
label "big_mem_mono_cpus"
tag "$file_id"
input:
tuple val(file_id), path(fasta)
output:
tuple val(fasta.simpleName), path("${fasta.simpleName}.*"), emit: index
file "*_bwa_report.txt", emit: report
script:
"""
bwa index -p ${fasta.simpleName} ${fasta} \
&> ${fasta.simpleName}_bwa_report.txt
"""
}
process mapping_fastq {
container = "${container_url}"
label "big_mem_multi_cpus"
tag "$file_id"
input:
tuple val(file_id), path(reads)
tuple val(index_id), path(index)
output:
tuple val(file_id), path("*.bam"), emit: bam
tuple val(file_id), path("${pair_id}_bwa_report.txt"), emit: report
script:
"""
bwa mem -t ${task.cpus} \
${index_id} ${reads[0]} ${reads[1]} 2> \
${file_id}_bwa_report.txt | \
sambamba view -S -t ${task/cpus} -f bam - > ${file_id}.bam
"""
}
......@@ -24,3 +24,21 @@ PicardCommandLine MarkDuplicates \
picard_${bam.baseName}.log
"""
}
process index_fasta {
container = "${container_url}"
label "big_mem_mono_cpus"
tag "$file_id"
input:
tuple val(file_id), file(fasta)
output:
tuple val(file_id), file("*.dict"), emit: index
script:
"""
PicardCommandLine CreateSequenceDictionary \
REFERENCE=${fasta} \
OUTPUT=${fasta.simpleName}.dict
"""
}
version = "1.11"
container_url = "lbmc/samtools:${version}"
process index_fasta {
container = "${container_url}"
label "big_mem_mono_cpus"
tag "$file_id"
input:
tuple val(file_id), path(fasta)
output:
tuple val(file_id), path("*.fai"), emit: index
script:
"""
samtools faidx ${fasta}
"""
}
process filter_bam {
container = "${container_url}"
label "big_mem_multi_cpus"
......
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