Snakemake运行docker container命令

时间:2018-05-25 10:07:37

标签: docker snakemake

我有这个问题我有一个程序的docker镜像,我想使用该代码运行一些命令。 我有从输入到命令脚本的规则

rule create_maf:
    input:
        read="/data3/CALL/{sample}.muetect2.vt.clean.vcf.gz"
    output:
        out="/data3/MAF/{sample}.muetect2.vt.vep.maf"
   params:
         db_ens= "/mnt/mpwork/vep/homo_sapiens:/mnt/homo_sapiens",
         fst= "/mnt/homo_sapiens/87_GRCh37/Homo_sapiens.GRCh37.75.dna.primary_assembly.fa",
         tumor_id="{sample}",
         data="/home/data/TEST/:/data3"
    script:
        "scripts/report.py"

这是report.py

from subprocess import run

if len(snakemake.input) > 1:
    cmd = ['docker','container','run', '-it','-v', snakemake.params['db_ens'], '-v', snakemake.params['data'],'vcf2maf_87','perl','vcf2maf.pl','--input-vcf',snakemake.input['read'],snakemake.output['out'],'--vep-data', '/mnt/','--ref-fasta',snakemake.params['fst'],'--tumor-id',snakemake.params['tumor_id'],'--normal-id','NORMAL']
    run(cmd)

内部/ home / data / TEST /我有所有样本的CALL ......他们没有找到问题....

0 个答案:

没有答案