fasta 파일을 Biostrings 패키지를 사용해서 R의 data frame 형태로 바꾸어 보자
library("Biostrings")
fastaFile= readDNAStringSet("..ASVs_rep.fastq.txt")
seq_name = names(fastaFile)
sequence = paste(fastaFile)
df <- data.frame(seq_name, sequence)
write.csv(df, "../ASV_table.csv", row.names = F)
반대로 dataframe을 fasta로 저장하는 방법은 아래 링크 참고 바람
| 출처
반응형