Issues with bcftools
0
0
Entering edit mode
8 months ago
hgrsce • 0

Hi, I can’t seem to solve this issue. I produced .gz versions of my .vcf files, and I indexed them, however bcftools merge does not seem to work any way I try. Example:

bcftools merge -l list.txt > merged.vcf

Produces:

Failed to open example.final_comp.vcf: unknown file type

  Thank you for your help.

Bcftools • 624 views
ADD COMMENT
0
Entering edit mode

Could you try by giving the absolute paths to VCFs in your list.txt?

ADD REPLY
0
Entering edit mode

I tried that, also. The issue persists.

ADD REPLY
0
Entering edit mode

check example.final_comp.vcf is a non-empty valid VCF.

furthermore, VCF files must be compressed and indexed

ADD REPLY
0
Entering edit mode

It is not empty, and it's valid. I also compressed and indexed the files. The same issue still occurs.

ADD REPLY
0
Entering edit mode

show me output of the following command:

cat list.txt | while read F; do file "${F}" && bcftools view "${F}" > /dev/null; done
ADD REPLY

Login before adding your answer.

Traffic: 1835 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6