This function converts a AAStringSet
into an ape
DNAbin
.
aastring2aabin(aa)
AAStringSet
[mandatory]
An object of class DNAbin
## define two cds sequences
cds1 <- Biostrings::DNAString("ATGCAACATTGC")
cds2 <- Biostrings::DNAString("ATG---CATTGC")
cds1.cds2.aln <- c(Biostrings::DNAStringSet(cds1),
Biostrings::DNAStringSet(cds2))
## convert into AAbin
#aastring2aabin(cds2aa(cds1.cds2.aln))
cds1.cds2.aln |> cds2aa() |> aastring2aabin()
#> 2 amino acid sequences in a list
#>
#> All sequences of the same length: 4
#>