У меня есть файл fasta с несколькими последовательностями ДНК, которые я хочу преобразовать в объект DNAStringSet. Вот мои команды:
myfasta<-read.fasta('myfasta.fasta', as.string = TRUE)
set = NULL
for (i in 1:length(myfasta)){
set = c(set, myfasta[[i]][[1]], collapse="")
set}
myset <- DNAStringSet(set)
Это прекрасно работает для меня, проблема в том, чтобы сохранить имя каждой последовательности, я пытался включить имя в аргумент «collapse», но затем функцию DNAStringSet выдает ошибку.
Знаете ли вы, как я могу сохранить имена последовательностей фаста в объекте DNAStringSet?
Спасибо!