Skip to content
Snippets Groups Projects
Commit df439707 authored by juergen's avatar juergen
Browse files

optimising handling properties of sequences in clustalw wrapper

git-svn-id: https://dng.biozentrum.unibas.ch/svn/openstructure/trunk@2171 5a81b35b-ba03-0410-adc8-b2c5c5119f08
parent 9a68affc
No related branches found
No related tags found
No related merge requests found
...@@ -16,17 +16,14 @@ def ClustalW(seq1, seq2, clustalw=None): ...@@ -16,17 +16,14 @@ def ClustalW(seq1, seq2, clustalw=None):
out) out)
ps=subprocess.Popen(command, shell=True, stdout=subprocess.PIPE) ps=subprocess.Popen(command, shell=True, stdout=subprocess.PIPE)
ps.stdout.readlines() ps.stdout.readlines()
taln=io.LoadAlignment(out) aln=io.LoadAlignment(out)
aln=seq.CreateAlignment() aln.SetSequenceOffset(0,seq1.GetSequenceOffset())
if (seq1.HasAttachedView() and seq2.HasAttachedView()): aln.SetSequenceOffset(1,seq2.GetSequenceOffset())
seq00=taln.GetSequence(0).Copy() if seq1.HasAttachedView():
seq00.AttachView(seq1.GetAttachedView().Copy()) aln.AttachView(0,seq1.GetAttachedView().Copy())
seq01=taln.GetSequence(1).Copy() if seq2.HasAttachedView():
seq01.AttachView(seq2.GetAttachedView().Copy()) aln.AttachView(1,seq2.GetAttachedView().Copy())
aln.AddSequence(seq00)
aln.AddSequence(seq01)
else:
aln=taln.Copy()
#~ temp_dir.Cleanup() #~ temp_dir.Cleanup()
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment