BAM to Bedgraph creating an empty file
0
0
Entering edit mode
19 days ago
irebekah.c • 0

Hi,

I downloaded a fastq file from a paper and was trying to create a Bedgraph file.

After mapping using bismark, I got a 12Gb bam file. And after deduplication, I used bismark_methylation_extractor to create the bedgraph file. There was no error but the Bedgraph file was only 40b and showed nothing. Also, thousands of temp files were created but not deleted.

Here is my code:

bismark --non_directional --genome ref/ trim/filename_trimmed.fq --multicore 2 -o map/
deduplicate_bismark --bam map/filename_trimmed_bismark_bt2.bam 
bismark_methylation_extractor --gzip --bedGraph filename_trimmed_bismark_bt2.deduplicated.bam

Could someone help me please?

Bisulfite-seq bismark • 299 views
ADD COMMENT
0
Entering edit mode

You have a space in your reference path. Not sure it's intended?

ADD REPLY
0
Entering edit mode

Hi, Thank you for replying. I have about 2.5TB free space there and I think that is enough? And I also tried with mapped-only reads or limited multicore but it doesn’t work anyway.

This is how it reported:

Now writing methylation information for file >>filename.deduplicated.mapped.txt.gz<< to individual files for each chromosome

Failed to open filehandle: Too many open files at /usr/local/bin/bismark2bedGraph line 278, <IN> line 44246607. Finished BedGraph conversion ...

ADD REPLY

Login before adding your answer.

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