Dear Modeller users,
I'm trying to force a local alignment while using salign function in mod9.15, and it fails with error:
"modeller.ModellerError: align_l_251E> MATRIX_OFFSET probably has an invalid value: 0.0000"
My script ( mod-3nnt-prim.cons-aln.py ):
from modeller import * env = environ() aln = alignment(env) m = model(env, file='./tmp-NmArZQ.pdb', model_segment=('FIRST:A', 'LAST:B')) aln.append_model(m, atom_files='./tmp-NmArZQ.pdb', align_codes='3nnt_x') aln.append( file='./tmp-q4VfPF.fasta', alignment_format='FASTA', align_codes='prim.cons') aln.salign(local_alignment=True) aln.write(file='mod-3nnt-prim.cons.pir')
The output:
Traceback (most recent call last): File "mod-3nnt-prim.cons-aln.py", line 8, in ? aln.salign(local_alignment=True) File "/usr/lib/modeller9.15/modlib/modeller/alignment.py", line 395, in salign nsegm, matrix_offset_3d, break_break_bonus) _modeller.ModellerError: align_l_251E> MATRIX_OFFSET probably has an invalid value: 0.0000
The error comes from a call to a function:
retval = func(self.modpt, ... , local_alignment, .... , nsegm, matrix_offset_3d, break_break_bonus)
that is defined as:
func = _modeller.mod_salign
If I remove the option "local_alignment=True", it runs flawlessly.
Any idea about what can be happening?
Thanks in advance,
David.