Skip to content
Snippets Groups Projects
Commit d1818e3b authored by Tobias Schmidt's avatar Tobias Schmidt
Browse files

added nopgap flag to clustalw binding

parent d711f981
Branches
Tags
No related merge requests found
...@@ -3,7 +3,7 @@ from ost import settings, io, seq, LogError ...@@ -3,7 +3,7 @@ from ost import settings, io, seq, LogError
import os import os
import subprocess import subprocess
def ClustalW(seq1, seq2=None, clustalw=None, keep_files=False): def ClustalW(seq1, seq2=None, clustalw=None, keep_files=False, nopgap=False):
clustalw_path=settings.Locate(('clustalw', 'clustalw2'), clustalw_path=settings.Locate(('clustalw', 'clustalw2'),
explicit_file_name=clustalw) explicit_file_name=clustalw)
...@@ -32,6 +32,8 @@ def ClustalW(seq1, seq2=None, clustalw=None, keep_files=False): ...@@ -32,6 +32,8 @@ def ClustalW(seq1, seq2=None, clustalw=None, keep_files=False):
command='%s -infile="%s" -output=fasta -outfile="%s"' % (clustalw_path, command='%s -infile="%s" -output=fasta -outfile="%s"' % (clustalw_path,
temp_dir.files[0], temp_dir.files[0],
out) out)
if nopgap:
command+=" -nopgap"
ps=subprocess.Popen(command, shell=True, stdout=subprocess.PIPE) ps=subprocess.Popen(command, shell=True, stdout=subprocess.PIPE)
ps.stdout.readlines() ps.stdout.readlines()
aln=io.LoadAlignment(out) aln=io.LoadAlignment(out)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment