import sys
RS=sys.argv[1]

ref=open('/pbld/netapp/database/human/hg19/hg19.fa','r')
line=ref.readline().rstrip()
while line:
	data=[]
	if '>' in line:
		rfh=open(line[1:]+'_'+RS+'.txt','w')
		chro=line[1:]
		count=0
		line=ref.readline().rstrip()
		continue
	elif '>' not in line[0]:
		line=line.upper()
		line2=ref.readline().rstrip()
		if '>' not in line2:
			line=(line+line2).upper()
		if len(line)==100:
			for i in range(0,50):
				if line[i:i+len(RS)]==RS:
					rfh.write(str(count+i+1)+'\n')
		elif len(line)==50:
			for i in range(0,50-RS+1):
			        if line[i:i+len(RS)]==RS:
			                rfh.write(str(count+i+1)+'\n')
		line=line2
		count+=50