diff --git a/scripts/create_metagene_profile.py b/scripts/create_metagene_profile.py
index 6e95ce41d61690b3f449e0e6d15faec0c56084ea..7b56569f6a8895b9e62f5ececa29ce9b8d670997 100755
--- a/scripts/create_metagene_profile.py
+++ b/scripts/create_metagene_profile.py
@@ -77,7 +77,7 @@ from tempfile import NamedTemporaryFile
 from collections import defaultdict
 # from itertools import chain
 from sqlite3 import OperationalError
-from yaml import load as yload
+from yaml import safe_load as yload
 # https://pythonhosted.org/gffutils/
 from gffutils import FeatureDB, create_db
 # from deeptools import heatmapper
@@ -254,7 +254,7 @@ PARAMETER_INFO = "\n".join(DEFAULT_PARAMETERS.keys())
 #                    **extra_parameters):
 #     """Combine information from bigwig files *bigwig_filenames* and bed file
 #     *bed_filename*.
-# 
+#
 #     If *plot_filename* is set, write the corresponding meta profile
 #     in this file.
 #     """
@@ -279,7 +279,7 @@ PARAMETER_INFO = "\n".join(DEFAULT_PARAMETERS.keys())
 #                 val = [None] * nb_samples
 #         hm_params[key] = val
 #     heatm.parameters = hm_params
-# 
+#
 #     if plot_filename is not None:
 #         print(f"plotting profile to {plot_filename}")
 #         prof_params = keyfilter(is_prof_param, parameters)