bioroebe 0.10.80 → 0.11.32

data/lib/bioroebe/shell/add.rb

@@ -1,108 +0,0 @@
-#!/usr/bin/ruby -w
-# Encoding: UTF-8
-# frozen_string_literal: true
-# =========================================================================== #
-# require 'bioroebe/shell/add.rb'
-# =========================================================================== #
-module Bioroebe
-
-class Shell < ::Bioroebe::CommandlineApplication
-
-  # ========================================================================= #
-  # === add_to_end
-  #
-  # This method will invoke append() if something has to be appended.
-  # ========================================================================= #
-  def add_to_end(i)
-    add_to_start_or_end(i, :end)
-  end; alias add add_to_end # === add (add tag)
-
-  # ========================================================================= #
-  # === add_to_start
-  # ========================================================================= #
-  def add_to_start(i)
-    add_to_start_or_end(i, :start)
-  end; alias left_add add_to_start # === left_add
-
-  # ========================================================================= #
-  # === add_to_start_or_end
-  #
-  # This method can either add to the start or to the end.
-  #
-  # The default is to append to the nucleotide sequence.
-  #
-  # We can input a number - in this case, we simply add these many
-  # nucleotides onto the main string.
-  # ========================================================================= #
-  def add_to_start_or_end(
-      i                 = '',
-      append_or_prepend = :append
-    )
-    i = i.join('') if i.is_a? Array
-    i = i.dup
-    old_length = i.size
-    case i.to_s # Easier start/stop entries.
-    # ======================================================================= #
-    # === Add a start codon.
-    # ======================================================================= #
-    when 'start',
-         'START',
-         /ORF/i,
-         ':start'
-      i = ::Bioroebe.start_codon? # Used to be hardcoded -> 'ATG'
-    # ======================================================================= #
-    # === Add a stop codon.
-    # ======================================================================= #
-    when /^stop$/i,
-         ':stop'
-      i = ::Bioroebe.stop_codons?.sample
-    end
-    i.upcase!
-    case append_or_prepend
-    # ======================================================================= #
-    # === :prepend
-    # ======================================================================= #
-    when :prepend,
-         :to_start,
-         :start
-      if i =~ /^\d+$/ # If input is only numbers.
-        erev 'Only numbers were given: Prepending '+
-              sfancy(i.to_s)+rev+
-             ' random nucleotides to the main string now.'
-        i.to_i.times { prepend(add_nucleotide) }
-      else
-        erev "Prepending #{sfancy(i)}#{rev} to the main string."
-        prepend(i)
-      end
-    # ======================================================================= #
-    # === :append
-    # ======================================================================= #
-    when :append,
-         :to_end,
-         :end
-      if i =~ /^\d+$/ # If input is only numbers.
-        erev "Only numbers were given: Adding #{sfancy(i.to_s)}#{rev}"\
-             " random nucleotides to the main string now."
-        i.to_i.times { append(return_random_nucleotide) }
-      else
-        if only_nucleotides?(i)
-          msg = "Adding #{sfancy(i)}#{rev}"
-          msg = msg.dup if msg.frozen?
-          if ::Bioroebe.is_a_stop_codon? i
-            msg << ' (a stop codon)'
-          end
-          msg << ' to the main string.'
-          erev msg
-        end
-        append(i)
-      end
-    end
-    new_length = string?.size.to_s
-    unless new_length.to_i == old_length.to_i
-      erev "The new length of the main string is now: "\
-           "#{simp(new_length)}#{rev}."
-    end
-    show_dna_sequence
-  end
-
-end; end

data/lib/bioroebe/shell/assign.rb

@@ -1,360 +0,0 @@
-#!/usr/bin/ruby -w
-# Encoding: UTF-8
-# frozen_string_literal: true
-# =========================================================================== #
-# require 'bioroebe/shell/assign.rb'
-# =========================================================================== #
-module Bioroebe
-
-class Shell < ::Bioroebe::CommandlineApplication
-
-  # ========================================================================= #
-  # === set_dna_sequence                          (assign tag, assigning tag)
-  #
-  # This method can be used to set/assign the main DNA string.
-  #
-  # It is a fairly long method, mostly because it will do lots of
-  # additional tasks, way aside from setting/assigning to a DNA
-  # sequence only.
-  #
-  # As of June 2016, the method will also keep a backup of the
-  # generated sequence in a local file as well. This will allow
-  # us to "replay" the given sequence on startup of the shell.
-  #
-  # Note that as of Jun 2016, we will chop off any '"' found in
-  # the input String.
-  # ========================================================================= #
-  def set_dna_sequence(
-      i          = nil,
-      be_verbose = true,
-      do_upcase  = :check_for_config_value_here
-    )
-    if is_the_main_sequence_frozen?
-      report_that_the_main_sequence_is_frozen
-      return
-    end
-    case do_upcase
-    # ======================================================================= #
-    # === :check_for_config_value_here
-    # ======================================================================= #
-    when :check_for_config_value_here,
-         :default
-      if @config and @config.respond_to?(:upcase_nucleotides)
-        do_upcase = @config.upcase_nucleotides
-      else
-        do_upcase = false
-      end
-    # ======================================================================= #
-    # === :do_not_upcase
-    # ======================================================================= #
-    when :do_not_upcase
-      do_upcase = false
-    end
-    if do_upcase
-      shall_we_upcase = do_upcase if @config
-    end
-    # ======================================================================= #
-    # === Sanitize the second variable next
-    # ======================================================================= #
-    case be_verbose
-    # ======================================================================= #
-    # === :be_verbose
-    # ======================================================================= #
-    when :be_verbose
-      be_verbose = true
-    # ======================================================================= #
-    # === :do_not_upcase
-    # ======================================================================= #
-    when :do_not_upcase,
-         :no_upcase,
-         :no_upcasing,
-      shall_we_upcase = false
-      be_verbose = true # Sync this.
-    # ======================================================================= #
-    # === :be_silent
-    # ======================================================================= #
-    when :be_silent,
-         :be_quiet
-      be_verbose = false
-    end
-    if i.nil?
-      # ===================================================================== #
-      # === Hande nil value for the variable i next:
-      #
-      # Here in this entry point, we will handle if the first input
-      # argument is nil.
-      #
-      # If @internal_hash[:misc_sequence] is NOT nil, then it will be used
-      # once, then cleared. Otherwise, we will use dna_sequence?.
-      # ===================================================================== #
-      if @internal_hash[:misc_sequence]
-        i = @internal_hash[:misc_sequence]
-        @internal_hash[:misc_sequence] = nil # Set it to nil again.
-      else
-        i = dna_sequence?
-      end
-    end
-    # ======================================================================= #
-    # === Convert Array into String next
-    # ======================================================================= #
-    i = i.join if i.is_a? Array
-    i = i.to_s # Work on a String past this point here.
-    if i.include? ' ' # Sanitize the input a little bit.
-      i = i.strip.delete(' ')
-    end
-    if i.include? '|' # I think we will not need '|' in our main string.
-      i.delete!('|')
-    end
-    # ======================================================================= #
-    # === If the input has only numbers
-    # ======================================================================= #
-    if i =~ /^\d+$/
-      i = random_dna_sequence(i)
-    end
-    # ======================================================================= #
-    # === Handle Fasta format next
-    # ======================================================================= #
-    if i.start_with?('>') and i.include?(N)
-      erev 'It seems as if you have a FASTA sequence there, as it starts with '\
-           'a > character.'
-      erev 'We will omit this first section (the header) though.'
-      i[0 .. i.index(N)] = '' # This does the chop-off action.
-    end
-    i.delete!('-') if i.include? '-'
-    i.delete!('_') if i.include? '_' # Don't want '_' characters.
-    i.delete!('?') if i.include? '?'
-    i.delete!('#') if i.include? '#' # This was added at 08.05.2016.
-    i.delete!('"') if i.include? '"' # This was added at 02.06.2016.
-    i.gsub!(/\^C/,'') if i.include? '^C'
-    i = i.to_s.dup
-    case i
-    # ======================================================================= #
-    # === :GFP
-    # ======================================================================= #
-    when ':GFP' # ← This means the default GFP sequence.
-      i = return_default_GFP_sequence
-    # ======================================================================= #
-    # === random
-    # ======================================================================= #
-    when /^random$/i,
-         '' # As of 04.01.2015, we will try '' too here.
-      i = random :do_not_show_the_string
-    # ======================================================================= #
-    # === multiline
-    # ======================================================================= #
-    when 'multiline' # Handle multiline input.
-      erev 'Multiline input detected.'
-      erev 'Finish by issuing __ on an empty line. (These are 2 "_" tokens.)'
-      i  = $stdin.gets('__')
-      # ===================================================================== #
-      # Next, sanitize it a bit if it starts with a '>' identifier.
-      # ===================================================================== #
-      if i.strip.start_with? '>'
-        i = i[(i.index(N)+1)..-1]
-      end
-    else
-      if i
-        # =================================================================== #
-        # === Next upcase the input if the flag was set
-        # =================================================================== #
-        i.upcase! if shall_we_upcase
-        case mode?
-        # =================================================================== #
-        # === :dna
-        # =================================================================== #
-        when :dna,
-             :rna
-          dna_sequence_object?.set_mode(mode?)
-        # =================================================================== #
-        # === :aminoacids
-        # =================================================================== #
-        when :aminoacids
-          set_aminoacids(i)
-        end
-      else # else input is nil.
-        erev 'Missing Input to assign_sequence(). Please input '\
-             'your DNA string now (Hitting the enter key will '\
-             'finish this):'
-        i = $stdin.gets.chomp
-        i = i.upcase if shall_we_upcase
-      end
-    end
-    i = '' unless i
-    # ======================================================================= #
-    # === Remove possible newlines in the given input
-    # ======================================================================= #
-    i.delete!(N)   if i.include? N
-    i.delete!('_') if i.include? '_'
-    # ======================================================================= #
-    # Is faster than .gsub()# if i
-    #   i.upcase! # Since as of Feb 2016. Nope, disabled again as of May 2016.
-    # end
-    # ======================================================================= #
-    if File.exist?(i) # Assume a file was given. Read its content then.
-      if i.end_with? '.yml'
-        i = YAML.load_file(i)
-      else
-        i = File.readlines(i).reject {|line|
-          line.start_with? '>'
-        }.join
-      end
-    end
-    # ======================================================================= #
-    # Also get rid of numbers.
-    # ======================================================================= #
-    i = i.dup if i.frozen?
-    i.gsub!(/\d/,'')
-    # ======================================================================= #
-    # The DNA Sequence will be stored in the Sequence object directly.
-    # The nil is required for now, because otherwise set_input() will
-    # upcase the input by default.
-    # ======================================================================= #
-    if dna_sequence_object?.nil?
-      _ = Bioroebe::Sequence.new(i.dup) { :is_DNA }
-      @internal_hash[:array_dna_sequence] << _
-    else
-      # ===================================================================== #
-      # Assign it here.
-      # ===================================================================== #
-      dna_sequence_object?.set_input(i.dup, nil) # We can re-use the old @sequence object.
-    end
-    # ======================================================================= #
-    # === :coding_area
-    #
-    # Whenever the DNA sequence is modified like that, we will also reset
-    # the "coding_area" entry in the internal Hash to the initial nil state.
-    # ======================================================================= #
-    @internal_hash[:coding_area] = nil
-    # ======================================================================= #
-    # === Report that we will perform an assignment next
-    # ======================================================================= #
-    if be_verbose # Check for verbosity
-      case mode?
-      when :dna
-        if i.empty?
-          erev 'Nothing was assigned yet.'
-        else
-          erev 'Setting DNA sequence to ('+
-                sfancy(i.size.to_s)+
-                rev+' nucleotides):'
-          show_string
-        end
-      end
-    end
-    # ======================================================================= #
-    # === Set the Xorg Buffer next
-    # ======================================================================= #
-    if @config and @config.respond_to?(:additionally_set_xorg_buffer) and
-       @config.additionally_set_xorg_buffer
-      begin
-        set_xorg_buffer(i)
-      rescue NoMethodError; end
-    end
-    unless i.empty? # We do not save empty nucleotide Strings.
-      result = i.to_s.dup
-      # ===================================================================== #
-      # We will save into a file next.
-      # ===================================================================== #
-      this_file = "#{log_dir?}dna_string.yml"
-      dna_sequence_object?.save_sequence_to_this_file(this_file) # Save it here.
-      # ===================================================================== #
-      # Since as of Jun 2016, we will also create a new RNA string.
-      #
-      # This has been disabled as of December 2021. It was too strange
-      # to automatically create a RNA string.
-      # ===================================================================== #
-      # if @internal_hash[:array_rna_sequence]
-      #   @internal_hash[:array_rna_sequence].last.set_sequence(result)
-      # end
-      # ===================================================================== #
-      # === Automatically assign the aminoacid sequence as well
-      # ===================================================================== #
-      _ = translate_dna_into_aminoacid(result)
-      set_aminoacids(_, nil, :be_quiet)
-      return result
-    end
-  end; alias set_dna_string           set_dna_sequence # === set_dna_string
-       alias set_string               set_dna_sequence # === set_dna_sequence
-       alias set_main_sequence        set_dna_sequence # === set_main_sequence (set_main_sequence tag)
-       alias set_sequence             set_dna_sequence # === set_sequence (set_sequence tag)
-       alias set_dna                  set_dna_sequence # === set_dna
-       alias set_assign               set_dna_sequence # === set_dna
-       alias assign_sequence          set_dna_sequence # === assign_sequence
-       alias assign_dna_sequence      set_dna_sequence # === assign_dna_sequence
-       alias assign_this_dna_sequence set_dna_sequence # === assign_this_dna_sequence
-       alias assign                   set_dna_sequence # === assign
-
-  require 'bioroebe/aminoacids/create_random_aminoacids.rb'
-  # ========================================================================= #
-  # === set_aminoacids               (assign protein tag, set aminoacids tag)
-  #
-  # Assign a protein sequence here. The first argument shall be the
-  # aminoacid sequence in question.
-  #
-  # Usage example:
-  #
-  #   assign_protein FLIMVSPTAYHQNKDECWRGX*
-  #
-  # ========================================================================= #
-  def set_aminoacids(
-      i                    = :random,
-      how_many_to_generate = 1000, # Will only be applied if the input is :random or :generate
-      be_verbose           = true
-    )
-    # ======================================================================= #
-    # === Sanitize the third argument next
-    # ======================================================================= #
-    case be_verbose
-    when :be_silent,
-         :be_quiet
-      be_verbose = false
-    end
-    i = :random unless i
-    case i
-    # ======================================================================= #
-    # === :random
-    # ======================================================================= #
-    when :random,
-         :generate 
-      i = Bioroebe.create_random_aminoacids(how_many_to_generate) # Generate them here.
-    end
-    i = i.join if i.is_a? Array
-    if i =~ /^\d+$/ # only numbers
-      i = Bioroebe.create_random_aminoacids(i)
-    end
-    if be_verbose
-      erev "Now assigning aminoacid sequence to: #{sfancy(i.to_s)}"
-    end
-    if i.is_a? String # Convert it into a Sequence object here.
-      i = ::Bioroebe.sequence(i) { :aminoacid }
-    end
-    # ======================================================================= #
-    # Do the assignment next.
-    # ======================================================================= #
-    if @internal_hash[:array_aminoacid_sequence].empty?
-      @internal_hash[:array_aminoacid_sequence] << i
-    else
-      @internal_hash[:array_aminoacid_sequence].pop
-      @internal_hash[:array_aminoacid_sequence] << i
-    end 
-  end; alias assign_aminoacid_sequence set_aminoacids # === assign_aminoacid_sequence
-       alias assign_protein_sequence   set_aminoacids # === assign_protein_sequence
-       alias set_aminoacid_sequence    set_aminoacids # === set_aminoacid_sequence
-
-  # ========================================================================= #
-  # === report_that_a_string_must_be_assigned_first
-  # ========================================================================= #
-  def report_that_a_string_must_be_assigned_first
-    erev 'No sequence has been assigned yet. Please first "'+
-         sfancy('assign')+rev+'" a string.'
-  end
-
-  # ========================================================================= #
-  # === assign_sequence2
-  # ========================================================================= #
-  def assign_sequence2(i)
-    i = i.join(' ') if i.is_a? Array
-    @internal_hash[:sequence2] = Bioroebe::Sequence.new(i)
-  end
-
-end; end