Suppose you've got a local BLAST (specifically, blastp) instance and a "fasta reference dataset", and you just wanna align (using blastp) a given query protein sequence against sequences from that dataset whose lengths are the same (or higher) than query's length. Can you set any option when running blastp so that it performs this behavior? If so, which option would be, and how to activate it?
•
written
by
inaciogmedeiros •
1 answer
Not exactly what you asked for, but UniProt's BLAST (http://www.uniprot.org/blast) allows you to filter your results in a post-processing step. You can "map to UniProtKB" and then add any search criteria, e.g. sequence length (range), by using the advanced search. However, this just allows you to retrieve/see the sequences, and doesn't maintain the scoring/ranking by BLAST similarity score/e-value etc.
Log in to answer this question.
-
Changing Temporary File Location for Local BLASTP to nr Database
written by katieostrouchov •
-
BLAST e-value between two sequences vs BLAST e-value between sequence and data base
written by nyanovsky •
-
What should be the best parameter to perform blastp
written by Nelo •
-
conducting individual blastp for multiple protein sequences
written by rb77 •
-
BLASTP local vs website
written by timothy.kirkwood •
-
Changing to new BLAST+ local program - database error
written by daisyrowley •
-
Set up a blast remote server
written by ostrokach
-
standalone blastp: increasing word size extremely slows down the search
written by aleksanderczeszyk •
-
Blastp command line to get only full matches of query sequence
written by rlee03 •
-
Get Query Coverage From Local Blast
written by bioinf80 •
why don't you just filter before building the blastdb or filter the blast output ?
I thought this too, and wondered if BLAST could do this "automatically" for me.
Include query length, subject length and alignment length in tabular output and you can do it post blastp with awk with very little effort.
blastp -helpfor how to include this stuff in blast output..