vcf2maf流程中断后要删除vep.vcf文件

vcf2maf批量把vcf转化成maf文件,使用的是本人写的pipeline,会自动记录哪些步骤已run,中断后重新开始不需要再run
通过多进程方法,能同时运行多个pipeline,还能提示一个步骤运行了多少时间。

因为挂新硬盘的原因,关机挂载硬盘,转移数据,重新运行脚本,按理说在关机前被中断的5个vcf2maf过程会重新开始,然而发现这5个vcf2maf的时间显著低于前面样本同步骤,且生成的maf大小明显偏低。

关机,挂硬盘折腾了一番,仍然如此,头大。

后来想着干脆全部样本重新跑试试 ,作备份时,把待转的vcf文件目录下,在vcf2maf过程中由vep生成的vep.vcf文件全部转移掉,再怀着试一试的想法重新运行vcf2maf,发现时间应该正常了!!

猜想:可能是vcf2maf调用 vep前,如果检测到同目录下有vcf文件对应的vep.vcf文件,就认为已经完成vep注释过程,vcf2maf会直接从vep.vcf文件转成maf文件。起初被中断的5个样本,并没有完成相应的vep注释就直接进入生成maf过程,从而出现一开始描述的问题。

解决方法:及时把生成的vep.vcf文件转移或删除

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

友情链接更多精彩内容