Question

我开始使用snakemake。我设法定义了一些我可以独立运行但不在工作流程中的规则。也许问题在于他们有无关的输入和输出。

我目前的工作流程是这样的：

configfile: './config.yaml'

rule all:
    input: dynamic("task/{job}/taskOutput.tab")
rule split_input:
     input: "input_fasta/snp.fa"
     output: dynamic("task/{job}/taskInput.fa")
     shell:
     "rm -Rf tasktmp task; \
     mkdir tasktmp task; \
     split -l 200 -d {input} ./tasktmp/; \
     ls tasktmp | awk '{{print \"mkdir task/\"$0}}' | sh; \
     ls tasktmp | awk '{{print \"mv ./tasktmp/\"$0\" ./task/\"$0\"/taskInput.fa\"}}' | sh"
rule task:
     input: "task/{job}/taskInput.fa"
     output: "task/{job}/taskOutput.tab"
     shell: "cp {input} {output}"
rule make_parameter_file:
     output:
    "par/parameters.txt
     shell:
    "rm -Rf par;mkdir par; \
    echo \"\
minimumFlankLength=5\n\
maximumFlankLength=200\n\
alignmentLengthDifference=2\
allowedMismatch=4\n\
allowedProxyMismatch=2\n\
allowedIndel=3\n\
ambiguitiesAsMatch=1\n\" \
    > par/parameters.txt"
rule build_target:
    input:
       "./my_target"
    output:
       touch("build_target.done")
    shell:
       "build_target -template format_nt -source {input} -target my_target"

如果我这样称呼：

snakemake -p -s snakefile

前三条规则正在执行，其他规则未执行。

我可以通过将其指定为参数来运行最后一条规则。

snakemake -p -s snakefile build_target

但我不知道我怎么能跑完所有。

非常感谢有关如何解决此问题的任何建议。

Answer 1

默认情况下，snakemake只执行snakefile的第一条规则。这是规则全部。为了生成规则所有的输入dynamic("task/{job}/taskOutput.tab")，它需要运行以下两个规则task和split_input，所以它确实如此。

如果你想要运行其他规则，你应该将它们的输出放在规则中，例如：

rule all:
    input: 
        dynamic("task/{job}/taskOutput.tab"),
        "par/parameters.txt",
        "build_target.done"

规则没有被snakemake拾取

1 个答案: