##################################
#                                #
# Last modified 7/1/2009         # 
#                                #
# Georgi Marinov                 #
#                                # 
##################################

import sys
import string

from commoncode import *

def run():

    if len(sys.argv) < 3:
        print 'usage: python %s label SICERfilename bedfilename [-color r,g,b]' % sys.argv[0]
        sys.exit(1)

    factorlabel = sys.argv[1]
    SICERfilename = sys.argv[2]

    color = '0,0,0'
    if '-color' in sys.argv:
        color = sys.argv[sys.argv.index('-color') + 1]

    outfile = open(sys.argv[3], 'w')
    outfile.write('track name=%s visibility=4 color=%s\n' % (factorlabel, color))

    listoflines = open(SICERfilename)
    lineslist = listoflines.readlines()

    for line in lineslist:
        if line[0] == '#':
            continue
        fields=line.split('\n')[0].split('\t')
        outfile.write('%s\t%s\t%s\t%s\n' % (fields[0], fields[1], fields[2], fields[3]))
#        outfile.write(line)
    outfile.close()

run()