I'm trying to search homologs for a query sequence in a list of genomes given their accesion numbers obtained from ref seq nucleotides database. I have no results and I don't understand how to find the homologs of the protein in a genome. Please help me! Thank you.

The next one is the code:

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

# Read or define your query protein sequence
query_sequence = "MAAGVK"

# Define the genome accession numbers
genome_accession_numbers = ["NC_000913.3", "NC_060200.1"]  # Example accession numbers for E. coli and S. cerevisiae

# Perform BLASTP search against each genome
for accession in genome_accession_numbers:
    # Build the entrez query for the specific genome
    entrez_query = f"refseq_genomic[filter] AND {accession}[accession]"

    # Perform BLASTP search for the specific genome
    result_handle = NCBIWWW.qblast("blastp", "refseq_select", query_sequence, entrez_query=entrez_query)

    # Parse the results
    blast_records = NCBIXML.parse(result_handle)

    # Process each record
    for record in blast_records:
        for alignment in record.alignments:
            for hsp in alignment.hsps:
                # Access relevant information
                print("Hit ID:", alignment.hit_id)
                print("Hit Description:", alignment.hit_def)
                print("Alignment Length:", hsp.align_length)
                print("Bit Score:", hsp.bits)
                print("E-Value:", hsp.expect)
                print("Alignment Sequence:", hsp.query)
                print("Hit Sequence:", hsp.sbjct)
                print("\n")

# Close the result handle
result_handle.close()