我有这个问题我有一个程序的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"
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 ......他们没有找到问题....