Question: stand alone psi-blast
0
gravatar for crystalrose
5.6 years ago by
crystalrose40
China
crystalrose40 wrote:

I'm using psi-blast (blast version 2.2.30+),this is my command:

psiblast -query seq.fasta -db db -num_threads 15 -out outfile.txt -outfmt 6

I'm not sure than if I have got the right result.It seems that it only run one iteration. But if I add the parameter: -num_interations 2, warnings appear: 

Frequency ratios for PSSM are all zeros, frequency ratios for BLOSUM62 will be used during traceback in composition based statistics

I don't know where is the mistake and how can PSSM are all zeros. any help will be appreciated!

ADD COMMENTlink modified 3.2 years ago by 1992minto0 • written 5.6 years ago by crystalrose40
0
gravatar for RamRS
5.6 years ago by
RamRS28k
Houston, TX
RamRS28k wrote:

PSSM scores might be all zeros because no domains were found in the input sequences. Could you give us more details on how the db was created?

ADD COMMENTlink written 5.6 years ago by RamRS28k

Thanks for your answer.I create the db using my own fasta file:mekeblastdb -in seq.fasta -dbtype prot -out db -parse_seqids.The fasta file here is the same as psiblast query file.I want to classify these protein seqeuences ,which can't be classified into pfam ,by psiblast.

ADD REPLYlink written 5.6 years ago by crystalrose40

Are you sure running PSI-BLAST against the same set of sequences as both subject and query is a good idea? What is your ultimate aim?

ADD REPLYlink written 5.6 years ago by RamRS28k

Actually , I aim to classify protein families for some protein suquences.I think I should change my strategy and  search the known-function-sequence db. Do you think this can be work? 

ADD REPLYlink written 5.6 years ago by crystalrose40

You might wanna just query against nr or SwissProt. Classification is usually done by checking homology to proteins with known function. If your query is all unknowns and you're not looking to internally cluster them, the DB should be nr or SwissProt.

ADD REPLYlink written 5.6 years ago by RamRS28k

Thanks for your help!I sovled this problem by adding an parameter "-comp_based_stats 0" and got some conserved results. 

ADD REPLYlink written 5.6 years ago by crystalrose40
0
gravatar for 1992minto
3.2 years ago by
1992minto0
1992minto0 wrote:

enter image description herehi ! here , I get the two matrix of PSSM , but I don't their meaning .As shown in the follow picture

enter image description here

ADD COMMENTlink modified 3.2 years ago • written 3.2 years ago by 1992minto0

You have broken links. Please open a new question and add a link to this question in there so people can help you out.

ADD REPLYlink written 3.2 years ago by RamRS28k
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: 1667 users visited in the last hour