diff --git a/modelling/pymod/_afdb_modelling.py b/modelling/pymod/_afdb_modelling.py
index 08073ee49a67511b53aeff8b61cf705df1fb91e9..243abc35bab7c02ed2029f2f6e51a7ad49a6d4e4 100644
--- a/modelling/pymod/_afdb_modelling.py
+++ b/modelling/pymod/_afdb_modelling.py
@@ -548,7 +548,7 @@ def AFDBTPLSearch(fs_server, pentamatch, trg_seq, pentamatch_n = 100,
         aln = seq.alg.LocalAlign(trg_seq, seq.CreateSequence("A", omf_s),
                                  seq.alg.BLOSUM62)[0]
         if seq.alg.SequenceIdentity(aln) >= seqid_thresh:
-            bfactors = omf.GetBFactors("A")
+            bfactors = omf.GetAvgBFactors("A")
             summed_bfac = 0.0
             current_pos = aln.GetSequence(1).offset
             for col in aln: