diff --git a/traitar/hmmer2filtered_best.py b/traitar/hmmer2filtered_best.py
index 514b95f..9829e46 100755
--- a/traitar/hmmer2filtered_best.py
+++ b/traitar/hmmer2filtered_best.py
@@ -47,7 +47,7 @@ def aggregate_domain_hits(filtered_df, out_f):
     #sort by gene identifier and Pfam
     with open(out_f, 'w') as out_fo:
         ps.DataFrame(filtered_df.columns).T.to_csv(out_f, sep = "\t", index = False, header = False, mode = 'a')
-        filtered_df.sort(columns = ["target name", "query name"], inplace = True)
+        filtered_df.sort_values(by = ["target name", "query name"], inplace = True)
         if filtered_df.shape[0] > 0:
             current_max = filtered_df.iloc[0,] 
         else: