This function shows the region used slot from a
DNAStringSet
or an AAStringSet
metadata
information.
regionused(seq)
DNAStringSet
or AAStringSet
[mandatory]
population names from metadata
## load example sequence data
data("hiv", package="MSA2dist")
## create mask
mask1 <- IRanges::IRanges(start=c(11,41,71), end=c(20,50,80))
## use mask
hiv.region <- hiv |> cds2aa() |> string2region(mask=mask1)
#(hiv.region |> slot("metadata"))$regionUsed
hiv.region |> regionused()
#> IRanges object with 4 ranges and 0 metadata columns:
#> start end width
#> <integer> <integer> <integer>
#> [1] 1 10 10
#> [2] 21 40 20
#> [3] 51 70 20
#> [4] 81 91 11