已知转录本test.ID，根据biomart包批量获得ensembl数据库中对应序列信息，并以fa格式输出

该方法需要魔法，我已进行了调试，可以成功运行

转录本test.ID样式(ENSMUST00000040608.4)，以xlsx形式储存

输入数据格式

library(Biostrings)
library(biomaRt)
library(openxlsx)
library(seqinr)

rm(list = ls())
genelist <- read.xlsx("待测基因testID汇总.xlsx")
genelist <- genelist$test_id
#移除genelist的小数点
genelist<-unlist(lapply(genelist, FUN = function(x) {return(strsplit(x, split = ".",fixed = T)[[1]][1])}))

mart <- useMart("ensembl", dataset="mmusculus_gene_ensembl") # 小鼠的数据库
Ensemble_to_seq <- function(x) {cat("Getting sequence for gene", x, "\n")  # 添加打印语句
                                seq = biomaRt::getSequence(id = x, 
                                type = "ensembl_transcript_id", 
                                seqType = "cdna",