Class: Bio::MAF::FASTARangeReader

Inherits:
Object
  • Object
show all
Defined in:
lib/bio/maf/tiler.rb

Constant Summary collapse

GT =
'>'.getbyte(0)

Instance Attribute Summary collapse

Instance Method Summary collapse

Constructor Details

#initialize(fspec) ⇒ FASTARangeReader



190
191
192
193
194
195
196
197
198
199
200
201
202
# File 'lib/bio/maf/tiler.rb', line 190

def initialize(fspec)
  if fspec.respond_to? :seek
    @f = fspec
  else
    reader_class = if fspec =~ /.gz$/
                     Zlib::GzipReader
                   else
                     File
                   end
    @f = reader_class.open(fspec)
  end
  position_at_start
end

Instance Attribute Details

#fObject (readonly)

Returns the value of attribute f



188
189
190
# File 'lib/bio/maf/tiler.rb', line 188

def f
  @f
end

#posObject (readonly)

Returns the value of attribute pos



188
189
190
# File 'lib/bio/maf/tiler.rb', line 188

def pos
  @pos
end

Instance Method Details

#position_at_startObject



206
207
208
209
210
# File 'lib/bio/maf/tiler.rb', line 206

def position_at_start
  first = f.readline
  raise "expected FASTA comment" unless first =~ /^>/
  @pos = 0
end

#read_interval(z_start, z_end) ⇒ Object



212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
# File 'lib/bio/maf/tiler.rb', line 212

def read_interval(z_start, z_end)
  if z_start < pos
    position_at_start
  end
  data = ''
  region_size = z_end - z_start
  in_region = false
  f.each_line do |line_raw|
    if line_raw.getbyte(0) == GT
      raise "unexpected description line: #{line_raw.inspect}"
    end
    line = line_raw.strip
    end_pos = pos + line.size
    if (! in_region) && pos <= z_start && z_start < end_pos
      offset = z_start - pos
      end_offset = [(offset + region_size), line.size].min
      data << line.slice(offset...end_offset)
      in_region = true
    elsif in_region
      need = region_size - data.size
      raise "should not happen: region #{region_size}, data #{data.size}, need #{need}" if need < 0
      if need > line.size
        data << line
      else
        # last line
        data << line.slice(0, need)
        break
      end
    end
    @pos = end_pos
  end
  return data
end