This function returns a DNAStringSet
reduced by all
sites containing any gaps ("-", "+", ".") or missing ("N") sites.
globalDeletion(dna)
DNAStringSet
[mandatory]
DNAStringSet
## define two cds sequences
cds1 <- Biostrings::DNAString("ATGCAACATTGC")
cds2 <- Biostrings::DNAString("ATG---CATTGC")
cds1.cds2.aln <- c(Biostrings::DNAStringSet(cds1),
Biostrings::DNAStringSet(cds2))
globalDeletion(cds1.cds2.aln)
#> DNAStringSet object of length 2:
#> width seq
#> [1] 9 ATGCATTGC
#> [2] 9 ATGCATTGC