starskyzheng / panpop

Application of pan-genome for population
MIT License
87 stars 8 forks source link

cal_aug_dp未输出{sample}.depth.txt.gz #57

Closed quanc1989 closed 3 months ago

quanc1989 commented 3 months ago

您好! 想请问一下在mergeSV3_pop.py的流程中,rule cal_aug_dp这一部分,在output中输出了如下的文件:

dpfile = '05_merge_samples/04.dps/{sample}.depth.txt.gz'

但是在实际运行的shell中,并未输出该文件:

echo '{wildcards.sample} {input.bam}' > {output.list_packpg} && \
perl {workflow.basedir}/scripts/cal_range_depth_aug.pl --vcfposs {input.poss} --list_packpg {output.list_packpg} --outdir {params.outdir} --threads {params.realthreads} --input_format bam --output_avgdpinfo {output.avgdp} --ref_only >> {log} 2>&1

由于该文件会在下一步rule fill_aug_dp中被作为输入使用,所以想请教一下该文件的内容及格式?

非常感谢!

quanc1989 commented 3 months ago

抱歉,我发现该文件会在04.dps文件夹下自动生成,并不需要指定文件名。