data(bcrdata)
bcr_data <- bcrdata$data
clonal_family <- bcr_data %>%
seqCluster(seqDist(bcr_data), .fixed_threshold = 3) %>%
repGermline(.threads = 1) %>%
repAlignLineage(.min_lineage_sequences = 2, .align_threads = 2, .nofail = TRUE) %>%
repClonalFamily(.threads = 1, .nofail = TRUE)
# This condition can be omitted; it prevents the example from crashing
# when ClustalW or PHYLIP are not installed
if (!("step_failure_ignored" %in% class(clonal_family))) {
vis(clonal_family[["full_clones"]][["TreeStats"]][[2]])
}
Run the code above in your browser using DataLab