variant_bc: Calculates non-reference base calls

In dieterich-lab/JACUSA2helper: Post-processing for JACUSA2 output

Calculates non-reference base calls

Description

Calculates non-reference base call counts from a vector bc and ref(ference) bases. Make sure the combination of base calls and reference bases <= 2 alleles.

Usage

variant_bc(ref, bc)

Arguments

ref	vector of reference bases..
bc	vector of base calls.

Value

vector of most frequent base call counts.

Examples

ref <- c("A", "A")
bc <- c("A", "AG")
# show non-reference base calls
variant_bc(ref, bc)

dieterich-lab/JACUSA2helper documentation built on March 1, 2023, 12:09 a.m.