|
0
|
1 png("/sonas-hs/bsr/hpc/data/yjin/test_BAMqc/exp/test1/figs/smp2.clipping_profile.png",width=500,height=500,units="px")
|
|
|
2 read_pos=c(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99)
|
|
|
3 count=c(9897,8881,8055,7532,7071,6644,6233,5790,5380,4983,4673,4419,4164,3892,3642,3403,3180,2938,2701,2456,2234,2000,1788,1577,1358,1179,998,799,620,447,306,184,92,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,77,152,275,404,546,700,889,1068,1241,1434,1632,1854,2085,2269,2480,2702,2931,3160,3397,3636,3894,4147,4415,4708,5084,5476,5908,6362,6820,7294,7783,8595,9576)
|
|
|
4 plot(read_pos,1-(count/80603),pch=20,xlab="Position of reads",ylab="Mappability",col="blue")
|
|
|
5 dev.state=dev.off()
|
