bedtools segmentation fault
Hi all,
I have a multi-sample VCF file and the Kaviar VCF file. I'm doing:
bedtools intersect -sorted -a Family5_sorted.vcf.gz -b Kaviar-150812-ISB-hg19.vcf.gz -wb
It prints 9795 results for chromosome 1 and then dies with a Segementation fault (core dumped). The line that it dies on doesn't seem to be any different than previous lines. I've also tried reversing the a and b file and it dies on a different line (but close to where it dies the other way avb is 627982 and bva is 625369). How do I troubleshoot what's going on?
Thanks,
Denise
• 4,790 views
•
link
updated
by
Ram
•
written
by
denise.mauldin •
•
written
by
denise.mauldin •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Visualizing a mapped SNPs/INDELs on chromosomes of chicken
written by mohsamir2016 •
-
Opening a VCF.bgz file in python
written by Eliza •
-
Biopython byte positions are not compatible with bgzip
written by b10hazard •
-
Segmentation fault (core dumped) during vg surject
written by LilithElina •
-
Samtools segmentation fault when extracting exons from bam file.
written by jstevenson2256 •
-
Problem with bedtools subtract
written by zion22
-
pysam segementation fault with multiprocessing
written by wolffj •
-
bcftools call segmentation fault
written by AP
-
Please help me with samtools
written by seta
-
BedTools IntersectBed resulting in Segmentation Fault with BedFile
written by ali •
Normally this error is associated with low memory (RAM). Try some of the suggested solutions in this post (https://www.biostars.org/p/124620/).
Had same problem (as Ashutosh mentioned it was caused by low memory). Either use more powerful machine or try looping through chromosomes.