import sys
dfh=open(sys.argv[1],'r')
rfh=open(sys.argv[1]+'_no_intra_less.txt','w')
rfh2=open(sys.argv[1]+'_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[3]:
		if abs(int(line[1])-int(line[4]))<1000:
			intra_less+=1
			rfh2.write(line[0]+'\t'+line[1]+'\t'+line[3]+'\t'+line[4]+'\n')
		if abs(int(line[1])-int(line[4]))>=1000:
			intra_more+=1
			rfh.write(line[0]+'\t'+line[1]+'\t'+line[3]+'\t'+line[4]+'\n')
	if line[0]!=line[3]:
		inter+=1
		rfh.write(line[0]+'\t'+line[1]+'\t'+line[3]+'\t'+line[4]+'\n')

print "Total_mapped_reads intra_less, Intra_more, Inter: "+'\t'+str(Total)+'\t'+str(intra_less)+'\t'+str(intra_more)+'\t'+str(inter)