获取chr信息的功能

Question

我正在尝试创建一个函数来获取给定基因组中给定chr的大小。按照我到目前为止所做的尝试。

chrx <- as.character("chr1")
Bchrs <- function(x){
  library(BSgenome.Ggallus.UCSC.galGal4)
  genome <- BSgenome.Ggallus.UCSC.galGal4
  boundaries <- NULL
  boundaries$Chrl <- length(genome$x)
  boundaries$Chrn <- 'x'
  return(boundaries)
}

使用此功能我可以期待：

chrx <- "chr1" 
Bchrs(chrx)

  

$ Chrl [1] 195276750

     

$ Chrn [1]＆＃34; chr1＆＃34;

但是我收到了这个错误：

Bchrs(chrx)

  

x [[name]]出错：没有这样的序列

我会感激任何光线使它有效。