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 'fuck' not in line[6:20] and line[6:20].count('NA')<4:
		data=[]
		for a in line[6:20]:
			if a!='NA':
				data.append(float(a))
		if min(data)>=0.01:
			#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]!='fuck':
				if line[20]=='.':
					met_nonCpGi_unbound+=1
				if line[20]!='.':
					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)<=0.005:
			#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]!='fuck':
				if line[20]=='.':
					unmet_nonCpGi_unbound+=1
				if line[20]!='.':
					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)