[Biopython] qblast against RefSeq

Wibowo Arindrarto w.arindrarto at gmail.com
Sun Dec 9 04:59:07 UTC 2012


Hi Miguel,

If I recall correctly, you need to specify which 'refseq' database
you're using (there are 'refseq_rna' and 'refseq_protein'). Try
running the search again but using 'refseq_protein' as the database
name.

Hope that helps :),
Bow

On Sun, Dec 9, 2012 at 4:53 AM, Miguel Romero <phyrexian.kavu at gmail.com> wrote:
> Hi!
> Im learning how to run BLAST over the internet and Im interested in
> blasting against the RefSeq database. For now I have managed to blast
> against the nr database with:
>
> result_handle = NCBIWWW.qblast("blastp", "nr", sequence_string)
>
> but I dont know how to search against RefSeq, if I enter:
>
> result_handle = NCBIWWW.qblast("blastp", "refseq", sequence_string)
>
> the blast record is empty! Is it possible to blast against the RefSeq
> database? What should I type in the second argument?
>
> Thank you!
>
> Miguel Romero
>
> --
> [Theropoda is my profession]
> _______________________________________________
> Biopython mailing list  -  Biopython at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/biopython



More information about the Biopython mailing list