snakemake wrapper Gatk/Mutect2

Viewed 25

In the Snakemake Wrapper documentation

rule mutect2_bam:
    input:
        fasta="genome/genome.fasta",
        map="mapped/{sample}.bam",
    output:
        vcf="variant_bam/{sample}.vcf",
        bam="variant_bam/{sample}.bam",
    message:
        "Testing Mutect2 with {wildcards.sample}"
    threads: 1
    resources:
        mem_mb=1024,
    log:
        "logs/mutect_{sample}.log",
    wrapper:
        "v1.10.0/bio/gatk/mutect"

Is just only one input bam file, but if I have to use two bam files. What should I do?

0 Answers
Related