R 讀入fastq文件

##--構建函數--## seq_import <- function( file ){ seq <- readLines(file) # 讀入序列，每一個元素存入一行 seq <- seq[seq != ""] # 去除空行 is.anno <- regexpr("^>", seq, perl=T) # 正則匹配（regular expression）註釋行,是註釋行爲1，不然爲-1

>>阅读原文<<