png('/sonas-hs/bsr/hpc/data/yjin/test_BAMqc/exp/test1/figs/smp1.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(200,400,561,726,832,920,1032,1085,1102,1152,1242,1245,1193,1258,1260,1306,1238,1246,1259,1290,1318,1304,1374,1410,1391,1477,1438,1445,1468,1418,1457,1443,1382,1404,1422,1489,1478,1445,1467,1484,1486,1446,1485,1474,1468,1433,1448,1466,1466,1472,1493,1446,1465,1507,1489,1536,1628,1639,1606,1638,1612,1612,1593,1631,1672,1715,1767,1782,1832,1881,1866,1877,1938,1992,1969,1940,1936,1936,1960,1988,1949,1964,1955,1988,1984,2020,2055,2110,2140,2100,2106,2128,2135,2150,2146,2062,1995,1875,1549,982,284)
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()