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

Peter biopython at maubp.freeserve.co.uk
Fri Aug 28 14:38:36 UTC 2009


bar tomas wrote:
> 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?

Yes, you can use a species specific blast database or include an
Entrez query, see for example:
http://lists.open-bio.org/pipermail/biopython/2009-June/005215.html

Peter




More information about the Biopython mailing list