diff --git a/src/dft/dft.cpp b/src/dft/dft.cpp
index 4c2ca64..587d8c3 100644
--- a/src/dft/dft.cpp
+++ b/src/dft/dft.cpp
@@ -3,16 +3,21 @@
 #include "storm-dft/builder/ExplicitDFTModelBuilderApprox.h"
 #include "storm/settings/SettingsManager.h"
 #include "storm-dft/settings/modules/DFTSettings.h"
+#include "storm-dft/storage/dft/DFTIsomorphism.h"
 
 // Thin wrapper for model building using one formula as argument
 template<typename ValueType>
-std::shared_ptr<storm::models::sparse::Model<ValueType>> buildModelFromJsonDft(std::string const& jsonDft) {
+std::shared_ptr<storm::models::sparse::Model<ValueType>> buildModelFromJsonDft(std::string const& jsonDft, bool symred) {
     // Build DFT
     storm::parser::DFTJsonParser<ValueType> parser;
     storm::storage::DFT<ValueType> dft = parser.parseJson(jsonDft);
     // Build model
     std::map<size_t, std::vector<std::vector<size_t>>> emptySymmetry;
     storm::storage::DFTIndependentSymmetries symmetries(emptySymmetry);
+    if (symred) {
+        auto colouring = dft.colourDFT();
+        symmetries = dft.findSymmetries(colouring);
+    }
     storm::builder::ExplicitDFTModelBuilderApprox<ValueType> builder(dft, symmetries, true);
     typename storm::builder::ExplicitDFTModelBuilderApprox<ValueType>::LabelOptions labeloptions({}, true);
     builder.buildModel(labeloptions, 0, 0.0);
@@ -24,6 +29,6 @@ void define_dft(py::module& m) {
             storm::settings::addModule<storm::settings::modules::DFTSettings>();
         }, "Initialize Storm-dft");
     // Build model
-    m.def("build_sparse_model_from_json_dft", &buildModelFromJsonDft<double>, "Build the model", py::arg("jsonDft"));
-    m.def("build_sparse_parametric_model_from_json_dft", &buildModelFromJsonDft<storm::RationalFunction>, "Build the parametric model", py::arg("jsonDft"));
+    m.def("build_sparse_model_from_json_dft", &buildModelFromJsonDft<double>, "Build the model", py::arg("jsonDft"), py::arg("symred") = false);
+    m.def("build_sparse_parametric_model_from_json_dft", &buildModelFromJsonDft<storm::RationalFunction>, "Build the parametric model", py::arg("jsonDft"), py::arg("symred") = false);
 }