@entry = Seq.new
end
+ test "Seq.new with differnet length SEQ and SCORES raises" do
+ assert_raise(SeqError) { Seq.new("test", "ATCG", "dna", "hhh") }
+ end
+
test "Seq.new_bp returns correctly" do
record = {:SEQ_NAME => "test", :SEQ => "ATCG", :SEQ_TYPE => :dna, :SCORES => "hhhh"}
seq = Seq.new_bp(record)
assert_equal(1, seq1.hamming_distance(seq2))
end
+ test "#hamming_distance with ambiguity codes return correctly" do
+ seq1 = Seq.new("test1", "ATCG")
+ seq2 = Seq.new("test2", "atng")
+
+ assert_equal(1, seq1.hamming_distance(seq2))
+ assert_equal(0, seq1.hamming_distance(seq2, ambiguity: true))
+ end
+
test "#edit_distance returns correctly" do
seq1 = Seq.new("test1", "ATCG")
seq2 = Seq.new("test2", "tgncg")