Annotatr: include extra columns
library(annotatr) read_annotations(con = "test.bed", genome = 'xxx', name = 'test', format = 'bed', extraCols=c(gene_id="character", symbol="character", tx_id="character")) print(annotatr_cache$get('xxx_custom_test')) test.bed chr1 199 240 gene1 abc mm10 chr2 200 400 gene2 def mm10