bix at sendu.me.uk
Wed Aug 9 04:43:58 EDT 2006
Bernd Brandt wrote:
> "frac_identical and frac_conserved can still be very wrong"
> They were wrong with hmm report parsing (hmmsearch). The bioperl-live
> (CVS 14 July 2006) returned 0 for both fracions. I will check it with
> the newest CVS and send a small test script.
For both hmmsearch and hmmpfam parsing, those are deliberately set to 0.
They needn't be; there is a midline to the hsp-like-thing that shows
identical and conserved positions. I suppose the reason they're set to 0
currently is that its not a real alignment against a set sequence, but
against a model. Still, I think it's reasonable to give an answer in
both cases if you can get an answer from the original output files.
I'm in the process of writing new hmmer parsers, and will allow non-0
answers for those methods.
(Its not really possible to alter the existing parser because it only
allows one hsp per hit, preventing a correct answer even if it didn't
set to 0)
More information about the Bioperl-l