From 97b6b28d423443b6a703db81390f357644b06f67 Mon Sep 17 00:00:00 2001
From: Gabriel Studer <gabriel.studer@unibas.ch>
Date: Thu, 2 Jul 2020 16:03:40 +0200
Subject: [PATCH] avoid hardcoded evalue thresh when building hhblits profile
---
modules/bindings/pymod/hhblits.py | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/modules/bindings/pymod/hhblits.py b/modules/bindings/pymod/hhblits.py
index 861308e06..e8db627ad 100644
--- a/modules/bindings/pymod/hhblits.py
+++ b/modules/bindings/pymod/hhblits.py
@@ -569,11 +569,12 @@ class HHblits:
full_nrdb = os.path.join(os.path.abspath(os.path.split(nrdb)[0]),
os.path.split(nrdb)[1])
# create MSA
- opts = {'cpu' : 1, # no. of cpus used
- 'n' : 1} # no. of iterations
+ opts = {'cpu' : 1, # no. of cpus used
+ 'n' : 1, # no. of iterations
+ 'e' : 0.001} # evalue threshold
opts.update(options)
opt_cmd, _ = _ParseOptions(opts)
- hhblits_cmd = '%s -e 0.001 -i %s -oa3m %s -d %s %s' % \
+ hhblits_cmd = '%s -i %s -oa3m %s -d %s %s' % \
(self.hhblits_bin, self.filename, a3m_file, full_nrdb,
opt_cmd)
--
GitLab