I am trying to get a list of proteins in a genomic region surrounding a protein of interest identified by a protein accession number. I am using biopython and the entrez module. First I start off by getting a list of protein accession from GenBank and now I want to get an idea of the genomic region surrounding the protein.
I am using
Entrez.efetch(db="protein", id=rec, rettype="ipg", retmode="text")
In order to get the nucleotide accession number and the start / stop sites of the protein CDS.
My question is, how do I then download a genbank file representing a 30kb region around the CDS.
Can someone point me in the right direction.