line 164: 15342 Segmentation fault (core dumped) ${cmd} ${options} - 0<&0
1
0
Entering edit mode
9 months ago
Niharika • 0

command I ran:

vcf2bed-megarow < my_sample.vcf.gz --max-mem 127G --sort-tmpdir temp > my_sample.bed

I am using version: 2.4.41 of vcf2bed size of compressed vcf : 48.9 mb

vcf2bed megarow bedops • 785 views
ADD COMMENT
1
Entering edit mode

what is vcf2bed-megarow ?

ADD REPLY
0
Entering edit mode

It is part of BEDOPS.

ADD REPLY
2
Entering edit mode
9 months ago

You can't pass a gzipped file to vcf2bed (or vcf2bed-megarow, or any of the convert2bed format translators).

You should extract the gzip file and move the options to the command:

vcf2bed-megarow --max-mem 127G --sort-tmpdir temp < <(gunzip -c my_sample.vcf.gz) > my_sample.bed

Or the equivalent:

gunzip -c my_sample.vcf.gz | vcf2bed-megarow --max-mem 127G --sort-tmpdir temp - > my_sample.bed

You can almost certainly reduce this parameter: --max-mem 127G to something much much lower. It should be less than available system memory, in any case.

ADD COMMENT

Login before adding your answer.

Traffic: 1503 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