# NOT RUN {
# create a nucleotide sequence with a process attached
s<-NucleotideSequence(string="ATGGCCA",processes=list(list(JC69())))
# plot total rates in range 1:4
plot(s,1:4)
# plot all total rates
plot(s)
# }
Run the code above in your browser using DataLab