##################################
#                                #
# Last modified 04/08/2015       # 
#                                #
# Georgi Marinov                 #
#                                # 
##################################

import sys

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

def run():

    if len(sys.argv) < 2:
        print 'usage: python %s config outfilename [-minSeqLen number]' % sys.argv[0]
        print '\tconfig format: label <tab> filename <tab> fieldID'
        sys.exit(1)

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

    outfile = open(outputfilename, 'w')

    lineslist  = open(config)
    i=0
    for line in lineslist:
        if line.startswith('#'):
            continue
        fields = line.split('\t')
        label = fields[0]
        file = fields[1]
        fieldID = int(fields[2])
        i = 0
        lines = open(file)
        for ll in lines:
            i+=1
            seq = ll.strip().split('\t')[fieldID]
            outline = '>' + label + '_' + str(i)
            outfile.write(outline + '\n')
            outfile.write(seq + '\n')

    outfile.close()

run()

