[Biopython] Biopython local blastn query

Peter Cock p.j.a.cock at googlemail.com
Tue Jul 30 17:08:16 UTC 2013


On Tue, Jul 30, 2013 at 6:02 PM, Ara Kooser <ghashsnaga at gmail.com> wrote:
> This will sound like a silly question. I found the nt.nal file that lists
> all the databses. How do I call the alias from biopython?
>
> I thought it would be something like this:
>
> nt = "/Users/arakooser/blast/db/nt.nal"
>
>  result = NcbiblastnCommandline(task="megablast",query="-", db=nt,
>                                    outfmt=5, perc_identity=100,
> out="temp.xml",
>                                    max_hsps_per_subject=1, num_alignments=1)
>
> But that throws an error letting me know that nothing was returned.
>
> ara

Just as a string in quotes, "nt",

NcbiblastnCommandline(task="megablast", query="-", db="nt", ...)

Peter



More information about the Biopython mailing list