GATK做vcf合并的模块应该是CombineVariants吧,具体你可以看一下这篇文章:https://www.omicsclass.com/article/1732,里面有写如何做多个vcf的合并
样本多建议每个样本先生成GVCF,再合并。
你贴的haplotypecaller命令生成的是VCF不是GVCF,下面的代码合并会会报错:
大样本 call SNP 生成GVCF这部分内容重测序课程里面有,你看看的:点击学习:https://bdtcd.xetlk.com/s/1VQOjQ
如果觉得我的回答对您有用,请随意打赏。你的支持将鼓励我继续创作!