Retrieve all annotated features in a GenBank record
0
0
Entering edit mode
10 months ago
al-ash ▴ 150

Is edirect supporting retrieval of all annotated features from a record?

Example: I can retrieve all features annotated as "gene" in GenBank record of a mitochondrial genome sequence :

efetch -db nucleotide -id KP026275  -format gene_fasta

However, the retrieved multifasta is missing the features not annotated as "gene", in this case two ribosomal RNAs. I would like to retrieve all features.

I can think of one solution, where I would first retrieve the feature table, then parse the genomic coordinates of all features and use them to extract the sequences of all the features, but this feels like re-inventing something which might be already implemented in edirect. Thanks for any hints & suggestions!

efetch edirect genbank feature annotations • 291 views
ADD COMMENT
1
Entering edit mode

Don't think there is a way to get those directly without parsing the feature table.

ADD REPLY

Login before adding your answer.

Traffic: 1343 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6