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 in the NCBI website, still can't get the same results... Can someone help me?

The code I'm using is the following.

from Bio.Blast import NCBIWWW
from Bio.Blast import NCBIXML
record = '''GGGGCCCTTCTTCTAGATTACCCATGCATGTTATTGAATATAAAGAAAAGACGAACAGCTCAGTAACTCT
TAGTGTAAAGATGTATTAGGACTGTTAACTATGTTAAATACAGATAATAATAATACAATAAGAATAAGAT
CTATCAGTTAGTTGTTAAGGTAATGGCTTAACAAGACAGTGACGGATAACGGTATTACTTTGTAATATTC
CGGAGAAGGAGCCTGAGAGATTGCTACTAAGTCTAAGGATTGCAGCAGGGGCTAAACTTGACCTATGATA
TATTATTGAGGCAGTTATGAGTAGTATTTTATAATTATTGTAATATTGTAAGTATATATTACAAGATGAA
TCGGAGGGCAAATCGAGTGCCAACAGCCGGGGGAATGCGGGGGGCGATAGTGTGTGTGATGATTGATGGA
GGGTGGGAGTCTGGAGGAGGAATGAGGATAAGGATGGGAGGGTATATTGGGGGGTTAGGAGAGAGACGAA
ATGTGATAACCCTAACTGGATGAACCGAAGCGAAAGCTATATTCTTTATTGTTTTATTAGAACAAGGTCG
TTTGTTTTGAGGTTTGAAGATGCTCAGATTCCTC'''
result_handle = NCBIWWW.qblast("blastn", "nr/nt", record,
                               megablast=False,
                               expect=10,
                               word_size=11,
                               nucl_reward=2,
                               nucl_penalty=-3,
                               gapcosts="5 2",
                               short_query=True,
                               entrez_query='microsporidians (taxid:6029)')
for res in NCBIXML.parse(result_handle):
    for d in res.descriptions:
        print(d.title,d.score,d.e)