Question: Psiblast Warning: Composition-Based Score Adjustment
4
gravatar for Stef
8.6 years ago by
Stef40
Stef40 wrote:

Warning: lcl|Query_1 1DCH:A|PDBID|CHAIN|SEQUENCE: Warning: Composition-based score adjustment conditioned on sequence properties and unconditional composition-based score adjustment is not supported with PSSMs, resetting to default value of standard composition-based statistics

Can someone help me as to what this warning message means? I mean the results that I get are pretty fine, so until now I just ignored it. But slowly I'm wondering if I'm doing a mistake.

I run PSI-Blast with following parameters

./psiblast -query 1cbyA.seq -db ./blastdb -matrix PAM30 -numiterations 5 -numalignments 3000 -num_descriptions 3000 -gapopen 32767 -gapextend 32767 -evalue 10 -outfmt 7 -out ./output.txt

blast statistics • 5.0k views
ADD COMMENTlink modified 2.5 years ago by ChuanD0 • written 8.6 years ago by Stef40

Hello,Stef!I'm facing the same problem as you met 3 years ago. Did you find an appropriate method to solve it?I just set the parameter -numiteration to be 2,but still the warning appears.

ADD REPLYlink written 5.5 years ago by crystalrose40
5
gravatar for Michael Kuhn
8.6 years ago by
Michael Kuhn5.0k
EMBL Heidelberg
Michael Kuhn5.0k wrote:

Composition-based score adjustment adapts the substitution matrix to the context of the sequences to increase performance, see Yu & Altschul, 2005. PSI-BLAST, on the other hand, generates profiles for each position, so there is no substitution matrix to adjust after the first iteration. Therefore, the warning makes sense, but you don't need to worry about it.

I'm not sure why the warning is triggered, your params seem fine (although I can't judge why you force an ungapped alignment, this seems to contradict the aim of finding more remote homologs with PSI-BLAST).

ADD COMMENTlink written 8.6 years ago by Michael Kuhn5.0k

Hm but I do give him as parameter a substitution matrix. If it's only used for the first iteration I guess in the second it throws the warning or? So that would be the reason(still kind of stupid to throw it though)

Well for my application I can only use ungapped alignments so I have to forfeit the advantages of finding remote homologs with gaps since the are useless to me. I mean if another application exists that can do fast sequence alignments without gaps I guess it would fit my needs better, but I don't know of any, so I just force BLAST and HHSearch to no gaps.

Thanks

ADD REPLYlink written 8.6 years ago by Stef40
0
gravatar for ChuanD
2.5 years ago by
ChuanD0
ChuanD0 wrote:

I meet the same problem. Did you solve it?

ADD COMMENTlink modified 2.5 years ago • written 2.5 years ago by ChuanD0
0
gravatar for ChuanD
2.5 years ago by
ChuanD0
ChuanD0 wrote:

I added "-comp_based_stats 1" option, which avoid to print the warning information.

ADD COMMENTlink written 2.5 years ago by ChuanD0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1719 users visited in the last hour