worked on sam - added 1 to phred score check in seq.rb

[biopieces.git] / code_ruby / lib / maasha / seq.rb

diff --git a/code_ruby/lib/maasha/seq.rb b/code_ruby/lib/maasha/seq.rb
index a3fe78905159ab4edeee58fcd0067ee4aca2988c..3b950aef6c3ae6326b65d4a2217f0eef4d672976 100644 (file)
--- a/code_ruby/lib/maasha/seq.rb
+++ b/code_ruby/lib/maasha/seq.rb
@@ -335,7 +335,7 @@ class Seq
   def convert_phred2illumina!
     self.qual.gsub!(/./) do |score|
       score_phred  = score.ord - SCORE_PHRED
-      raise SeqError, "Bad Phred score: #{score} (#{score_phred})" unless (0 .. 40).include? score_phred
+      raise SeqError, "Bad Phred score: #{score} (#{score_phred})" unless (0 .. 41).include? score_phred
       score_illumina = score_phred + SCORE_ILLUMINA
       score          = score_illumina.chr
     end