Contents

require 'fileutils'

module FastaRead
  class Sequence
    include Methadone::Main
    include Methadone::CLILogging

    attr_reader :separate_lines

    def initialize(fasta, chromosome, cstart, cend)
      @fasta =  fasta.gsub("\n", "")
      @chromosome = ">chr#{chromosome}"
      @cstart = cstart.to_i
      @cend =   cend.to_i
    end

    def process
      @fasta[real_cstart..real_cend]
    end

    def real_cstart
      cstart = @cstart
      cstart -= 1 unless @cstart == 0
      cstart + @chromosome.length
    end

    def real_cend
      @cend - 1 + @chromosome.length
    end
  end
end

Version data entries

1 entries across 1 versions & 1 rubygems

Version	Path
fasta_read-1.1.0	lib/fasta_read/sequence.rb