Question: Plotting large file size of bedtools genomecov output
0
gravatar for nitinra
4 months ago by
nitinra10
nitinra10 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 4 months ago by Pierre Lindenbaum133k • written 4 months ago by nitinra10
3
gravatar for Pierre Lindenbaum
4 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum133k 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 4 months ago by Pierre Lindenbaum133k
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: 1178 users visited in the last hour