Home
Reading
Searching
Subscribe
Sponsors
Statistics
Posting
Contact
Spam
Lists
Links
About
Hosting
Filtering
Features Download
Marketing
Archives
FAQ
Blog
 
Gmane
From: Jason Stajich <jason.stajich <at> gmail.com>
Subject: Re: Modifying Blast plus parameters
Newsgroups: gmane.comp.lang.perl.bio.general
Date: Tuesday 6th May 2014 17:37:10 UTC (over 2 years ago)
Or ggsearch or ssearch in fasta package


Jason Stajich
[email protected]


On Mon, May 5, 2014 at 5:58 AM, Paul Cantalupo 
wrote:

> 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  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
> > [email protected]
> > http://lists.open-bio.org/mailman/listinfo/bioperl-l
> >
> _______________________________________________
> Bioperl-l mailing list
> [email protected]
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
 
CD: 17ms