Hi,

Apologies for the naive post but I'm wondering if there is a way of exporting information I have on a gene acquired from ensembldb to a GenBank file?

I have the following code, where I fetch the sequence and annotation of the Actin beta gene:

library(EnsDb.Hsapiens.v86)
library(dplyr)

Hs_edb <- EnsDb.Hsapiens.v86
Hs_dna <- getGenomeTwoBitFile(Hs_edb)

ACTB_db <- Hs_edb %>% 
  ensembldb::filter(filter = GeneNameFilter("ACTB")) %>% 
  ensembldb::filter(filter = ~tx_biotype == "protein_coding")
ACTB_gene <- genes(ACTB_db)
ACTB_seq <- getSeq(Hs_dna, ACTB_gene)

This gives me a GRanges object with the features and a DNAStringSet object with the sequence. What I would like is to export the sequence and features into a .gb or .gbk file. The reason I'd like to do this is that many of my colleagues use different forms of software to view sequence information and are not skilled with R. The only commonly interchangeable format that most pieces of software are able to parse is GenBank.

Thank you.

I would be surprised if there were any functionality to do this. I can see a use case for parsing a .gbk file, but what's the use case for generating one? If you want the β-actin .gbk file, why wouldn't you just go to NCBI and download it, particularly since what you will be providing your colleagues is just a small proportion of what is actually in a .gbk file?