X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=code_ruby%2Flib%2Fmaasha%2Fgenbank.rb;h=c3515ecbe326c03e65abb834b1ff84b5413dc0e2;hb=c4b49c5ce1ed3b46ae37e6ebd73c21d67d6d4810;hp=dcdfbe32593ec5c883f2ae9036e555ebe3963c52;hpb=0ed40d0be2b5f54a778f8158557160efe90c818b;p=biopieces.git

diff --git a/code_ruby/lib/maasha/genbank.rb b/code_ruby/lib/maasha/genbank.rb
index dcdfbe3..c3515ec 100644
--- a/code_ruby/lib/maasha/genbank.rb
+++ b/code_ruby/lib/maasha/genbank.rb
@@ -1,4 +1,4 @@
-# Copyright (C) 2007-2010 Martin A. Hansen.
+# Copyright (C) 2007-2011 Martin A. Hansen.
 
 # This program is free software; you can redistribute it and/or
 # modify it under the terms of the GNU General Public License
@@ -22,6 +22,7 @@
 
 # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
 
+require 'maasha/locator'
 require 'maasha/seq'
 require 'maasha/filesys'
 require 'pp'
@@ -80,7 +81,7 @@ class Genbank < Filesys
     i = @entry.size
     j = i - 1
 
-    while @entry[j] and @entry[j] !~ /^[A-Z]/
+    while @entry[j] and @entry[j] =~ /^\s+\d/
       j -= 1
     end
 
@@ -232,130 +233,4 @@ class GenbankFeatures
   end
 end
 
-
-# Error class for all exceptions to do with Genbank/EMBL/DDBJ feature table locators.
-class LocatorError < StandardError; end
-
-class Locator
-  attr_accessor :locator, :seq, :subseq
-
-  def initialize(locator, seq)
-    @locator = locator
-    @seq     = seq
-    @subseq  = Seq.new(nil, "", "dna")
-    parse_locator(locator)
-  end
-
-  def strand
-    if @locator.match("complement")
-      return "-"
-    else
-      return "+"
-    end
-  end
-
-  def s_beg
-    if @locator =~ /(\d+)/
-      return $1.to_i - 1
-    end
-  end
-
-  def s_end
-    if @locator.reverse =~ /(\d+)/
-      return $1.reverse.to_i - 1
-    end
-  end
-
-  private
-
-  # Method that uses recursion to parse a locator string from a feature
-  # table and fetches the appropriate subsequence. the operators
-  # join(), complement(), and order() are handled.
-  # the locator string is broken into a comma separated lists, and
-  # modified if the parens donnot balance. otherwise the comma separated
-  # list of ranges are stripped from operators, and the subsequence are
-  # fetched and handled according to the operators.
-  # SNP locators are also dealt with (single positions).
-  def parse_locator(locator, join = nil, comp = nil, order = nil)
-    intervals = locator.split(",")
-
-    unless balance_parens?(intervals.first)   # locator includes a join/comp/order of several ranges
-      case locator
-      when /^join\((.*)\)$/
-        locator = $1
-        join     = true
-      when /^complement\((.*)\)$/
-        locator = $1
-        comp     = true
-      when /^order\((.*)\)$/
-        locator = $1
-        order    = true
-      end
-
-      parse_locator(locator, join, comp, order)
-    else
-      intervals.each do |interval|
-        case interval
-        when /^join\((.*)\)$/
-          locator = $1
-          join    = true
-          parse_locator(locator, join, comp, order)
-        when /^complement\((.*)\)$/
-          locator = $1
-          comp    = true
-          parse_locator(locator, join, comp, order)
-        when /^order\((.*)\)$/
-          locator = $1
-          order   = true
-          parse_locator(locator, join, comp, order)
-        when /^[<>]?(\d+)[^\d]+(\d+)$/
-          int_beg = $1.to_i - 1
-          int_end = $2.to_i - 1
-
-          newseq = Seq.new(nil, @seq.seq[int_beg...int_end], "dna")
-
-					unless newseq.seq.nil?
-						newseq.revcomp if comp
-
-						@subseq.seq << (order ? " " + newseq.seq : newseq.seq)
-					end
-        when /^(\d+)$/
-          pos = $1.to_i - 1
-
-          newseq = Seq.new(nil, @seq.seq[pos], "dna")
-
-					unless newseq.seq.nil?
-          	newseq.revcomp if comp 
-
-          	@subseq.seq << (order ? " " + newseq.seq : newseq.seq)
-					end
-        else
-          $stderr.puts "WARNING: Could not match locator -> #{locator}";
-          @subseq.seq << ""
-        end
-      end
-    end
-
-    return @subseq
-  end
-
-  def balance_parens?(locator)
-    parens = 0
-
-    locator.each_char do |char|
-      case char
-      when '(' then parens += 1
-      when ')' then parens -= 1
-      end
-    end
-
-    if parens == 0
-      return true
-    else
-      return false
-    end
-  end
-end
-
-
 __END__