import sys
fn=sys.argv[1]
dfh=open(fn,'r')
dfh.readline()
rfh_met=open('methylated_selex.txt','w')
rfh_unmet=open('unmethylated_selex.txt','w')
met_CpGi_bound=0
met_CpGi_unbound=0
met_nonCpGi_bound=0
met_nonCpGi_unbound=0
unmet_nonCpGi_bound=0
unmet_nonCpGi_unbound=0
unmet_CpGi_bound=0
unmet_CpGi_unbound=0
for i in dfh:
	line=i.split()
	if 'NA' not in line[1:15]:
		data=map(lambda x:float(x),line[1:15])
		if min(data)>=80:
			rfh_met.write(line[0]+'\t'+line[0]+'\t'+'\t'.join(map(lambda x:str(x/100),data))+'\t'+str(100*int(line[-1]))+'\n')
			if line[-1]=='0':
				if line[-3]=='0':
					met_nonCpGi_unbound+=1
				if line[-3]=='1':
					met_nonCpGi_bound+=1
			if line[-1]=='1':
				if line[-3]=='0':
					met_CpGi_unbound+=1
				if line[-3]=='1':
					met_CpGi_bound+=1
		if max(data)<=20:
			rfh_unmet.write(line[0]+'\t'+line[0]+'\t'+'\t'.join(map(lambda x:str(x/100),data))+'\t'+str(100*int(line[-1]))+'\n')
			if line[-1]=='0':
				if line[-3]=='0':
					unmet_nonCpGi_unbound+=1
				if line[-3]=='1':
					unmet_nonCpGi_bound+=1
			if line[-1]=='1':
				if line[-3]=='0':
					unmet_CpGi_unbound+=1
				if line[-3]=='1':
					unmet_CpGi_bound+=1
print 'met_CpGi_bound: '+str(met_CpGi_bound)
print 'met_CpGi_unbound: '+str(met_CpGi_unbound)
print 'unmet_CpGi_bound: '+str(unmet_CpGi_bound)
print 'unmet_CpGi_unbound: '+str(unmet_CpGi_unbound)
print 'met_nonCpGi_bound: '+str(met_nonCpGi_bound)
print 'met_nonCpGi_unbound: '+str(met_nonCpGi_unbound)
print 'unmet_nonCpGi_bound: '+str(unmet_nonCpGi_bound)
print 'unmet_nonCpGi_unbound: '+str(unmet_nonCpGi_unbound)