Question: Counting a number of mobile elements in a list of genomes using BioPython
gravatar for little_more
12 months ago by
little_more20 wrote:

Suppose I have a big list of assembly IDs:

['GCA_000005845.2',  'GCA_000006925.2',  'GCA_000007405.1',  'GCA_000007445.1', ...]

And I need to count a number of mobile elements in each chromosome in all assemblies. How could I do this using BioPython? It seems that I need to download all genbank files and count mobile elements by parsing each one of them. Unfortunately, I can't find a function to download a full genbank file given an assembly ID.

I tried to do this since the direct search for "GCA_000005845.2" in assembly db doesn't work: = "" # Always tell NCBI who you are
handle = Entrez.esearch(db="assembly", term="GCA_000005845.2")
record =
handle1=Entrez.esummary(db="assembly",id=record['IdList'], report="full")

But I don't know what to do next as the result does not contain GenBank ID (U00096.3). So is there a way to do this with BioPython or maybe someone knows another simple way?

biopython genome • 292 views
ADD COMMENTlink written 12 months ago by little_more20
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1593 users visited in the last hour