[Biopython] (no subject)

Matthias Schade matthiasschade.de at googlemail.com
Sat Apr 14 11:41:33 UTC 2012


Hello everyone,


I would like to run a blastn-query of a small nucleotide-sequence against a genome. The code works already, but my queries are still slow and mostly ineffective, so I would like to ask:

Is there a way to tell the blastn-algorithm that once a 'perfect match' has been found it can stop and send back the results?

Background: I am interested in only the first full match because I would like to design a nucleotide-probe which -if possible- has no(!) known match in a host-genome, neither in RNA nor DNA. Actually, I would reject all perfect-matches and all single-mismatches but allow every sequence with two or more mismatches.

Currrently, I use this line of code with seq_now being about 15-30 nt long:
result_handle = NCBIWWW.qblast("blastn", "nt", seq_now, entrez_query="Canis familiaris[orgn]")


I am still new to this. Thank you for your help and input,

Matt




More information about the Biopython mailing list