Skip to contents

Calculates non-reference base calls

Source: R/sub-ratio.R
variant_bc.Rd

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)
#> [1] ""  "G"