> basecomp<-function(inseq)
#Calculate base composition of an input sequence
#inseq is a vector of integers
{
f<-rep(0,4) #vector to store base counts
for(i in 1:4)
{
f[i]<-length(inseq[inseq==i])/length(inseq)
}
return(f)
}