Generate Coverage Plot in R from Depth Data

Here is a simple way to take coverage data (coverage.depth) and make a plot in R to visualize it.  All you need is a column with base pair coordinates (V1) and a column with respective depth (V2) and you are all set to plug it in the below R code.

	# Read Data
	data=read.table("coverage.depth", as.is=T, header=F)
	# Chromosomal Plot
	plot(data$V1/1000000,data$V2, las=1, pch=20, col="white", xlab="Chromosomal Position (Mb)", ylab="Depth")
	arrows(data$V1/1000000,0,data$V1/1000000,data$V2, length=0, col="gray")
	abline(h=0)

view raw plot_depth.R hosted with ❤ by GitHub