#!/usr/bin/Rscript 
a <- readLines("starts.txt")
c <- unlist(splitstr(a, ""))
d <- data.frame(table(c))
#pdf("freqrel.pdf", width=10, height=5)
barplot(xtabs(Freq ~ c, d)/sum(xtabs(Freq ~ c, d)), main="Relative frequencies of amino acids in proteins")
#dev.off()