|
0
|
1 png("/sonas-hs/bsr/hpc/data/yjin/test_BAMqc/exp/test1/figs/smp0.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(10290,9294,8453,7955,7448,6951,6512,6050,5602,5219,4931,4682,4396,4114,3827,3540,3283,3032,2781,2555,2301,2073,1849,1638,1411,1212,1012,799,634,485,341,201,104,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,108,203,338,496,660,807,995,1203,1420,1604,1824,2020,2220,2460,2668,2901,3139,3379,3631,3895,4137,4426,4691,4976,5369,5810,6205,6686,7174,7678,8224,9087,10073)
|
|
|
4 plot(read_pos,1-(count/77686),pch=20,xlab="Position of reads",ylab="Mappability",col="blue")
|
|
|
5 dev.state=dev.off()
|
