I want to perform some statistics on the lengths of DNA-Sequences inside NCBI Genbank,
for example, I want to draw a histogram of these lengths.
So I am asking is there a way to retrieve the lengths of all DNA-Sequences inside the Genbank? Can it be done using Python or R, or is this information already available via NCBI Genbank?
You can retrieve sequence start and end for your queries. And by using both, you can get the length of sequence which can easily be plotted by any program.
However, in my personal opinion, it would be not so useful plot or analysis if you are just trying to do for everything present in huge database such as Genbank. Anyway, it your choice. But sure it can be done.
Belal Medhat
Download all sequences from GenBank FTP. Use infoseq or another script to calculate the length of sequences.
Abhijeet Singh Thank you so much for your answer, yes I think using BioPython library I can download or retrieve all of the DNA-Sequences then I can calculate the length from the FASTA file as you said.
Sureshkumar v. Thank you so much for your answer, yes I googled the tool Infoseq, and there is a GUI package for it.
- Badges
- Science method