png('/sonas-hs/bsr/hpc/data/yjin/test_BAMqc/exp/test1/figs/smp2.geneBodyCoverage.png',width=500,height=500,units='px')
x=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,100)
y=c(172,322,489,620,693,723,766,838,861,939,1020,1028,988,1055,1122,1087,1031,1082,1111,1119,1184,1144,1167,1176,1145,1143,1129,1185,1192,1207,1225,1185,1213,1241,1268,1313,1240,1204,1190,1235,1269,1295,1328,1308,1229,1220,1252,1273,1272,1295,1283,1335,1335,1334,1377,1385,1434,1411,1407,1424,1420,1450,1458,1483,1511,1560,1583,1627,1669,1710,1777,1792,1812,1825,1802,1862,1857,1900,1880,1874,1867,1859,1869,1859,1851,1897,1883,1971,2030,1992,1971,2017,2030,1983,1995,1959,1896,1787,1461,980,288)
smoothsp = smooth.spline(x,y,spar=0.35)
plot(smoothsp,type="l",col="blue",xlab="Percentile of Gene Body (5'->3')",ylab="Number of read",xlim=c(0,100))
dev.state = dev.off()