add optional fasta argument to ensemblvep

This commit is contained in:
Nicolas Vannieuwkerke 2022-06-09 11:45:27 +02:00
parent 438f6c5281
commit 5ccf6fbcc9
4 changed files with 14 additions and 3 deletions

View file

@ -13,6 +13,7 @@ process ENSEMBLVEP {
val species
val cache_version
path cache
path fasta
path extra_files
output:
@ -27,6 +28,8 @@ process ENSEMBLVEP {
def args = task.ext.args ?: ''
def prefix = task.ext.prefix ?: "${meta.id}"
def dir_cache = cache ? "\${PWD}/${cache}" : "/.vep"
def reference = fasta ? "--fasta $fasta" : ""
"""
mkdir $prefix
@ -34,6 +37,7 @@ process ENSEMBLVEP {
-i $vcf \\
-o ${prefix}.ann.vcf \\
$args \\
$reference \\
--assembly $genome \\
--species $species \\
--cache \\

View file

@ -36,6 +36,11 @@ input:
type: file
description: |
path to VEP cache (optional)
- fasta:
type: file
description: |
reference FASTA file (optional)
pattern: "*.{fasta,fa}"
- extra_files:
type: tuple
description: |

View file

@ -10,5 +10,7 @@ workflow test_ensemblvep {
file(params.test_data['sarscov2']['illumina']['test_vcf'], checkIfExists: true)
]
ENSEMBLVEP ( input, "WBcel235", "caenorhabditis_elegans", "104", [], [] )
fasta = file(params.test_data['sarscov2']['genome']['genome_fasta'], checkIfExists: true)
ENSEMBLVEP ( input, "WBcel235", "caenorhabditis_elegans", "104", [], fasta, [] )
}

View file

@ -1,7 +1,7 @@
- name: ensemblvep test_ensemblvep
command: nextflow run ./tests/modules/ensemblvep -entry test_ensemblvep -c ./tests/config/nextflow.config -c ./tests/modules/ensemblvep/nextflow.config
command: nextflow run ./tests/modules/ensemblvep -entry test_ensemblvep -c ./tests/config/nextflow.config -c ./tests/modules/ensemblvep/nextflow.config
tags:
- ensemblvep
files:
- path: output/ensemblvep/test.ann.vcf
- path: output/ensemblvep/test.summary.html
- path: output/ensemblvep/test.summary.html