This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Plot The Coverage Of A Region In The Genome

HI,

I would like to plot the coverage distribution of the targeted region and the genome. For this i have used genomeCoverageBed and coverageBed functions from Bedtools and got the output from the two functions.

Could some one suggest a way or Rscript to plot 1. coverage distribution for target region (X-axis:depths=1,5,10... and Yaxis:fraction of target bases covered)

  1. coverage for genome (X-axis:depths=1,5,10... and Yaxis:fraction of target bases covered)

from the Bedtools outptut?

coverage next-gen
• 3,980 views • link
updated by Sean Davis Sean Davis • written by hellbio hellbio

how exactly your data (to be R input) looks like?

• 0 views • link
written by Pavel Senin Pavel Senin

1 answer

See the GenomicRanges package for dealing with aligned reads and regions. In particular, you may find the coverage() function useful.

• 0 views • link
written by Sean Davis Sean Davis

But, working in R with huge datasets needs huge memory and time. So is there any alternate way?

• 0 views • link
written by hellbio hellbio

Log in to answer this question.

More posts like this