Question: bcftools call segmentation fault
1
gravatar for AP
2.9 years ago by
AP90
AP90 wrote:

Hi,

I am currently trying to run bcftools call on a large bcf file and I keep receiving the "Segmentation fault" message. I use the latest version of bcftools (1.3.1). I am using: bcftools call -v -c -f gq in.bcf > out.bcf

I successfully used bcftools call in the past on various files and never had this problem before. The only difference here is that I slightly changed my workflow and now filter for depth directly on my bcf file, not on my vcf. Here is my workflow

  1. samtools mpileup on a set of bam files
  2. vcftools --bcf file.bcf --min-meanDP 5 --max-meanDP 25 --recode-bcf --out file_depth_filtered
  3. bcftools call -v -c -f gq file_depth_filtered.bcf > out.bcf

In the past, I use to make a vcf file from my bcf right after step 1 and filtered directly on my vcf file. This worked just fine. Also, using bcftools call on my none-filtered bcf file works well. So basically, there is something wrong in my filtered file.

The error I receive from my HPC is "line 11: 3165 Segmentation fault (core dumped) bcftools call -v -c -f gq file_depth_filtered,bcf > out.bcf

What do you think could be causing this?

Thanks!

error call bcftools segfault • 1.4k views
ADD COMMENTlink written 2.9 years ago by AP90

The error I receive from my HPC is "line 11: 3165 Segmentation fault (core dumped) bcftools call -v -c -f gq file_depth_filtered.bcf > out.bcf"

No comma before bcf in "file_depth_filtered.bcf"

ADD REPLYlink written 2.9 years ago by AP90
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: 925 users visited in the last hour