添加链接
link管理
链接快照平台
  • 输入网页链接,自动生成快照
  • 标签化管理网页链接

I have used snippy package for SNP calling and later on used those for pyseer-GWAS analyses. But it seems that this command of bedtools cannot read the vcf files.

bcftools merge -m none -0 -O z *.vcf.gz > merged.vcf.gz
I am getting this error:
bcftools: Relink /root/miniconda3/bin/../lib/./libgfortran.so.4' with /lib/x86_64-linux-gnu/librt.so.1' for IFUNC symbol `clock_gettime'
Failed to open 100010367.vcf.gz: could not load index

Is there anyone having this kind of experience with pyseer before? Is there any problem with bedtools or vcf files?

Regards

After your reply, I have used other packages to generate indexed vcf file.

vcfcombine *.vcf > merged.vcf (I guess there was some issue with bedtools)

bgzip -c merged.vcf > snps.vcf.gz

tabix -p vcf snps.vcf.gz

Thanks for pointing out the direction.

Actually there was some issue with the bedtools/bcftools compiling code and installation. Previous wrapper tool also used bedtools for indexing and compression. While re-indexing I figured out this point and I used different tool to combine, compress and indexing. Also thanks for explaining the bgzip core concept.