calculate % of identical residues — calc_ident_mat • bstringr

calculate % of identical residues

calc_ident_mat(bstrobj, gap_chr = "-")

Arguments

bstrobj	bstrobj
gap_chr	gap character

Examples

temp <- bstr(c("ATGC--CG", "ACGC-GGG", "--------", "AT"))
round(calc_ident_mat(temp), 2)
#>           no name 1 no name 2 no name 3 no name 4
#> no name 1        NA        NA        NA        NA
#> no name 2      0.57        NA        NA        NA
#> no name 3      0.00      0.00        NA        NA
#> no name 4      0.25      0.12         0        NA