I'm looking to train a profile HMM to try and recognize different domains in a protean sequence. I'm sure there are library's to do this and was wondering what the best one to use would be. I'd prefer something in python or with python bindings. Thanks