]> git.donarmstrong.com Git - biopieces.git/blobdiff - code_ruby/lib/maasha/seq.rb
fixed indels in methods mask_seq_hard/soft
[biopieces.git] / code_ruby / lib / maasha / seq.rb
index d22bfd1c9a1628fe83c72f2566f5663a11ab8d68..e4b0a947bd5460c5b0a12aa69c63bc84c4b43775 100644 (file)
@@ -457,6 +457,7 @@ class Seq
     na_seq  = NArray.to_na(self.seq, "byte")
     na_qual = NArray.to_na(self.qual, "byte")
     mask    = (na_qual - SCORE_ILLUMINA) < cutoff
+    mask   *= na_seq.ne("-".ord)
 
     na_seq[mask] = 'N'.ord
 
@@ -475,6 +476,7 @@ class Seq
     na_seq  = NArray.to_na(self.seq, "byte")
     na_qual = NArray.to_na(self.qual, "byte")
     mask    = (na_qual - SCORE_ILLUMINA) < cutoff
+    mask   *= na_seq.ne("-".ord)
 
     na_seq[mask] ^= ' '.ord