This function shows the position slot from a
DNAStringSet
or an AAStringSet
metadata
information.
getpos(seq)
DNAStringSet
or AAStringSet
[mandatory]
GenomicRanges
information from metadata
## load example sequence data
data(iupac, package="MSA2dist")
## add position
iupac <- iupac |> addpos2string(chrom="chr1", start=1, end=1000)
#(iupac |> slot("metadata"))$GRanges
iupac |> getpos()
#> GRanges object with 1 range and 0 metadata columns:
#> seqnames ranges strand
#> <Rle> <IRanges> <Rle>
#> [1] chr1 1-1000 *
#> -------
#> seqinfo: 1 sequence from an unspecified genome; no seqlengths