Snakemake规则

Question

我想用snakemake制作一个生物信息学管道，我在谷歌上搜索了它，阅读了文档和其他东西，但我仍然不知道如何让它工作。

下面是我的一些原始数据文件。

原始数据/010__bua_1.fq.gz，原始数据/010__bua_2.fq.gz

原始数据/11_15_ap_1.fq.gz，原始数据/11_15_ap_2.fq.gz

...they都是成对的文件。)

这是我的align.snakemake

from os.path import join


STAR_INDEX = "/app/ref/ensembl/human/idx/"
SAMPLE_DIR = "Rawdata"
SAMPLES, = glob_wildcards(SAMPLE_DIR + "/{sample}_1.fq.gz")

R1 = '{sample}_1.fq.gz'
R2 = '{sample}_2.fq.gz'


rule alignment:
    input:
        r1 = join(SAMPLE_DIR, R1),
        r2 = join(SAMPLE_DIR, R2),
    params:
        STAR_INDEX = STAR_INDEX
    output:
        "Align/{sample}.bam"
    message:
        "--- mapping STAR---"
    shell:"""
        mkdir -p Align/{wildcards.sample}
        STAR --genomeDir {params.STAR_INDEX} --readFilesCommand zcat  --readFilesIn {input.r1} {input.r2}  --outFileNamePrefix Align/{wildcards.sample}/log
"""

就是这个。我通过"snakemake -np -s align.snakemake“运行这个文件，得到了这个错误。

WorkflowError:目标规则不能包含通配符。请指定具体文件或不带通配符的规则。

很抱歉我问了这个问题，虽然有很多人用得很好。任何帮助都是非常有价值的。对不起，我的英语不好。

附注:我阅读了官方文档和教程，但仍然一无所知。

Answer 1

哦，我说过了。这是我对我的问题的回答，有些人可能需要一些帮助。

  from os.path import join


STAR_INDEX = "/app/ref/ensembl/human/idx/"
SAMPLE_DIR = "Rawdata"
SAMPLES, = glob_wildcards(SAMPLE_DIR + "/{sample}_1.fq.gz")

R1 = '{sample}_1.fq.gz'
R2 = '{sample}_2.fq.gz'


rule all:
    input:
        expand("Align/{sample}/Aligned.toTranscriptome.out.bam", sample=SAMPLES)

rule alignment:
    input:
        r1 = join(SAMPLE_DIR, R1),
        r2 = join(SAMPLE_DIR, R2)
    params:
        STAR_INDEX = STAR_INDEX
    output:
        "Align/{sample}/Aligned.toTranscriptome.out.bam"
    threads:
        8
    message:
        "--- Mapping STAR---"
    shell:"""
        mkdir -p Align/{wildcards.sample}
        STAR --genomeDir {params.STAR_INDEX} --outSAMunmapped Within --outFilterType BySJout --outSAMattributes NH HI AS NM MD --outFilterMultimapNmax 20 --outFilterMismatchNmax 999  --outFilterMismatchNoverLmax 0.04  --alignIntronMin 20  --alignIntronMax 1000000  --alignMatesGapMax 1000000  --alignSJoverhangMin 8  --alignSJDBoverhangMin 1  --sjdbScore 1  --runThreadN {threads}  --genomeLoad NoSharedMemory  --outSAMtype BAM Unsorted  --quantMode TranscriptomeSAM  --outSAMheaderHD \@HD VN:1.4 SO:unsorted  --readFilesCommand zcat  --readFilesIn {input.r1} {input.r2}  --outFileNamePrefix Align/{wildcards.sample}/log
"""