Different conservation result for the protein P17931

[FlorenceChu]

Describe the bug
when i use the conservation_hmm_based.py to calculate the ic of protein P17931, the results of the residues whose value are -1000 in prankweb, are not negative but some normal value.

To Reproduce
command:
python conservation_hmm_based.py AF-P17931-F1-model_v4.fasta /hmmer-3.3.2/uniref50.fasta . AF-P17931-F1-model_v4.txt

Expected behavior
residue1-116 should be -1000.00 (in the prankweb)

Screenshots
the result i got:

[skodapetr]

@FlorenceChu Can you please specify which version/branch are you using to run the script?

[FlorenceChu]

Sorry I have just noticed this email. I found that it was my fault. When I set mask_output True, this problem is solved.  But I meet another problem. When I use a protein with multiple chains, there will be an error (attached photo). Please tell me what should I do with the multiple-chain protein? Thank you very much. And the code I’m using is https://github.com/cusbg/prankweb/blob/main/conservation/hmm_based/conservation_hmm_based.py 发自我的iPhone

…

------------------ Original ------------------ From: Petr Škoda ***@***.***> Date: Tue,Dec 31,2024 2:16 AM To: cusbg/prankweb ***@***.***> Cc: FlorenceChu ***@***.***>, Mention ***@***.***> Subject: Re: [cusbg/prankweb] Different conservation result for the proteinP17931 (Issue #152) @FlorenceChu Can you please specify which version/branch are you using to run the script? — Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you were mentioned.Message ID: ***@***.***>

[FlorenceChu]

This is the error report