硬编码输出,不扩展Snakefile

时间:2019-03-20 23:24:36

标签: output directory snakemake expansion

我的Snakefile如下:

SAMPLES, = glob_wildcards("data/{sample}_R1.fq.gz")
rule all: 
    input:
        expand("samtools_sorted_out/{sample}.raw.snps.indels.g.vcf", sample=SAMPLES), 
        expand("samtools_sorted_out/combined_gvcf")
rule combine_gvcf:
    input: "samtools_sorted_out/{sample}.raw.snps.indels.g.vcf"
    output:directory("samtools_sorted_out/combined_gvcf")
    params: gvcf_file_list="gvcf_files.list",
            gatk4="/storage/anaconda3/envs/exome/share/gatk4-4.1.0.0-0/gatk-package-4.1.0.0-local.jar"
    shell:"""
        java -DGATK_STACKTRACE_ON_USER_EXCEPTION=true \
            -jar {params.gatk4} GenomicsDBImport \
            -V {params.gvcf_file_list} \
            --genomicsdb-workspace-path {output}
       """

以空转测试时,出现错误:

RuleException in line 335 of /data/yifangt/exomecapture/Snakefile:
Wildcards in input, params, log or benchmark file of rule combine_gvcf cannot be determined from output files:
'sample'

我需要在两个地方提供帮助:

  1. {output}是将由shell部件创建的文件夹;
  2. {output}文件夹是命令行手动要求的硬编码(其内容事先未知)。

与{input}相比,问题似乎是{output}没有扩展。 我应该如何处理这种情况?非常感谢!

0 个答案:

没有答案