https://github.com/fjruizruano/ngs-protocols
Raw File
Tip revision: 39a091d1fa569a7fc717ac73c4b3de07f0a1204d authored by fjruizruano on 03 August 2023, 11:48:27 UTC
adding gfa2fas.py and extract_gfa.py
Tip revision: 39a091d
join_rm_list.py
#!/usr/bin/python

import sys

file1 = open(sys.argv[1]).readlines()
file2 = open(sys.argv[2]).readlines()

dict1 = {}
dict2 = {}
dicta = {}

for el in file1:
	sp = el.split("\t")
	dict1[sp[0]] = sp[1][:-1]

for el in file2:
	sp = el.split("\t")
	dict2[sp[0]] = sp[1][:-1]

elements = list(set(dict1.keys() + dict2.keys()))

for el in elements:
	dicta[el] = []
	if el in dict1:
		dicta[el].append(dict1[el])
	else:
		dicta[el].append("NA NA")

	if el in dict2:
		dicta[el].append(dict2[el])
	else:
		dicta[el].append("NA NA")

dicta2 = sorted(dicta.iterkeys())

out = open("salida.txt", "w")
out.write("Name\t%s \t%s \n" % (sys.argv[1], sys.argv[2]))

for el in dicta2:
	out.write("%s\t%s\t%s\n" % (el, dicta[el][0], dicta[el][1]))

out.close()
back to top