基于图形泛基因组文件call取sv的时候,运行命令nohup vg call new_index.giraffe.gbz -r new_index.snarls -k D23.pack -a -A -t 30 -s D23 > D23.vcf 2>call.o &,显示done,但生成的vcf文件里面没有变异位点,只有个头,如截图中所示。这种问题怎么解决呢?

attachments-2025-12-apOQNOCb695287a6b7ece.jpg

请先 登录 后评论

1 个回答

Ti Amo

生成到 call.o里面去了。如果是这种 输出文件用>指定的,写成sh脚本去运行

请先 登录 后评论