[Biopython] How to restrict blast query to proteins of a certain species

bar tomas bartomas at gmail.com
Fri Aug 28 11:27:20 EDT 2009


Hi,
I just tried using a taxonomic id:
The following both give the same result:


result_handle = NCBIWWW.qblast(*"blastp"*, *"nr"*,
fasta_string,entrez_query=*"house mouse[orgn]"*)

result_handle = NCBIWWW.qblast(*"blastp"*, *"nr"*,
fasta_string,entrez_query=*"txid10090[orgn]"*)

(entrez_query needs to be after the non keyword arguments or else the
compiler complains)




On Fri, Aug 28, 2009 at 4:19 PM, Kelly F Oakeson <kelly.oakeson at utah.edu>wrote:

> Hell list,
> I would like to do something similar, but I would like to limit my
> blast search to just the microbial Taxonomic ID.
>
> Thanks,
>
> Kelly Oakeson
> kelly.oakeson at utah.edu
>
>
>
>
>  On Aug 28, 2009, at 8:17 AM, bar tomas wrote:
>
> > Hi,
> >
> > I'd like to use blast to measure homology  between protein sequences
> > of a
> > species and polyketide sequences from a database.
> > I've been looking in the Biopython tutorial(p.73) in the section about
> > blast.
> > I'd like to do something similar,  like this:
> >
> > result_handle = NCBIWWW.qblast("blastp", "genbank", record.format
> > ("fasta"))
> >
> > Is it possible to add an option to restrict the search to genbank
> > records
> > that correspond to a given species?
> >
> > Thanks very much
>  > _______________________________________________
> > Biopython mailing list  -  Biopython at lists.open-bio.org
> > http://lists.open-bio.org/mailman/listinfo/biopython
>
>


More information about the Biopython mailing list