Complement a DNA sequence

Question

Suppose I have a DNA sequence. I want to get the complement of it. I used the following code but I am not getting it. What am I doing wrong ?

s=readline()
AT         


        

Answer 1

The Bioconductor package Biostrings has many useful functions for this sort of operation. Install once:

source("http://bioconductor.org/biocLite.R")
biocLite("Biostrings")

then use

library(Biostrings)
dna = DNAStringSet(c("ATCTCGGCGCGCATCGCGTACGCTACTAGC", "ACCGCTA"))
complement(dna)