This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How Can I Plot A Region'S Depth As Follow? While Tools Should I Use?

X is the position in genome, and Y is the read depth . while tools can do this work ?

depth-of-coverage plot

2 answers

It really depends on how comfortable you are in using plotting tools. You can use R or gnuplot and can make pretty cool figures. You can bin the genomic locations so that you can cover the whole chromosome in one plot. I generally use circos (http://circos.ca/tutorials/course/) and plotting tools from R. If you want to do it in an easy way you can use:

  1. IGV (http://www.broadinstitute.org/igv/).
  2. Seqmonk (http://www.bioinformatics.babraham.ac.uk/projects/seqmonk/)
  3. Grid Package in R (http://biobeat.wordpress.com/2012/12/18/plot-genome-coverage-graphs-using-grid-package-in-r/)

PS: You should wait for other people's reply too. They may suggest something really cool. I am not into plotting things.

You can try Gviz R/bioconductor https://www.bioconductor.org/packages/devel/bioc/vignettes/Gviz/inst/doc/Gviz.html

It cab be used to to plot Sequencing Depth directly from bam file.

enter image description here

Best wishes,

-Ahmed Alhendi

Log in to answer this question.