Calculates reference base to base call (base substitution).

Calculates and formats base changes for ref and bc. All elements in ref must be all length 1.

Usage

base_sub(bases, ref)

Arguments

bases

vector of base calls or tibble of base call counts.

ref

vector of reference bases.

Value

vector of formatted base call substitutions.

Examples

bases <- c("AG", "G", "C", "G")
ref <- c("A", "A", "C", "A")
base_sub(ref, bases)
#> [1] "AG->A"     "G->A"      "no change" "G->A"