import sys
import numpy
bin_size=sys.argv[2]
percent=float(sys.argv[3])
dfh=open(sys.argv[1],'r')
rfh=open(sys.argv[1]+'_'+bin_size+'_bp_bin_'+sys.argv[3]+'percentile.txt','w')
rfh.write('000000000000_Bin\t'+dfh.readline())
pre=''
temp_datas=[]
for i in range(0,14):
	temp_datas.append([])
bin_size=int(bin_size)
for i in dfh:
	line=i.split()
	index=int(line[0].split('_')[0])
	if 'NA' not in i:
		if index/bin_size!=pre:
			if pre!='':
				ID=str(pre*bin_size)
				for a in range(0,12-len(str(pre*bin_size))):
					ID='0'+ID
				rfh.write(ID+'\t'+'\t'.join(map(lambda x:str(numpy.percentile(x,percent)),temp_datas))+'\n')
			pre=index/bin_size
			temp_datas=[]
			for a in range(0,14):
				temp_datas.append([])
		data=map(lambda x:float(x.split('(')[0]),line[1:])
		map(lambda x,y:x.append(y),temp_datas,data)