# NOT RUN {
library(phyclust, quiet = TRUE)
X <- seq.data.toy$org
X[, 5] <- .nucleotide$nid[.nucleotide$code == "-"]
ret <- phyclust(X, 2)
Mu.GAPs <- prune.Mu(X, ret$class.id, ret$Mu)
ret$Mu[, 5]
Mu.GAPs[, 5] # Replace by GAPs.
# }
Run the code above in your browser using DataLab