Yes thank you! I've been looking for a site like that.
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 e-value I chose) in the xml. I've been looking all over but can't find how to make that number higher. So my question is, how can you show more results from BLAST using biopython. I'm using NCBIWWW.qblast from BIO.BLAST.
from Bio.Blast import NCBIWWW
File = "MIF"
fasta_string = open(File+".fasta").read()
result_handle = NCBIWWW.qblast("blastp", "nr", fasta_string)
Thanks, Niek
3 answers
Hello Niek,
Daniel pointed out an email to the Biopython mailing list which answers your question - basically try reading the built in documentation. This is also mentioned in the Biopython Tutorial, quote:
For more about the optional BLAST arguments, we refer you to the NCBI’s own documentation, or that built into Biopython:
from Bio.Blast import NCBIWWW
help(NCBIWWW.qblast)
Also audyyy pointed out some old (but apparently still relevant) documentation from Andrew Dalke's homepage.
You said
I've been looking all over but can't find how to make that number higher.
so you probably looked at some of these bits of documentation (right?). Could you suggest how it could be made clearer for newcomers? Feedback does help to improve the docs.
Thanks,
Peter
Log in to answer this question.
Which blast service are you using?
NCBIWWW.qblast
File = "MIF" fasta_string = open(File+".fasta").read() result_handle = NCBIWWW.qblast("blastp", "nr", fasta_string)
NCBIWWW.qblast
NCBIWWW.qblast
I editted it in my post.