| 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