Question: Plotting large file size of bedtools genomecov output
0
gravatar for nitinra
5 weeks ago by
nitinra0
nitinra0 wrote:

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!

ADD COMMENTlink modified 5 weeks ago by Pierre Lindenbaum131k • written 5 weeks ago by nitinra0
3
gravatar for Pierre Lindenbaum
5 weeks ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum131k wrote:

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 COMMENTlink written 5 weeks ago by Pierre Lindenbaum131k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2189 users visited in the last hour