本工作目录下的sample_dna.txt文件的DNA序列弄互补链

rm(list = ls())
library(Biostrings)
data_dna = readLines("sample_data.txt")
data_dna
DNAString(data_dna)
reverseComplement(DNAString(data_dna))
toString(reverseComplement(DNAString(data_dna)))
data_dna_a =toString(reverseComplement(DNAString(data_dna)))
data_dna_a
write(data_dna_a,file = "data_sample.txt")