0

我想使用snakemake一个接一个地运行多个规则。但是，当我运行这个脚本时，bam_list 规则出现在 samtools_markdup 规则之前，并给我一个错误，它找不到输入文件，这些文件显然还没有生成。如何解决这个问题呢？

rule all:
    input: 
        expand("dup/{sample}.dup.bam", sample=SAMPLES)
        "dup/bam_list"

rule samtools_markdup:
    input:
        sortbam ="rg/{sample}.rg.bam"
    output:
        dupbam = "dup/{sample}.dup.bam"
    threads: 5
    shell:
        """
        samtools markdup -@ {threads} {input.sortbam} {output.dupbam}
        """

rule bam_list:
    output:
         outlist = "dup/bam_list"
    shell:
         """
         ls dup/*.bam > {output.outlist}
         """
4

1 回答 1

1

Snakemake 遵循您想要的指示，dup/bam_list并且无需任何输入即可生成。我认为你的意思是：

rule all:
    input: 
        "dup/bam_list"

rule samtools_markdup:
    input:
        sortbam ="rg/{sample}.rg.bam"
    output:
        dupbam = "dup/{sample}.dup.bam"
    threads: 5
    shell:
        """
        samtools markdup -@ {threads} {input.sortbam} {output.dupbam}
        """

rule bam_list:
    input: 
        expand("dup/{sample}.dup.bam", sample=SAMPLES)
    output:
        outlist = "dup/bam_list"
    shell:
         """
         ls dup/*.bam > {output.outlist}
         """

现在 bam_list 将等到所有 samtools_markdup 作业完成。顺便说一句，我希望 dup_list 的内容与 相同expand("dup/{sample}.dup.bam", sample=SAMPLES)，因此如果您稍后在工作流中使用该文件，您可能只使用扩展输出。

于 2020-10-20T19:51:22.683 回答