Method: BioDSL::Trim#quality_trim

Defined in:
lib/BioDSL/seq/trim.rb

#quality_trim(min_qual, min_len = 1) ⇒ Object

Method to progressively trim a Seq object sequence from both ends until a run of min_len residues with quality scores above min_qual is encountered.



79
80
81
82
83
84
85
86
87
88
89
90
 
# File 'lib/BioDSL/seq/trim.rb', line 79

def quality_trim(min_qual, min_len = 1)
  check_trim_args(min_qual, min_len)

  pos_right = trim_right_pos_c(@qual, length, min_qual, min_len,
                               Seq::SCORE_BASE)
  pos_left  = trim_left_pos_c(@qual, length, min_qual, min_len,
                              Seq::SCORE_BASE)

  pos_left = pos_right if pos_left > pos_right

  self[pos_left...pos_right]
end