在snakemake通配符中使用正则表达式

时间:2017-01-04 14:07:23

标签: snakemake

我在蛇形通配符中使用正则表达式,但我遇到了一个我不明白的错误。

在这个缩短的例子中,它起作用:

rule graphviz:
    input: "{graph}.dot"
    output: "{graph}.{ext,(pdf|png)}"
    shell: "dot -T{wildcards.ext} -o {output} {input}"

在这个例子中,它没有:

## This is working
rule fastqc:
    input: "{reads}.fastq"
    output: "{reads}_fastqc/{sample}_fastqc.html"
    shell:"fastqc --format fastq {input}"

## This is not working
rule fastqc:
    input: "{reads}.{ext,(bam|fastq)}"
    output: "{reads}_fastqc/{sample}_fastqc.html"
    shell:"fastqc --format {wildcards.ext} {input}"

我附上了我收到的错误消息的屏幕截图。谢谢你的帮助。

Here is a screencap of the error message.

0 个答案:

没有答案