--- a +++ b/retrieve_genome_sequence.py @@ -0,0 +1,17 @@ +from Bio import Entrez + +def retrieve_genome_sequence(genome_id): + Entrez.email = 'your_email@example.com' # Set your email address here + handle = Entrez.efetch(db='nucleotide', id=genome_id, rettype='fasta', retmode='text') + record = handle.read() + handle.close() + return record + +# Example usage +genome_id = 'NC_045512' # Replace with the actual genome ID +sequence = retrieve_genome_sequence(genome_id) + +# Output the DNA sequence in a markdown code block +print("```") +print(sequence) +print("```")