Skip to content
Snippets Groups Projects
Commit b958bfac authored by Riccardo Pellarin's avatar Riccardo Pellarin
Browse files

print colored positions to check consistency of residue mapping

parent e6f3cfde
No related branches found
No related tags found
No related merge requests found
File added
...@@ -27,7 +27,23 @@ for csv in ["Resultat.b.csv"]: ...@@ -27,7 +27,23 @@ for csv in ["Resultat.b.csv"]:
#cldb_tmp.rename_proteins({"TssKHis6":"TssK","StrepTssF":"TssF","TssGFlag":"TssG","HATssE":"TssE"}) #cldb_tmp.rename_proteins({"TssKHis6":"TssK","StrepTssF":"TssF","TssGFlag":"TssG","HATssE":"TssE"})
cldb_tmp.set_name(csv) cldb_tmp.set_name(csv)
cldb.append_database(cldb_tmp) cldb.append_database(cldb_tmp)
cldb.align_sequence('pairwise_alignment_NeisseriaTo_3ja1n_0.08.txt') cldb.align_sequence('tools/pairwise_alignment_NeisseriaTo_3ja1n_0.08.txt')
renamedict={}
for xl in cldb:
try:
protname=xl[cldb.protein1_key].split("|")[2]
renamedict[xl[cldb.protein1_key]]=protname
except:
pass
try:
protname=xl[cldb.protein2_key].split("|")[2]
renamedict[xl[cldb.protein2_key]]=protname
except:
pass
cldb.rename_proteins(renamedict)
#print(cldb) #print(cldb)
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment