blast against genomes in biopython

Question

from Bio.Blast import NCBIXML
from Bio.Blast import NCBIWWW

result_handle = NCBIWWW.qblast(
    \"blastn\",
    \"nr\",
    \"CAC         


        

Answer 1

Your query isn't returning any results. The default parameters for blast are the cause. These parameters work better in this particular case of small length queries:

result_handle = NCBIWWW.qblast(
    "blastn",
    "nr",
    "CACTTATTTAGTTAGCTTGCAACCCTGGATTTTTGTTTACTGGAGAGGCC",
    megablast=False,
    expect=1000,
    word_size=7,
    nucl_reward=1,
    nucl_penalty=-3,
    gapcosts="5 2",
    entrez_query='Beutenbergia cavernae DSM 12333 [Organism]')

Particularly the expect parameter plays a major role here.