png('/sonas-hs/bsr/hpc/data/yjin/test_BAMqc/exp/test1/figs/smp0.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(194,396,581,734,832,892,1002,1069,1105,1181,1242,1183,1154,1288,1353,1327,1235,1250,1297,1373,1405,1401,1412,1430,1406,1438,1485,1524,1512,1499,1528,1496,1466,1408,1392,1430,1394,1369,1373,1408,1487,1466,1455,1493,1461,1429,1461,1496,1449,1476,1466,1447,1440,1473,1503,1581,1606,1587,1605,1584,1609,1688,1685,1694,1665,1732,1750,1802,1891,1926,1939,1951,1927,1902,1933,1938,1990,1996,1988,2008,2022,2049,2036,2017,2016,2032,2004,2015,2020,2070,2118,2144,2161,2117,2084,2056,1916,1797,1488,963,285)
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()