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

g2gtools: add incorporate_snp and incorporate_indel processes

parent 308f6c25
No related branches found
No related tags found
No related merge requests found
......@@ -8,9 +8,9 @@ process vci_build {
input:
tuple val(file_id), path(vcf)
tuple val(ref_id), path(fasta), path(fai), path(dict)
tuple val(ref_id), path(fasta)
output:
tuple val(file_id), path("*.chain"), emit: vci
tuple val(file_id), path("*.vci.gz"), emit: vci
script:
"""
g2gtools vcf2vci \
......@@ -18,6 +18,45 @@ g2gtools vcf2vci \
-f ${fasta} \
-i ${vcf} \
-s ${file_id} \
-o ${file_id}.vci
-o ${file_id}.vci.gz
"""
}
process incorporate_snp {
container = "${container_url}"
label "big_mem_multi_cpus"
tag "$file_id"
input:
tuple val(file_id), path(vci)
tuple val(ref_id), path(fasta)
output:
tuple val(file_id), path("${file_id}_snp.fasta"), path("${vci}"), emit: fasta
script:
"""
g2gtools patch \
-p ${task.cpus} \
-i ${fasta} \
-c ${vci} \
-o ${file_id}_snp.fasta
"""
}
process incorporate_indel {
container = "${container_url}"
label "big_mem_multi_cpus"
tag "$file_id"
input:
tuple val(file_id), path(fasta), path(vci)
output:
tuple val(file_id), path("${file_id}_snp_indel.fasta"), path("${vci}"), emit: fasta
script:
"""
g2gtools transform \
-p ${task.cpus} \
-i ${fasta} \
-c ${vci} \
-o ${file_id}_snp_indel.fasta
"""
}
\ No newline at end of file
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