reverseSeq            package:Biostrings            R Documentation

_R_e_v_e_r_s_e _S_e_q_u_e_n_c_e

_D_e_s_c_r_i_p_t_i_o_n:

     WARNING: The functions described in this man page have been
     deprecated in favor of 'reverse,XString-method' and
     'reverseComplement'.

     Functions to obtain the reverse and reverse complement of a
     sequence

_U_s_a_g_e:

     reverseSeq(seq)
     revcompDNA(seq)
     revcompRNA(seq)

_A_r_g_u_m_e_n_t_s:

     seq: Character vector.  For 'revcompRNA' and 'revcompDNA' the
          sequence should consist of appropriate letter codes:
          '[ACGUN]' and 'ACGTN', respectively.

_D_e_t_a_i_l_s:

     The function reverses the order of the constituent character
     strings of its argument.

_V_a_l_u_e:

     A character vector of the same length as 'seq'.

_A_u_t_h_o_r(_s):

     R. Gentleman, W. Huber, S. Falcon

_S_e_e _A_l_s_o:

     'alphabetFrequency', 'reverseComplement'

_E_x_a_m_p_l_e_s:

      w <-  c("hey there", "you silly fool")
      if (interactive()) {
        reverseSeq(w)  # deprecated (inefficient on large vectors)
      }
      reverse(BStringSet(w))  # more efficient

      w <- "able was I ere I saw Elba"
      if (interactive()) {
        reverseSeq(w)  # deprecated (inefficient on large vectors)
      }
      reverse(BStringSet(w))  # more efficient

      rna1 <- "UGCA"
      if (interactive()) {
        revcompRNA(rna1)  # deprecated (inefficient on large vectors)
      }
      reverseComplement(RNAString(rna1))  # more efficient

      dna1 <- "TGCA"
      if (interactive()) {
        revcompDNA(dna1)  # deprecated (inefficient on large vectors)
      }
      reverseComplement(DNAString(dna1))  # more efficient

      ## Comparing efficiencies:
      if (interactive()) {
        library(hgu95av2probe)
        system.time(y1 <- reverseSeq(hgu95av2probe$sequence))
        x <- DNAStringSet(hgu95av2probe$sequence)
        system.time(y2 <- reverse(x))

        system.time(y3 <- revcompDNA(hgu95av2probe$sequence))
        system.time(y4 <- reverseComplement(x))
      }

