Plotting large file size of bedtools genomecov output
1
0
Entering edit mode
3.8 years ago
nitinra ▴ 50

Hi everyone,

I generated the output for coverage using bedrolls genomecov using this command:

bedtools genomecov -ibam input.sorted.bam > output.coverage -bga

The resulting files are huge (range of 5-6gb/file). I do not think it is possible for me to use R to plot the histogram since it is so huge. Is there a way to plot it using other tools on command line? Thank you for your help.

Here is my output:

Chromosome_1    0   24  0
Chromosome_1    24  52  2
Chromosome_1    52  82  3
Chromosome_1    118 202 2
Chromosome_1    202 268 1

Thank you!

bedtools genomecov genome coverage depth • 1.5k views
ADD COMMENT
3
Entering edit mode
3.8 years ago

use bedtools mergewith option -o and -c to aggregate/reduce the size of the bed file.

https://bedtools.readthedocs.io/en/latest/content/tools/merge.html

ADD COMMENT

Login before adding your answer.

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