Hi,
I would like to plot the coverage distribution which gives the fraction of genome /target regions covered at different depths say 1,5,20,30,50....
An example of the plot is shown in the attachment which is generated from R package that consumes huge amount of time.
I would like to accomplish this in short time, and so i have used genomeCoverageBed and coverageBed from Bedtools and have the output coverage values from those two functions.
Could someone provide a script in R or any other language to generate a similar kind of plot show in the attachment?
I would like to plot the coverage distribution which gives the fraction of genome /target regions covered at different depths say 1,5,20,30,50....
An example of the plot is shown in the attachment which is generated from R package that consumes huge amount of time.
I would like to accomplish this in short time, and so i have used genomeCoverageBed and coverageBed from Bedtools and have the output coverage values from those two functions.
Could someone provide a script in R or any other language to generate a similar kind of plot show in the attachment?