I have a list of GI accession numbers that I obtained following a protein BLAST search. These are from different bacterial genomes. I've pasted some examples below.
108765232 108769511 108769512 108769513 108769516 108799092 108799093 108799094 108799097 108803989
For further analysis, I would like to get the genomic coordinates of each of them in their corresponding genomes. Is there a way to do this via NCBI E-utils or by any other means?
I'm currently trying out the following approach:
curl -s "http://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=protein&id=108769513&rettype=gp" | grep coded_by=
giving:
/coded_by="CP000384.1:2273218..2274069"
which I'm processing to use 2273218 and 2274069 as the genomic coordinates.
Let me know if this approach is right, and/or if there's a better way to do this.
Thank you.