当我对两个基因型文件位置取交集，并重新生成两个vcf：
$ bcftools view -R overlap.lst variant.filter.vcf.gz -Oz -o 300.vcf.gz
出现如下错误：
$ tabix 300.vcf.gz
[E::hts_idx_push] Unsorted positions on sequence #4: 29013869 followed by 29013853
tbx_index_build failed: 300.vcf.gz
原因是位点未排序。因此可知bcftools从vcf取位点子集时，并不会对位点进行排序。
解决如下：
zcat 300.vcf.gz |grep '^#'>header
zcat 300.vcf.gz |grep -v '^#' |sort -k1,1d -k2,2n > body
cat header body |bgzip -c >new.300.vcf.gz
tabix new.300.vcf.gz
hive 建表报错：ParseException - cannot recognize input near 'e...
使用tabix构建vcf的索引时，遇到失败的问题tbx_index_build failed: XXX.vcf....
mysql 1、多表连接2、子查询3、索引都有哪些，索引能解决什么问题 ，如何创建索引，建索引有什么原则，索引的底...
报错:TypeError: Cannot read property 'getNodeByValue' of un...
tabix 可以对NGS分析中常见格式的文件建立索引，从而加快访问速度，不仅支持VCF文件，还支持BED, GFF...
tabix 可以对NGS分析中常见格式的文件建立索引，从而加快访问速度，不仅支持VCF文件，还支持BED, GFF...
我们获得了bed文件之后，需要对它们建立索引，实现快速定位查找
mysql的建表一般要建索引，那索引与数据的结构是怎样的呢？这很重要，清楚它能解决一些场景问题。 数据库索引相当于...
欢迎关注"生信修炼手册"！ tabix 可以对NGS分析中常见格式的文件建立索引，从而加快访问速度，不仅支持VCF...
es版本6.8.4解决方式：对字段索引两次，一次索引分词（用于搜索）一次索引不分词（用于排序） 字符串类型排序报错...
本文标题：解决tabix建索引报错[E::hts_idx_push] Un
本文链接：https://www.haomeiwen.com/subject/szoordtx.html
网友评论