setwd("D:/ziliao/zhuanye/R bear")

设置工作目录，输入fasta文件所在文件夹

data <- readLines('data2.fasta')

导入序列，存为矩阵数据data

x=data

传导值给x

name <- data[grepl('>',data)]

提取各条序列名称，保存在name向量中

x[grepl('>',data)]="n"

使用特定的字符来分割个条序列

y=x[2]

for (i in 2:length(x))

{

y=paste(y, x[i], sep = "")

}

将序列转化成整体单一字符

y <- strsplit(y,'n')

分割各条序列

y <- unlist(y)

转化为向量格式

这个时候各条序列就保存在y向量的不同单元格中了