[Bioperl-l] Modifying Blast plus parameters
Paul Cantalupo
pcantalupo at gmail.com
Mon May 5 12:58:15 UTC 2014
Hi Oz,
SInce BLAST is a local alignment tool, perhaps you want something like
Needle http://www.ebi.ac.uk/Tools/psa/emboss_needle/nucleotide.html which
performs a global alignment?
You can add the '-perc_identity' option to the blastn command line like so:
blastn -perc_identity 80
to limit results to >=80% identity.
Paul
Paul Cantalupo
University of Pittsburgh
On Tue, Apr 29, 2014 at 12:49 PM, Oz <oz.sharabi at mail.huji.ac.il> wrote:
> Hi Paul,
> Thank you for your quick and informative reply.
> Yes, it is a local blast. About sequence coverage - I meant that I do not
> want the algorithm to split the query to two sub-sequences (or more) and
> then align each one at different location. Just to align it as one
> continous sequence.
> Also, I would like to make these restrictions during the Blastplus running
> and not to screen the un-wanted results afterwards..
> Since I'm a "green" user, could you please specify where I add the
> -perc_identity (I couldn't find example on the web)?
> Thank you very much,
> Oz
>
>
>
>
> --
> View this message in context:
> http://bioperl.996286.n3.nabble.com/Modifying-Blast-plus-parameters-tp17461p17478.html
> Sent from the Bioperl-L mailing list archive at Nabble.com.
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
More information about the Bioperl-l
mailing list