althonos / pyhmmer

Cython bindings and Python interface to HMMER3.
https://pyhmmer.readthedocs.io
MIT License
125 stars 12 forks source link

Fix transition probabilities produced by `Builder.build_msa` #14

Closed althonos closed 2 years ago

althonos commented 2 years ago

As reported via e-mail, the current code of Builder.build_msa was not producing HMMs consistent with hmmbuild.

The issue was occurring because the builder was always using gap-open and gap-extend probabilities to compute the transition probabilities instead of the alignment itself.