CombineGVCFs报错:htsjdk.samtools.SAMFormatException: Did not inflate expected amount

老师们好,我最近在合并300份样品的gVCF,但是遇到了如题目所示的问题。

我的重测序物种染色体较长,我将每条染色体按350000000bp+剩余长度进行了分段,并分染色体进行gVCF合并。目前,只有一条染色体出现如题所示的问题(长度为350000000bp),其余染色体均正常运行(包括几条350000000bp长度的染色体)。

我的GATK版本是4.3.0.0,运行脚本如下 gatk CombineGVCFs --java-options "-Xms100g -Xmx100g -XX:ParallelGCThreads=4" --use-jdk-inflater true --use-jdk-deflater true --reference /path/to/part.fasta $(for i in {001..300}; do echo "--variant /path/to/part_gvcf/${i}.g.vcf.gz " ;done) --intervals chr2_part1 --output /path/to/part_gvcf/chr2_part1.g.vcf.gz --tmp-dir /path/to/temp,其中“--use-jdk-inflater true --use-jdk-deflater true”是我查看了GATK官网后加上的,但是依然报错。

报错信息如下:

attachments-2024-01-m4UQgjkJ6594023b19e50.png麻烦各位老师帮我解决一下,非常感谢!

请先 登录 后评论

2 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

这个报错没遇到过,新版的GATK不推荐用这个方法了,用:GenomicsDBImport 导入数据库,再 GenotypeGVCFs

具体可以看课程:https://bdtcd.xetslk.com/s/1VQOjQ

请先 登录 后评论
大耿老师

破案了,是有个样品的gvcf文件和tbi索引在这条染色体上有问题,重新获取了之后运行正常了,谢谢老师

请先 登录 后评论
  • 2 关注
  • 0 收藏,1607 浏览
  • 大耿老师 提出于 2024-01-02 20:34

相似问题