输入与输出：作业成功完成，但是缺少一些输出文件

问题描述

我认为我的管道失败了，因为预期的全部规则输出与实际最终输出之间存在冲突。我相信snakemake正在等待文件kma/{sample}没有扩展名出现，相反，它正在获取包含multiext("kma/{sample}",".res",".aln",".fsa",".gz")的目录，而我很难使它们正常播放。

configfile: "config.yaml"

rule all:
    input:
        expand("kma/{sample}",sample = config["samples"])
        #multiext("kma/{sample}",".gz",sample = config["samples"])

rule seqtk_qualtiy_filter:
    input:
        lambda wildcards: "S5_Raw/" + config["samples"][wildcards.sample]
    output:
        temp("qtrim/{sample}.qtrim.fq")
    shell:
        "seqtk trimfq -b 0.01 {input} > {output}"

rule seqtk_clip:
    input:
        "qtrim/{sample}.qtrim.fq"
    output:
        temp("clip/{sample}.clip.fq")
    shell:
        "seqtk trimfq -b20 -L 350 {input} > {output}"

rule bbnorm:
    input:
        "clip/{sample}.clip.fq"
    output:
        "S5_processed/{sample}.norm.fq"
    shell:
        "bbnorm.sh in={input} out={output} target=100"

rule kma_map:
    input:
        "S5_processed/{sample}.norm.fq"
    params:
        ref = "ref/consensus.fasta"
    output:
        directory("kma/{sample}")
        #multiext("kma/{sample}",".gz")
    shell:
        "kma -i {input} -t_db {params.ref} -o {output}"

如果您以kma希望看到处理完成的方式运行该错误，则该错误

Waiting at most 5 seconds for missing files.
MissingOutputException in line 33 of /home/sean/Desktop/reo/antisera project/ReovirusS1AmpliconS5.smk:
Job completed successfully,but some output files are missing. Missing files after 5 seconds:
kma/BA8359-19
This might be due to filesystem latency. If that is the case,consider to increase the wait time with --latency-wait.
  File "/home/sean/.local/lib/python3.6/site-packages/snakemake/executors/__init__.py",line 544,in handle_job_success
  File "/home/sean/.local/lib/python3.6/site-packages/snakemake/executors/__init__.py",line 231,in handle_job_success
Shutting down,this might take some time.

我也尝试过增加延迟时间，但是由于从未真正创建预期的文件，因此等待多长时间都没有关系。

使用multiext函数时收到的错误

Error in rule kma_map:
    jobid: 11
    output: kma/BA8359-19.res,kma/BA8359-19.aln,kma/BA8359-19.fsa,kma/BA8359-19.frag.gz
    shell:
        kma -i S5_processed/BA8359-19.norm.fq -t_db ref/consensus.fasta -o kma/BA8359-19.res kma/BA8359-19.aln kma/BA8359-19.fsa kma/BA8359-19.frag.gz
        (one of the commands exited with non-zero exit code; note that snakemake uses bash strict mode!)

Shutting down,this might take some time.
Exiting because a job execution failed. Look above for error message

我相信会发生此错误，因为snakemake在-o标志中列出了所有文件类型，而kma则拒绝。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）