Method: BioDSL::Assemble#entry_overlap

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

#entry_overlap ⇒ BioDSL::Seq

Method to extract and upcase the overlapping part of an assembled pair.

Returns:



150
151
152
153
154
155
156
157
158
159
160
161
162
 
# File 'lib/BioDSL/seq/assemble.rb', line 150

def entry_overlap
  if @entry1.qual && @entry2.qual
    entry_overlap1 = @entry1[@offset1...@offset1 + @overlap]
    entry_overlap2 = @entry2[@offset2...@offset2 + @overlap]

    entry = merge_overlap(entry_overlap1, entry_overlap2)
  else
    entry = @entry1[@offset1...@offset1 + @overlap]
  end

  entry.seq.upcase!
  entry
end