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

import sys
import string
import math
from sets import Set

def run():

    if len(sys.argv) < 3:
        print 'usage: python %s input fieldID outputfilename' % sys.argv[0]
        sys.exit(1)
    
    input = sys.argv[1]
    chrFieldID = int(sys.argv[2])
    outfilename = sys.argv[3]

    outfile=open(outfilename,'w')

    listoflines = open(input)
    for line in listoflines:
        fields = line.strip().split('\t')
        if line.startswith('#'):
            continue
        chr=fields[chrFieldID].split(':')[0]
        left=fields[chrFieldID].split(':')[1].split('-')[0]
        right=fields[chrFieldID].split(':')[1].split('-')[1]
        outline = ''
        for i in range(chrFieldID):
            outline = outline + 'fields[i]' + '\t'
        outline = outline + chr + '\t' + left + '\t' + right + '\t'
        for i in range(chrFieldID+1,len(fields)):
            outline = outline + fields[i] + '\t'
        outfile.write(outline.strip() + '\n')

    outfile.close()
   
run()