bugfix

modelling/src/afdb.cc

@@ -118,6 +118,7 @@ void SeqToPentamerIndices(const String& seq, bool unique, std::vector<int>& indi
   }
   indices.resize(N-4);
   if(unique) {
+    std::sort(indices.begin(), indices.end());
     auto last = std::unique(indices.begin(), indices.end());
     indices.erase(last, indices.end());
   }