]> git.donarmstrong.com Git - biopieces.git/commitdiff
added reverse method to ruby Seq class
authormartinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Thu, 16 Sep 2010 10:08:39 +0000 (10:08 +0000)
committermartinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Thu, 16 Sep 2010 10:08:39 +0000 (10:08 +0000)
git-svn-id: http://biopieces.googlecode.com/svn/trunk@1090 74ccb610-7750-0410-82ae-013aeee3265d

code_ruby/Maasha/lib/seq.rb

index 2dd9c5fef3b3c378c47422d58efffa9e92aac82c..e00590c95254b4b3388e8be890ab66ed88f9eae4 100644 (file)
@@ -55,9 +55,9 @@ class Seq
     raise SeqError, "Missing seq_name" if self.seq_name.nil?
     raise SeqError, "Missing seq"      if self.seq.nil?
     record             = {}
-    record['SEQ_NAME'] = self.seq_name
-    record['SEQ']      = self.seq
-    record['SEQ_LEN']  = self.length
+    record[:SEQ_NAME] = self.seq_name
+    record[:SEQ]      = self.seq
+    record[:SEQ_LEN]  = self.length
     record
   end
 
@@ -80,6 +80,19 @@ class Seq
     ">#{seq_name}\n#{seq}\n"
   end
 
+  # Method to reverse complement sequence.
+  def reverse_complement
+    self.reverse
+    self.complement
+  end
+
+  alias revcomp reverse_complement
+
+  # Method to reverse the sequence.
+  def reverse
+    self.seq.reverse!
+  end
+
   # Method that complements sequence including ambiguity codes.
   def complement
     raise SeqError, "Cannot complement 0 length sequence" if self.length == 0