seqdist<-function(x,n) { dmat<-matrix(nrow=n,ncol=n) for(i in 1:n){ for (j in 1:n){ dmat[i,j]<-length(x[j,][x[j,]!=x[i,]])/length(x[1,])}} return(dmat) }