added qual clipping to ruby code

git-svn-id: http://biopieces.googlecode.com/svn/trunk@1276 74ccb610-7750-0410-82ae-013aeee3265d

diff --git a/code_ruby/Maasha/lib/seq.rb b/code_ruby/Maasha/lib/seq.rb
index 72ca1cd62c742b65bcb0787bb3ddb6cba3cde83f..de13ec21b034fe8a65074d8aaff8afbc17485121 100644 (file)
--- a/code_ruby/Maasha/lib/seq.rb
+++ b/code_ruby/Maasha/lib/seq.rb
@@ -336,7 +336,8 @@ class Seq
     pos = self.adaptor_locate_right(adaptor)
 
     if pos > 0
-      self.seq = self.seq[0 ... pos]
+      self.seq  = self.seq[0 ... pos]
+      self.qual = self.qual[0 ... pos] unless self.qual.nil?
     end
   end
 
@@ -347,7 +348,8 @@ class Seq
     pos = self.adaptor_locate_left(adaptor)
 
     if pos > 0
-      self.seq = self.seq[pos + 1 ... self.seq.length]
+      self.seq  = self.seq[pos + 1 ... self.seq.length]
+      self.qual = self.qual[pos + 1 ... self.qual.length] unless self.qual.nil?
     end
   end
 

diff --git a/code_ruby/Maasha/test/test_seq.rb b/code_ruby/Maasha/test/test_seq.rb
index 86687386e3403e5f03772a9f01d48f5e399f3356..2482422611ca53fbdf9a5cb5e61b3848011706e7 100755 (executable)
--- a/code_ruby/Maasha/test/test_seq.rb
+++ b/code_ruby/Maasha/test/test_seq.rb
@@ -290,12 +290,26 @@ class TestSeq < Test::Unit::TestCase
     assert_equal( "nnnnn", @entry.seq)
   end
 
+  def test_Seq_adaptor_clip_right_returns_correct_qual
+    @entry.seq  = "nnnnncgat"
+    @entry.qual = "abcdefghi"
+    @entry.adaptor_clip_right("cgat")
+    assert_equal( "abcde", @entry.qual)
+  end
+
   def test_Seq_adaptor_clip_left_returns_correct_sequence
     @entry.seq = "cgatnnnnn"
     @entry.adaptor_clip_left("cgat")
     assert_equal( "nnnnn", @entry.seq)
   end
 
+  def test_Seq_adaptor_clip_left_returns_correct_qual
+    @entry.seq  = "cgatnnnnn"
+    @entry.qual = "abcdefghi"
+    @entry.adaptor_clip_left("cgat")
+    assert_equal( "efghi", @entry.qual)
+  end
+
   def test_Digest_new_raises_on_bad_pattern_residue
     assert_raise(DigestError) { Digest.new(@entry, "X", 4) }
   end