From 1a49224a5a43505da023897a29fb24df81a18888 Mon Sep 17 00:00:00 2001 From: martinahansen Date: Tue, 1 Nov 2011 13:43:51 +0000 Subject: [PATCH] fixed unit tests for reverse and revcomp in seq.rb git-svn-id: http://biopieces.googlecode.com/svn/trunk@1607 74ccb610-7750-0410-82ae-013aeee3265d --- code_ruby/lib/maasha/seq.rb | 19 +++++++++++++++++++ code_ruby/test/maasha/test_seq.rb | 6 +++--- 2 files changed, 22 insertions(+), 3 deletions(-) diff --git a/code_ruby/lib/maasha/seq.rb b/code_ruby/lib/maasha/seq.rb index 037f99d..c1ce2da 100644 --- a/code_ruby/lib/maasha/seq.rb +++ b/code_ruby/lib/maasha/seq.rb @@ -239,6 +239,7 @@ class Seq def reverse self.seq.reverse! self.qual.reverse! if self.qual + self end # Method that complements sequence including ambiguity codes. @@ -327,6 +328,24 @@ class Seq self.subseq(start, length) end + def quality_trim(min) + end + + def quality_trim_right(min) + end + + def quality_trim_left(min) + end + + def quality_trim!(min) + end + + def quality_trim_right!(min) + end + + def quality_trim_left!(min) + end + # Method that returns the residue compositions of a sequence in # a hash where the key is the residue and the value is the residue # count. diff --git a/code_ruby/test/maasha/test_seq.rb b/code_ruby/test/maasha/test_seq.rb index 78a6543..29990b7 100755 --- a/code_ruby/test/maasha/test_seq.rb +++ b/code_ruby/test/maasha/test_seq.rb @@ -196,7 +196,7 @@ class TestSeq < Test::Unit::TestCase def test_Seq_reverse_returns_correctly @entry.seq = "ATCG" - assert_equal("GCTA", @entry.reverse) + assert_equal("GCTA", @entry.reverse.seq) end def test_Seq_complement_raises_if_no_sequence @@ -225,13 +225,13 @@ class TestSeq < Test::Unit::TestCase def test_Seq_reverse_complement_for_DNA_is_correct @entry.seq = 'ATCGatcg' @entry.type = 'dna' - assert_equal("cgatCGAT", @entry.reverse_complement) + assert_equal("cgatCGAT", @entry.reverse_complement.seq) end def test_Seq_reverse_complement_for_RNA_is_correct @entry.seq = 'AUCGaucg' @entry.type = 'rna' - assert_equal("cgauCGAU", @entry.reverse_complement) + assert_equal("cgauCGAU", @entry.reverse_complement.seq) end def test_Seq_hamming_distance_returns_correctly -- 2.39.2