Mercurial > repos > nilesh > rseqc

view test-data/output.clipping_profile.r @ 50:f242ee103277 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/lparsons/galaxy_tools/tree/master/tools/rseqc commit 91ad241aa3f34b70649d13a5f18611da7577a5ee
author lparsons
date Tue, 03 May 2016 16:36:57 -0400
parents 6b33e31bda10
children 09846d5169fa
line wrap: on
line source

pdf("output.clipping_profile.pdf")
read_pos=c(0,1,2,3,4,5,6,7,8,9,44,45,46,47,48,49,50)
count=c(16,12,11,8,6,5,1,1,1,1,1,2,2,2,3,4,4)
plot(read_pos,1-(count/40),col="blue",main="clipping profile",xlab="Position of reads",ylab="Mappability",type="b")
dev.off()