##################################
#                                #
# Last modified 2018/01/19       #
#                                #
# Georgi Marinov                 #
#                                # 
##################################

import sys
import string
import math
import random
from sets import Set
import time

def run():

    if len(sys.argv) < 4:
        print 'usage: python %s input fieldID seq5 seq3' % sys.argv[0]
        sys.exit(1)

    input = sys.argv[1]
    fieldID = int(sys.argv[2])
    seq1 = sys.argv[3]
    seq2 = sys.argv[4]

    linelist = open(input)
    for line in linelist:
        if line.startswith('#'):
            continue
        fields = line.strip().split('\t')
        seq = fields[fieldID]
        newseq = seq1 + seq + seq2
        outline = ''
        for i in range(fieldID):
            outline = outline + fields[i] + '\t'
        outline = outline + newseq + '\t'
        for i in range(fieldID+1,len(fields)):
            outline = outline + fields[i] + '\t'
        outline = outline.strip()
        print outline

run()
