import sys
dfh=open(sys.argv[1],'r')
rfh=open(sys.argv[1]+'no_intra_less.txt','w')
intra_less=0
intra_more=0
inter=0
Total=0
for i in dfh:
	Total+=1
	line=i.split()
	if line[0]==line[2]:
		if abs(int(line[1])-int(line[3]))<1000:
			intra_less+=1
		if abs(int(line[1])-int(line[3]))>=1000:
			intra_more+=1
			rfh.write(i)
	if line[0]!=line[2]:
		inter+=1
		rfh.write(i)
print "Total_mapped_reads intra_less, Intra_more, Inter: "+'\t'+str(Total)+'\t'+str(intra_less)+'\t'+str(intra_more)+'\t'+str(inter)