bedtools bamtofastq: Segmentation fault (core dumped)
1
0
Entering edit mode
2.1 years ago

I want to get all the reads from a bam file in a single fastq file. I am getting the message "Segmentation fault (core dumped)" while running

bedtools bamtofastq -i in.name.sorted.bam -fq out.fq

However, it works for the PE mode

bedtools bamtofastq -i in.name.sorted.bam -fq out.R1.fq -fq2 out.R2.fq

Bedtools version: v2.29.0

Any suggestion would be very helpful.

bedtools bamtofastq bam RNA-seq • 1.7k views
ADD COMMENT
0
Entering edit mode

It might be because your directory doesn't have enough space/memory or your program is trying to write in a file that is read only.

ADD REPLY
4
Entering edit mode
2.1 years ago

an alternative would be to try

samtools fastq

the bedtools bamtofastq also recommends that:

If you want to create a single, interleaved FASTQ file for paired-end data, you can just write both to /dev/stdout:

   bedtools bamtofastq -i x.bam -fq /dev/stdout -fq2 /dev/stdout > x.ilv.fq

Also, the samtools fastq command has more fucntionality and is a useful alternative.
ADD COMMENT
1
Entering edit mode

samtools fastq

worked perfectly.

ADD REPLY

Login before adding your answer.

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