# NOT RUN {
# create a sequence object
s<-Sequence(alphabets=list(NucleotideAlphabet()),string="AATTGCCCCCTTGG")
# get all Site states
getStates(s)
# get the states for a collection of sites
getStates(s,c(1:3,5,8))
# get states via virtual field
s$states
# }
Run the code above in your browser using DataLab