This function shows the population integer slot from a
DNAStringSet
or an AAStringSet
metadata
information.
popinteger(seq)
DNAStringSet
or AAStringSet
[mandatory]
population integer from metadata
## load example sequence data
data(iupac, package="MSA2dist")
iupac.aa <- iupac |> cds2aa(shorten = TRUE)
## create poplist
poplist <- list(FRA = grep("Mmd.FRA", names(iupac)),
GER = grep("Mmd.GER", names(iupac)),
IRA = grep("Mmd.IRA", names(iupac)),
AFG = grep("Mmm.AFG", names(iupac)))
iupac.aa <- iupac.aa |> addpop2string(poplist)
popinteger(iupac.aa)
#> $FRA
#> [1] 1 2 3 4 5 6 7 8
#>
#> $GER
#> [1] 9 10 11 12 13 14 15 16
#>
#> $IRA
#> [1] 17 18 19 20 21 22 23 24
#>
#> $AFG
#> [1] 25 26 27 28 29 30
#>