python10

# -*- coding: utf-8 -*-
#!/usr/bin/env python

import os
import sys

fileSQ = sys.argv[1] 
fileSD = sys.argv[2] 
outmergy = sys.argv[3]
outcon = open(outmergy,"w")

dictSQ = {}
with open(fileSQ,"r") as filSQliness:
    filSQlines = filSQliness.readlines()
    for filSQline in filSQlines:
        filSQlin = filSQline.strip().split("\t")
        keySQ = str(filSQlin[0]).strip().split("|")[0]
        dictSQ[keySQ] = ""
        for filSQli in filSQlin[1:]:
            dictSQ[keySQ] += str(filSQli) + "\t"

print(str(dictSQ))

with open(fileSD,"r") as fileSDliness:
    fileSDlines = fileSDliness.readlines()
    headlist = fileSDlines[0].strip().split(" ")
    conhead = ""
    for headlisti in headlist:
        conhead +=str(headlisti)+"\t"
    outcon.write("genename\tfit_value\tPredicted\tpos\t"+str(conhead)) 
    for fileSDline in fileSDlines[1:]:
        fileSDlin = fileSDline.strip().split(" ")
        keySD = str(fileSDlin[0])
        print("keySD: "+str(keySD))        
        if keySD in dictSQ.keys():
            dictscP = ""
            for fileSDli in fileSDlin[1:]:
                dictscP += str(fileSDli) + "\t"
            outcon.write(str(keySD)+"\t"+str(dictSQ[keySD])+str(dictscP)+"\n")

outcon.close()