Question: bedtools and SNPs for pyseer-GWAS analysis
0
gravatar for 3335098459
4 weeks ago by
333509845920
333509845920 wrote:

Hi

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

Awan

snp gene genome • 128 views
ADD COMMENTlink modified 4 weeks ago by ATpoint26k • written 4 weeks ago by 333509845920
2
gravatar for ATpoint
4 weeks ago by
ATpoint26k
Germany
ATpoint26k wrote:

Failed to open 100010367.vcf.gz: could not load index

Is the file indexed? Did you try to re-index?

ADD COMMENTlink written 4 weeks ago by ATpoint26k

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.

ADD REPLYlink modified 29 days ago • written 29 days ago by 333509845920
1

Glad it worked. Our of interest, what did you use to index before? The reason you need bgzip over standard gzip is that bgzip does blockwise compression so tools can random-access the file by only decompressing the relevant part of that file.

ADD REPLYlink modified 29 days ago • written 29 days ago by ATpoint26k

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.

ADD REPLYlink written 29 days ago by 333509845920
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1361 users visited in the last hour