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!