]> git.donarmstrong.com Git - biopieces.git/commitdiff
fixed indels in methods mask_seq_hard/soft
authormartinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Tue, 24 Jan 2012 08:34:23 +0000 (08:34 +0000)
committermartinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Tue, 24 Jan 2012 08:34:23 +0000 (08:34 +0000)
git-svn-id: http://biopieces.googlecode.com/svn/trunk@1734 74ccb610-7750-0410-82ae-013aeee3265d

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