This function shows the position slot from a DNAStringSet or an AAStringSet metadata information.

getpos(seq)

Arguments

seq

DNAStringSet or AAStringSet [mandatory]

Value

GenomicRanges information from metadata

See also

Author

Kristian K Ullrich

Examples

## 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