尝试通过 getPromoterSeq(GenomicFeatures 包)从 TxDb.Hsapiens.UCSC.hg19.knownGene 获取启动子序列时失败
library(TxDb.Hsapiens.UCSC.hg19.knownGene)
library(BSgenome.Hsapiens.UCSC.hg19)
library(Homo.sapiens)
promoter.gr = transcriptsBy(TxDb.Hsapiens.UCSC.hg19.knownGene, by='gene')[11338]
promoter.seq = getPromoteSeq(promoter.gr, Hsapiens, upstream=3000, downstream=0)
promoter.seq
结果如下:
DNAStringSetList of length 1
Error in .Call2("new_CHARACTER_from_XStringSet", x, xs_dec_lkup(x), PACKAGE='Biostrings'):
key 112 (char 'p') not in lookup table.
我搜索有关此错误的其他讨论。当人们试图加载序列文件时,似乎有非 T、C、G、A 序列,主要是序列文件中的“p”。虽然在这里我什至没有在本地提供序列,但它全部来自包。