diff --git a/modules/io/pymod/export_omf_io.cc b/modules/io/pymod/export_omf_io.cc
index fb98237b0ffbb2e5c6e5cd6d00ceceeefde080dd..43355a4a908f97ae3f567daba57f5f0e739c34c5 100644
--- a/modules/io/pymod/export_omf_io.cc
+++ b/modules/io/pymod/export_omf_io.cc
@@ -77,6 +77,7 @@ void export_omf_io() {
     .def("GetChainNames", &wrap_get_chain_names)
     .def("GetPositions", &OMF::GetPositions, return_value_policy<reference_existing_object>(),(arg("cname")))
     .def("GetBFactors", &OMF::GetBFactors, return_value_policy<reference_existing_object>(),(arg("cname")))
+    .def("GetAvgBFactors", &OMF::GetAvgBFactors,(arg("cname")))
     .def("GetSequence", &OMF::GetSequence, (arg("cname")))
   ;
 }
diff --git a/modules/io/src/mol/omf.cc b/modules/io/src/mol/omf.cc
index 973595161ed6295631bf6e26593d9984b81b7d1e..91fdca413267c42bf41582e4422a37149bdbb59c 100644
--- a/modules/io/src/mol/omf.cc
+++ b/modules/io/src/mol/omf.cc
@@ -3833,6 +3833,31 @@ const std::vector<Real>& OMF::GetBFactors(const String& cname) const {
   return it->second->bfactors;
 }
 
+std::vector<Real> OMF::GetAvgBFactors(const String& cname) const {
+  auto it = chain_data_.find(cname);
+  if(it == chain_data_.end()) {
+    throw ost::Error("Provided chain name not in OMF structure");
+  }
+  const std::vector<Real>& bfactors = it->second->bfactors;
+  const std::vector<int>& res_def_indices = it->second->res_def_indices;
+  std::vector<Real> avg_bfactors;
+  avg_bfactors.reserve(it->second->res_def_indices.size());
+  int current_atom_idx = 0;
+  for(auto i = res_def_indices.begin(); i != res_def_indices.end(); ++i) {
+    int size = residue_definitions_[*i].anames.size();
+    Real summed_bfac = 0.0;
+    for(int j = 0; j < size; ++j) {
+      summed_bfac += bfactors[current_atom_idx];
+      ++current_atom_idx;
+    }
+    if(size > 0) {
+      summed_bfac /= size;
+    }
+    avg_bfactors.push_back(summed_bfac);
+  }
+  return avg_bfactors;
+}
+
 String OMF::GetSequence(const String& cname) const {
   auto it = chain_data_.find(cname);
   if(it == chain_data_.end()) {
diff --git a/modules/io/src/mol/omf.hh b/modules/io/src/mol/omf.hh
index f5bb8b9966b0b5c74d9d114f745f212c16b52799..9503b1f541a928b1a3f3730a39fffc60bbab36c4 100644
--- a/modules/io/src/mol/omf.hh
+++ b/modules/io/src/mol/omf.hh
@@ -199,6 +199,8 @@ public:
 
   const std::vector<Real>& GetBFactors(const String& cname) const;
 
+  std::vector<Real> GetAvgBFactors(const String& cname) const;
+
   String GetSequence(const String& cname) const;
 
 private: