[Biopython] Biopython local blastn query

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


On Tue, Jul 30, 2013 at 5:10 PM, Ara Kooser <ghashsnaga at gmail.com> wrote:
> Peter,
>
>   Thanks for catching that! I missed that one. I also needed to upgrade to
> biopython 1.62b which I did.

Really? Maybe there was a BLAST wrapper update or something relevant?

> I still get one short sequence coming through.
>

BLAST e-value thresholds are not always the best approach to filtering...

> *General question*
> Hopefully one last question from me on this project. Can I query multiple
> blast databased in a single command? I have all the nt.xx downloaded and
> need to query each one to look for all my sequences.

There should be an nt.nal alias file so that you can just use "nt" as
the database name to search all of it.

Peter



More information about the Biopython mailing list