inst/doc/Sequences-Vignette.R

In dowser: B Cell Receptor Phylogenetics Toolkit

## ----eval=TRUE, warning=FALSE, message=FALSE----------------------------------
library(dowser)

data(ExampleClones)

# Collapse nodes with identical sequences. This will 
trees = collapseNodes(ExampleClones[1:2,])

# Plot trees with node ID numbers
plots = plotTrees(trees, tips="c_call", tipsize=2, node_nums=TRUE, labelsize=7)

plots[[1]]

sequence = getNodeSeq(trees, node=50, clone=3128)

print(sequence)