##################################
#                                #
# Last modified 10/15/2010        # 
#                                #
# Georgi Marinov                 #
#                                # 
##################################

import sys

try:
	import psyco
	psyco.full()
except:
	pass

def run():

    if len(sys.argv) < 2:
        print 'usage: python %s inputfilename outputfilename' % sys.argv[0]
        sys.exit(1)

    inputfilename = sys.argv[1]
    outputfilename = sys.argv[2]

    outfile = open(outputfilename, 'w')

    lineslist = open(inputfilename)
    i=0
    for line in lineslist:
        i+=1
        if i % 1000000 == 0:
            print i, 'lines processed'
        fields=line.strip().split('\t')
        try:
            chr=fields[0]
            start=fields[1]
            stop=fields[2]
            score=fields[3]
        except:
            print 'skipping line', line.strip()
            continue
        if '-' in fields[3]:
            outline=chr+'\t'+start+'\t'+stop+'\t'+score.split('-')[1]+'\n'
            outfile.write(outline)
        else:
            outline=chr+'\t'+start+'\t'+stop+'\t'+'-'+score+'\n'
            outfile.write(outline)

    outfile.close()

run()

