R 读入fastq文件

##--构建函数--## seq_import <- function( file ){ seq <- readLines(file) # 读入序列，每一个元素存入一行 seq <- seq[seq != ""] # 去除空行 is.anno <- regexpr("^>", seq, perl=T) # 正则匹配（regular expression）注释行,是注释行为1，不然为-1

>>阅读原文<<