我使用vcf2maf
作为snakemake
管道的一部分来注释变体
rule vcf2maf:
input:
vcf="vcfs/{sample}.vcf",
fasta=vep_fasta,
vep_dir=vep_dir
output:
"mafs/{sample}.maf"
conda:
"../envs/annotation.yml"
shell:
"""
vcf2maf.pl --input-vcf {input.vcf} --output-maf {output} \
--tumor-id {wildcards.sample}.tumor \
--normal-id {wildcards.sample}.normal \
--ref-fasta {input.fasta} --filter-vcf 0 \
--vep-data {input.vep_dir} --vep-path [need path]
"""
conda
环境有两个包:vcf2maf
和vep
vcf2maf
需要到vep
的路径才能正常运行,但我不确定如何访问vep
的路径,因为它存储在conda
环境中,该环境将具有用户特定的绝对路径。有没有一种简单的方法来获取vep
的路径,以便我可以参考它来查找--vep-path
您可以使用unix
which
命令,如下所示:变量
CONDA_PREFIX
包含当前conda环境的路径。所以你也可以这样做:相关问题 更多 >
编程相关推荐