[Biopython-dev] GI number for NcbiblastpCommandline

[Brown, Tom]

Currently using result_handle = NCBIWWW.qblast("blastp", "nr", blastGI) where blastGI = 113000 in the Biopython program and would like to convert it to a local blastp. Is there a way to specify the blastGi within NcbiblastpCommandline instead of having to provide a fasta file for blast?  What are my options.

from Bio.Blast.Applications import NcbiblastpCommandline
blastp_cline = NcbiblastpCommandline(query="sh3.fasta", db="nr", evalue=0.001, outfmt=5, out="sh3.xml")
stdout, stderr = blastp_cline()

Thanks

Tom


________________________________




Confidentiality Notice:

This e-mail, including all attachments, is for the sole use of the intended recipient(s) and may contain confidential and privileged information as defined under FERPA. Any unauthorized review, use, disclosure or distribution is prohibited unless specifically provided under the New Mexico Inspection of Public Records Act. If you are not the intended recipient, please contact the sender and destroy all copies of this message