This function is a fork from an internal function from
Biostrings
makePostalignedSeqs(x)
x
get internal function makePostalignedSeqs
## define two cds sequences
cds1 <- Biostrings::DNAString("ATGCAACATTGC")
cds2 <- Biostrings::DNAString("ATGCATTGC")
makePostalignedSeqs(Biostrings::pairwiseAlignment(
cds2aa(Biostrings::DNAStringSet(cds1)),
cds2aa(Biostrings::DNAStringSet(cds2))))
#> [[1]]
#> AAStringSet object of length 2:
#> width seq names
#> [1] 4 MQHC
#> [2] 4 M-HC
#>
#> [[2]]
#> IRanges object with 2 ranges and 0 metadata columns:
#> start end width
#> <integer> <integer> <integer>
#> [1] 1 4 4
#> [2] 1 3 3
#>
#> [[3]]
#> 4-letter BString object
#> seq: | ||
#>