Reference allele mismatch at chr1:1630781 .. REF_SEQ:'C' vs VCF:'T' Ref alleles mismatch between: [VC Unknown @ chr1:1630781 Q. of type=SNP alleles=[T*, C] attr={} GT=[[HG00472 T*/C]] ..

原因是待处理的vcf,他们的ref allele不一致,所以

解决方法是 bcftools norm  先统一 ref allele

如bcftools norm -c s -f /ldfssz1/ST_BIGDATA/USER/st_bigdata/Sentieon/reference_bigdatacompute/hg38_noalt_withrandom/hg38.fa -o /jdfssz1/ST_HEALTH/P21Z10200N0047/lizhichao/zbolt_test/data_input/1Kgenome/microarray_vcf/ALL.chip.omni_broad_sanger_combined.20140818.snps.genotypes.hg38.1_22.sort.vcf.reheaderzbolt.reref.gz -O z /jdfssz1/ST_HEALTH/P21Z10200N0047/lizhichao/zbolt_test/data_input/1Kgenome/microarray_vcf/ALL.chip.omni_broad_sanger_combined.20140818.snps.genotypes.hg38.1_22.sort.vcf.reheaderzbolt.gz && echo "re_ref done"

 

参数:-c s 很重要。

 

posted on   BioinformaticsMaster  阅读(353)  评论(0编辑  收藏  举报

相关博文:
阅读排行:
· 震惊!C++程序真的从main开始吗?99%的程序员都答错了
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 单元测试从入门到精通
· 上周热点回顾(3.3-3.9)
· winform 绘制太阳,地球,月球 运作规律

导航

< 2025年3月 >
23 24 25 26 27 28 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30 31 1 2 3 4 5
点击右上角即可分享
微信分享提示