Disambiguate {DECIPHER} | R Documentation |
Performs the inverse function of ConsensusSequence
by expanding any ambiguities present in sequences.
Disambiguate(myXStringSet)
myXStringSet |
A |
Ambiguity codes in the IUPAC_CODE_MAP
can be used to represent multiple nucleotides at a single position. Using these letters, multiple oligonucleotide permutations can be represented with a single ambiguous sequence. This function expands each sequence in the DNAStringSet
input into all of its permutations. Note that sequences with many ambiguities can result in a very large number of potential permutations.
A DNAStringSetList
or RNAStringSetList
with one element for each sequence in myXStringSet
.
Erik Wright eswright@pitt.edu
dna <- DNAStringSet(c("ACST", "NNN")) dna_list <- Disambiguate(dna) dna_list[[1]] dna_list[[2]] unlist(dna_list) rna <- RNAStringSet(c("ACGU", "AGAU")) # 2 permutations rna <- ConsensusSequence(rna) # "ASRU" Disambiguate(rna) # 4 permutations