@q_entry.seq.downcase!
@s_entry.seq.downcase!
- space = Space.new( 0, 0, @q_entry.length - 1, @s_entry.length - 1)
- kmer = 64
+ space = Space.new(0, 0, @q_entry.length - 1, @s_entry.length - 1)
+ kmer = 32
align_recurse(@q_entry.seq, @s_entry.seq, space, kmer)
matches_upcase
matches = matches_select_by_score(matches, space)
end
-# $stderr.puts "#{kmer} #{matches.size} #{space.q_dim} x #{space.s_dim}"
kmer /= 2
end
score_diag = match_score_diag(match, space)
match.score = score_length + score_diag
-
-
-# $stderr.puts "score_length: #{score_length} score_diag: #{score_diag} score: #{match.score}" # DEBUG
-# $stderr.puts match
end
matches.sort_by! { |match| match.score }