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 很重要。
本文来自博客园,作者:BioinformaticsMaster,转载请注明原文链接:https://www.cnblogs.com/koujiaodahan/p/16014731.html
posted on 2022-03-16 21:58 BioinformaticsMaster 阅读(353) 评论(0) 编辑 收藏 举报
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 震惊!C++程序真的从main开始吗?99%的程序员都答错了
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 单元测试从入门到精通
· 上周热点回顾(3.3-3.9)
· winform 绘制太阳,地球,月球 运作规律