# DNA 序列檔名更改與選擇
###### tags: `R教學` `序列整理相關`

練習檔案：[2920circular.fas](https://drive.google.com/file/d/18aYRZgqmqbhAeIkEkV8Pkx_NHBWQzuPM/view?usp=sharing)

rm(list=ls())
setwd(choose.dir())

#if (!requireNamespace("BiocManager", quietly = TRUE))
#  install.packages("BiocManager")
#BiocManager::install(version = "3.12")
#BiocManager::install("Biostrings")

library("Biostrings")

###序列讀取
sequence<- readDNAStringSet("2920circular.fas")

###檔名清洗
names(sequence) <- paste0("asp_", unlist(lapply(
  strsplit(names(sequence), "; "), 
  function(x){x[4]})))

###序列選擇
out.put.seq <- sequence[-grep("trn", names(sequence), value = F),]

###序列輸出
writeXStringSet (out.put.seq, "out.put.seq.fasta")