Invalid Api For Biopython'S Blast?

According to BioPython's API, Bio.Blast.NCBIWWW.qblast() should have a parameter of specifying megablast to be true or false. When I set my qblast as

result_handle = NCBIWWW.qblast("blastn", "nr", fasta_string, megablast=True)

then I get the error of

TypeError: qblast() got an unexpected keyword argument 'megablast'

What gives? Am I doing something wrong, or does BioPython not actually have that option?

biopython

• 3,142 views • link

updated 14.8 years ago by Peter

• written 14.8 years ago by Jesse J

Are you using the latest version of Biopython?

• 0 views • link

written 14.8 years ago by Niek De Klein

1 answer

The megablast argument was added to the qblast function in October 2010, so you need at least Biopython 1.56 which was released in November 2010. Please upgrade, the current release is Biopython 1.58 released August 2011.

• 2 views • link

written 14.8 years ago by Peter

Huh...turns out I do have an older version. Which is weird considering I installed Biopython only a few weeks ago.

• 0 views • link

written 14.8 years ago by Jesse J

If you got it on Linux via the distribution's package manager, they often don't have the latest release.

• 0 views • link

written 14.8 years ago by Peter

More posts like this

Can't get results using Biopyhon BLAST
written 4.6 years ago by Ric95 •

Hi everyone, I am writing a script to BLAST sequences using Biopython. Following previous posts like this, I tried to put all the constraints used …
How can I run a Biopython qblast search targeting a specific set of organisms?
written 6.2 years ago by sviatoslav.kendall

I want to use Biopython's qblast() command to query NCBI's databases but I want to limit my search to specific organisms. Looking at the [documentation …
How to write a Genbank Reference with Biopythons SeqIO?
written 7.4 years ago by Mathis DKZ •

Hello, I'm trying to create a Genbank file with the SeqIO Feature of Biopython. I am able to parse references of existing files with: for …
Making Annotated Chromosome Using Biopython
written 12.4 years ago by Janake

I am trying to recreate the annotated chromosome using biopython (<a href='http://biopython.org/DIST/docs/tutorial/Tutorial.html#sec345'>http://biopython.org/DIST/docs/tutorial/Tutorial.html#sec345</a>). I have a test code that would create one chromosome and one annotated …
Remove Uncultured Organisms From Blast Results With Biopython
written 13.2 years ago by t.poisot •

Hi, I'm trying to use biopython to blast a list of sequences and retrieve the best hit. This works fine. Now I'd like to filter …
How To Calculate Number Of Gap Windows And Their Length In A Multiple Alignment
written 13.9 years ago by Biojl

Hi, Anyone knows if there is an implemented function in biopython to calculate gap windows and their length in a multiple sequence alignment. I found …
Ncbi Blast Failing For Unknown Reasons
written 14.8 years ago by Jesse J

I have a Python script that runs BioPython's Web Blast function. We're using large fasta files, so the script breaks these files into smaller files …
Adding Consensus Numbering To An Antibody Amino Acid Sequence In Biopython
written 14.8 years ago by Jared Sampson •

I would like to add the consensus numbering, as determined by Andrew Smith's <a href="http://www.bioinf.org.uk/abs/abnum/">Abnum</a> utility, to a <a href="http://biopython.org/wiki/Main_Page">Biopython</a> Seq object. After browsing the …
How To Parse Psiblast Results Using Biopython And Blast-2.2.24+?
written 15.9 years ago by Niek De Klein

I'm trying to run a PSIBlast program which selects certain sequences out at every round before it does the next iteration. For this I need …
How To Show More Blast Results Using Biopython?
written 15.9 years ago by Niek De Klein

Hi, I'm using biopython to BLAST over the internet. However, it only saves 30 results (there are more than 30 results that are under the …