Question: How to get the maximum number of hits in BLAST 2.2.31+?
1
gravatar for Neha shri
3.4 years ago by
Neha shri30
United States
Neha shri30 wrote:

Hi ,

                I have a locally created database which contains thousands of protein sequences . I want to check the positive substitutions in every sequence with reference to the target protein  as random substitutions were made in every sequence. The format specifier i am using is "positive -Number of positive-scoring matches".

The command line used-

Blastp  -query  -db  -outfmt  "7 ssequid  mismatch  positive  ppos "

Query is " Whenever i am running BLAST , it generates only 500 hits, which is may be the default, How can i get the output for all the sequences.

 I am new to using command line BLAST. Any help will really be appreciated.


 

blast alignment • 1.6k views
ADD COMMENTlink modified 3.4 years ago by Protostome50 • written 3.4 years ago by Neha shri30
2
gravatar for Protostome
3.4 years ago by
Protostome50
European Union
Protostome50 wrote:

Use either of the following flags (taken from blastp -help)

 -max_target_seqs <Integer, >=1>
   Maximum number of aligned sequences to keep 
   Not applicable for outfmt <= 4
   Default = `500'

 

OR:

 -num_descriptions <Integer, >=0>

   Number of database sequences to show one-line descriptions for
   Not applicable for outfmt > 4
   Default = `500'
    * Incompatible with:  max_target_seqs
 -num_alignments <Integer, >=0>
   Number of database sequences to show alignments for
   Default = `250'

 

ADD COMMENTlink written 3.4 years ago by Protostome50

Thank you so much :) I am able to do it now.

Could you please also tell me a way to get the full fasta sequences in the output instead of only their sequence ids?

 

 

ADD REPLYlink written 3.4 years ago by Neha shri30
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: 1990 users visited in the last hour