jandot-bio 1.2.1

data/lib/bio/appl/bl2seq/report.rb

@@ -0,0 +1,334 @@
+#
+# = bio/appl/bl2seq/report.rb - bl2seq (BLAST 2 sequences) parser
+# 
+# Copyright:: Copyright (C) 2005 GOTO Naohisa <ng@bioruby.org>
+# License::   The Ruby License
+#
+#  $Id: report.rb,v 1.8 2007/04/05 23:35:39 trevor Exp $
+#
+# Bio::Bl2seq::Report is a NCBI bl2seq (BLAST 2 sequences) output parser.
+#
+# = Acknowledgements
+#
+# Thanks to Tomoaki NISHIYAMA <tomoakin __at__ kenroku.kanazawa-u.ac.jp> 
+# for providing bl2seq parser patches based on
+# lib/bio/appl/blast/format0.rb.
+#
+
+require 'bio/appl/blast/format0'
+
+module Bio
+class Blast
+
+  class Bl2seq
+
+    # Bio::Bl2seq::Report is a NCBI bl2seq (BLAST 2 sequences) output parser.
+    # It inherits Bio::Blast::Default::Report.
+    # Most of its methods are the same as Bio::Blast::Default::Report,
+    # but it lacks many methods.
+    class Report < Bio::Blast::Default::Report
+
+      # Delimiter of each entry. Bio::FlatFile uses it.
+      # In Bio::Bl2seq::Report, it it nil (1 entry 1 file).
+      DELIMITER = RS = nil
+
+      undef format0_parse_header
+      undef program, version, version_number, version_date,
+        message, converged?, reference, db
+
+      # Splits headers.
+      def format0_split_headers(data)
+        @f0query = data.shift
+      end
+      private :format0_split_headers
+
+      # Splits the search results.
+      def format0_split_search(data)
+        iterations = []
+        while r = data[0] and /^\>/ =~ r
+          iterations << Iteration.new(data)
+        end
+        if iterations.size <= 0 then
+          iterations << Iteration.new(data)
+        end
+        iterations
+      end
+      private :format0_split_search
+
+      # Stores format0 database statistics.
+      # Internal use only. Users must not use the class.
+      class F0dbstat < Bio::Blast::Default::Report::F0dbstat #:nodoc:
+        # Returns number of sequences in database.
+        def db_num
+          unless defined?(@db_num)
+            parse_params
+            @db_num = @hash['Number of Sequences'].to_i
+          end
+          @db_num
+        end
+
+        # Returns number of letters in database.
+        def db_len
+          unless defined?(@db_len)
+            parse_params
+            @db_len = @hash['length of database'].to_i
+          end
+          @db_len
+        end
+      end #class F0dbstat
+
+      # Bio::Bl2seq::Report::Iteration stores information about
+      # a iteration.
+      # Normally, it may contain some Bio::Bl2seq::Report::Hit objects.
+      #
+      # Note that its main existance reason is to keep complatibility
+      # between Bio::Blast::Default::Report::* classes.
+      class Iteration < Bio::Blast::Default::Report::Iteration
+        # Creates a new Iteration object.
+        # It is designed to be called only internally from
+        # the Bio::Blast::Default::Report class.
+        # Users shall not use the method directly.
+        def initialize(data)
+          @f0stat = []
+          @f0dbstat = Bio::Blast::Default::Report::AlwaysNil.instance
+          @hits = []
+          @num = 1
+          while r = data[0] and /^\>/ =~ r
+            @hits << Hit.new(data)
+          end
+        end
+
+        # Returns the hits of the iteration.
+        # It returns an array of Bio::Bl2seq::Report::Hit objects.
+        def hits; @hits; end
+
+        undef message, pattern_in_database, 
+          pattern, pattern_positions, hits_found_again,
+          hits_newly_found, hits_for_pattern, parse_hitlist,
+          converged?
+      end #class Iteration
+
+      # Bio::Bl2seq::Report::Hit contains information about a hit.
+      # It may contain some Bio::Blast::Default::Report::HSP objects.
+      # All methods are the same as Bio::Blast::Default::Report::Hit class.
+      # Please refer to Bio::Blast::Default::Report::Hit.
+      class Hit < Bio::Blast::Default::Report::Hit
+      end #class Hit
+
+      # Bio::Bl2seq::Report::HSP holds information about the hsp
+      # (high-scoring segment pair).
+      # NOTE that the HSP class below is NOT used because
+      # Ruby's constants namespace are normally statically determined
+      # and HSP object is created in Bio::Blast::Default::Report::Hit class.
+      # Please refer to Bio::Blast::Default::Report::HSP.
+      class HSP < Bio::Blast::Default::Report::HSP
+      end #class HSP
+
+    end #class Report
+  end #class Bl2seq
+
+end #class Blast
+end #module Bio
+
+######################################################################
+
+if __FILE__ == $0
+
+  Bio::FlatFile.open(Bio::Blast::Bl2seq::Report, ARGF) do |ff|
+  ff.each do |rep|
+
+  print "# === Bio::Blast::Bl2seq::Report\n"
+  puts
+  #@#print "  rep.program           #=> "; p rep.program
+  #@#print "  rep.version           #=> "; p rep.version
+  #@#print "  rep.reference         #=> "; p rep.reference
+  #@#print "  rep.db                #=> "; p rep.db
+  #print "  rep.query_id          #=> "; p rep.query_id
+  print "  rep.query_def         #=> "; p rep.query_def
+  print "  rep.query_len         #=> "; p rep.query_len
+  #puts
+  #@#print "  rep.version_number    #=> "; p rep.version_number
+  #@#print "  rep.version_date      #=> "; p rep.version_date
+  puts
+
+  print "# === Parameters\n"
+  #puts
+  #print "  rep.parameters        #=> "; p rep.parameters
+  puts
+  print "  rep.matrix            #=> "; p rep.matrix
+  print "  rep.expect            #=> "; p rep.expect
+  #print "  rep.inclusion         #=> "; p rep.inclusion
+  print "  rep.sc_match          #=> "; p rep.sc_match
+  print "  rep.sc_mismatch       #=> "; p rep.sc_mismatch
+  print "  rep.gap_open          #=> "; p rep.gap_open
+  print "  rep.gap_extend        #=> "; p rep.gap_extend
+  #print "  rep.filter            #=> "; p rep.filter
+  #@#print "  rep.pattern           #=> "; p rep.pattern
+  #print "  rep.entrez_query      #=> "; p rep.entrez_query
+  #puts
+  #@#print "  rep.pattern_positions  #=> "; p rep.pattern_positions
+  puts
+
+  print "# === Statistics (last iteration's)\n"
+  #puts
+  #print "  rep.statistics        #=> "; p rep.statistics
+  puts
+  print "  rep.db_num            #=> "; p rep.db_num
+  print "  rep.db_len            #=> "; p rep.db_len
+  #print "  rep.hsp_len           #=> "; p rep.hsp_len
+  print "  rep.eff_space         #=> "; p rep.eff_space
+  print "  rep.kappa             #=> "; p rep.kappa
+  print "  rep.lambda            #=> "; p rep.lambda
+  print "  rep.entropy           #=> "; p rep.entropy
+  puts
+  print "  rep.num_hits          #=> "; p rep.num_hits
+  print "  rep.gapped_kappa      #=> "; p rep.gapped_kappa
+  print "  rep.gapped_lambda     #=> "; p rep.gapped_lambda
+  print "  rep.gapped_entropy    #=> "; p rep.gapped_entropy
+  print "  rep.posted_date       #=> "; p rep.posted_date
+  puts
+
+  #@#print "# === Message (last iteration's)\n"
+  #@#puts
+  #@#print "  rep.message           #=> "; p rep.message
+  #puts
+  #@#print "  rep.converged?        #=> "; p rep.converged?
+  #@#puts
+
+  print "# === Iterations\n"
+  puts
+  print "  rep.itrerations.each do |itr|\n"
+  puts
+
+  rep.iterations.each do |itr|
+      
+  print "# --- Bio::Blast::Bl2seq::Report::Iteration\n"
+  puts
+
+  print "    itr.num             #=> "; p itr.num
+  #print "    itr.statistics      #=> "; p itr.statistics
+  #@#print "    itr.message         #=> "; p itr.message
+  print "    itr.hits.size       #=> "; p itr.hits.size
+  #puts
+  #@#print "    itr.hits_newly_found.size    #=> "; p itr.hits_newly_found.size;
+  #@#print "    itr.hits_found_again.size    #=> "; p itr.hits_found_again.size;
+  #@#if itr.hits_for_pattern then
+  #@#itr.hits_for_pattern.each_with_index do |hp, hpi|
+  #@#print "    itr.hits_for_pattern[#{hpi}].size #=> "; p hp.size;
+  #@#end
+  #@#end
+  #@#print "    itr.converged?      #=> "; p itr.converged?
+  puts
+
+  print "    itr.hits.each do |hit|\n"
+  puts
+
+  itr.hits.each_with_index do |hit, i|
+
+  print "# --- Bio::Blast::Bl2seq::Default::Report::Hit"
+  print " ([#{i}])\n"
+  puts
+
+  #print "      hit.num           #=> "; p hit.num
+  #print "      hit.hit_id        #=> "; p hit.hit_id
+  print "      hit.len           #=> "; p hit.len
+  print "      hit.definition    #=> "; p hit.definition
+  #print "      hit.accession     #=> "; p hit.accession
+  #puts
+  print "      hit.found_again?  #=> "; p hit.found_again?
+
+  print "        --- compatible/shortcut ---\n"
+  #print "      hit.query_id      #=> "; p hit.query_id
+  #print "      hit.query_def     #=> "; p hit.query_def
+  #print "      hit.query_len     #=> "; p hit.query_len
+  #print "      hit.target_id     #=> "; p hit.target_id
+  print "      hit.target_def    #=> "; p hit.target_def
+  print "      hit.target_len    #=> "; p hit.target_len
+
+  print "            --- first HSP's values (shortcut) ---\n"
+  print "      hit.evalue        #=> "; p hit.evalue
+  print "      hit.bit_score     #=> "; p hit.bit_score
+  print "      hit.identity      #=> "; p hit.identity
+  #print "      hit.overlap       #=> "; p hit.overlap
+
+  print "      hit.query_seq     #=> "; p hit.query_seq
+  print "      hit.midline       #=> "; p hit.midline
+  print "      hit.target_seq    #=> "; p hit.target_seq
+
+  print "      hit.query_start   #=> "; p hit.query_start
+  print "      hit.query_end     #=> "; p hit.query_end
+  print "      hit.target_start  #=> "; p hit.target_start
+  print "      hit.target_end    #=> "; p hit.target_end
+  print "      hit.lap_at        #=> "; p hit.lap_at
+  print "            --- first HSP's vaules (shortcut) ---\n"
+  print "        --- compatible/shortcut ---\n"
+
+  puts
+  print "      hit.hsps.size     #=> "; p hit.hsps.size
+  if hit.hsps.size == 0 then
+  puts  "          (HSP not found: please see blastall's -b and -v options)"
+  puts
+  else
+
+  puts
+  print "      hit.hsps.each do |hsp|\n"
+  puts
+
+  hit.hsps.each_with_index do |hsp, j|
+
+  print "# --- Bio::Blast::Default::Report::HSP (Bio::Blast::Bl2seq::Report::HSP)"
+  print " ([#{j}])\n"
+  puts
+  #print "        hsp.num         #=> "; p hsp.num
+  print "        hsp.bit_score   #=> "; p hsp.bit_score 
+  print "        hsp.score       #=> "; p hsp.score
+  print "        hsp.evalue      #=> "; p hsp.evalue
+  print "        hsp.identity    #=> "; p hsp.identity
+  print "        hsp.gaps        #=> "; p hsp.gaps
+  print "        hsp.positive    #=> "; p hsp.positive
+  print "        hsp.align_len   #=> "; p hsp.align_len
+  #print "        hsp.density     #=> "; p hsp.density
+
+  print "        hsp.query_frame #=> "; p hsp.query_frame
+  print "        hsp.query_from  #=> "; p hsp.query_from
+  print "        hsp.query_to    #=> "; p hsp.query_to
+
+  print "        hsp.hit_frame   #=> "; p hsp.hit_frame
+  print "        hsp.hit_from    #=> "; p hsp.hit_from
+  print "        hsp.hit_to      #=> "; p hsp.hit_to
+
+  #print "        hsp.pattern_from#=> "; p hsp.pattern_from
+  #print "        hsp.pattern_to  #=> "; p hsp.pattern_to
+
+  print "        hsp.qseq        #=> "; p hsp.qseq
+  print "        hsp.midline     #=> "; p hsp.midline
+  print "        hsp.hseq        #=> "; p hsp.hseq
+  puts
+  print "        hsp.percent_identity  #=> "; p hsp.percent_identity
+  #print "        hsp.mismatch_count    #=> "; p hsp.mismatch_count
+  #
+  print "        hsp.query_strand      #=> "; p hsp.query_strand
+  print "        hsp.hit_strand        #=> "; p hsp.hit_strand
+  print "        hsp.percent_positive  #=> "; p hsp.percent_positive
+  print "        hsp.percent_gaps      #=> "; p hsp.percent_gaps
+  puts
+
+  end #each
+  end #if hit.hsps.size == 0
+  end
+  end
+  end #ff.each
+  end #FlatFile.open
+
+end #if __FILE__ == $0
+
+######################################################################
+
+=begin
+
+= Bio::Blast::Bl2seq::Report
+
+    NCBI bl2seq (BLAST 2 sequences) output parser
+
+=end
+

data/lib/bio/appl/blast.rb

@@ -0,0 +1,351 @@
+#
+# = bio/appl/blast.rb - BLAST wrapper
+# 
+# Copyright::  Copyright (C) 2001,2008  Mitsuteru C. Nakao <n@bioruby.org>
+# Copyright::  Copyright (C) 2002,2003  Toshiaki Katayama <k@bioruby.org>
+# Copyright::  Copyright (C) 2006       Jan Aerts <jan.aerts@bbsrc.ac.uk>
+# License::    The Ruby License
+#
+# $Id: blast.rb,v 1.35 2008/04/15 13:54:39 ngoto Exp $
+#
+
+require 'net/http'
+require 'cgi' unless defined?(CGI)
+require 'bio/command'
+require 'shellwords'
+
+module Bio
+
+  # == Description
+  # 
+  # The Bio::Blast class contains methods for running local or remote BLAST
+  # searches, as well as for parsing of the output of such BLASTs (i.e. the
+  # BLAST reports). For more information on similarity searches and the BLAST
+  # program, see http://www.ncbi.nlm.nih.gov/Education/BLASTinfo/similarity.html.
+  #
+  # == Usage
+  #
+  #   require 'bio'
+  #   
+  #   # To run an actual BLAST analysis:
+  #   #   1. create a BLAST factory
+  #   remote_blast_factory = Bio::Blast.remote('blastp', 'SWISS',
+  #                                            '-e 0.0001', 'genomenet')
+  #   #or:
+  #   local_blast_factory = Bio::Blast.local('blastn','/path/to/db')
+  #
+  #   #   2. run the actual BLAST by querying the factory
+  #   report = remote_blast_factory.query(sequence_text)
+  #
+  #   # Then, to parse the report, see Bio::Blast::Report
+  #
+  # === Available databases for Bio::Blast.remote
+  #
+  #  ----------+-------+---------------------------------------------------
+  #   program  | query | db (supported in GenomeNet)
+  #  ----------+-------+---------------------------------------------------
+  #   blastp   | AA    | nr-aa, genes, vgenes.pep, swissprot, swissprot-upd,
+  #  ----------+-------+ pir, prf, pdbstr
+  #   blastx   | NA    | 
+  #  ----------+-------+---------------------------------------------------
+  #   blastn   | NA    | nr-nt, genbank-nonst, gbnonst-upd, dbest, dbgss,
+  #  ----------+-------+ htgs, dbsts, embl-nonst, embnonst-upd, epd,
+  #   tblastn  | AA    | genes-nt, genome, vgenes.nuc
+  #  ----------+-------+---------------------------------------------------
+  #
+  # == See also
+  #
+  # * Bio::Blast::Report
+  # * Bio::Blast::Report::Hit
+  # * Bio::Blast::Report::Hsp
+  #
+  # == References
+  # 
+  # * http://www.ncbi.nlm.nih.gov/blast/
+  # * http://www.ncbi.nlm.nih.gov/Education/BLASTinfo/similarity.html
+  # * http://blast.genome.jp/ideas/ideas.html#blast
+  #
+  class Blast
+
+    autoload :Fastacmd,     'bio/io/fastacmd'
+    autoload :Report,       'bio/appl/blast/report'
+    autoload :Default,      'bio/appl/blast/format0'
+    autoload :WU,           'bio/appl/blast/wublast'
+    autoload :Bl2seq,       'bio/appl/bl2seq/report'
+    autoload :RPSBlast,     'bio/appl/blast/rpsblast'
+
+    # This is a shortcut for Bio::Blast.new:
+    #  Bio::Blast.local(program, database, options)
+    # is equivalent to
+    #  Bio::Blast.new(program, database, options, 'local')
+    # ---
+    # *Arguments*:
+    # * _program_ (required): 'blastn', 'blastp', 'blastx', 'tblastn' or 'tblastx'
+    # * _db_ (required): name of the local database
+    # * _options_: blastall options \
+    # (see http://www.genome.jp/dbget-bin/show_man?blast2)
+    # *Returns*:: Bio::Blast factory object
+    def self.local(program, db, option = '')
+      self.new(program, db, option, 'local')
+    end
+
+    # Bio::Blast.remote does exactly the same as Bio::Blast.new, but sets
+    # the remote server 'genomenet' as its default.
+    # ---
+    # *Arguments*:
+    # * _program_ (required): 'blastn', 'blastp', 'blastx', 'tblastn' or 'tblastx'
+    # * _db_ (required): name of the remote database
+    # * _options_: blastall options \
+    # (see http://www.genome.jp/dbget-bin/show_man?blast2)
+    # * _server_: server to use (DEFAULT = 'genomenet')
+    # *Returns*:: Bio::Blast factory object
+    def self.remote(program, db, option = '', server = 'genomenet')
+      self.new(program, db, option, server)
+    end
+
+    # the method Bio::Blast.report is moved from bio/appl/blast/report.rb.
+    # only for xml format
+    def self.reports(input, parser = nil)
+      ary = []
+      input.each("</BlastOutput>\n") do |xml|
+        xml.sub!(/[^<]*(<?)/, '\1') # skip before <?xml> tag
+        next if xml.empty?          # skip trailing no hits
+        if block_given?
+          yield Report.new(xml, parser)
+        else
+          ary << Report.new(xml, parser)
+        end
+      end
+      return ary
+    end
+
+
+    # Program name (_-p_ option for blastall): blastp, blastn, blastx, tblastn
+    # or tblastx
+    attr_accessor :program
+
+    # Database name (_-d_ option for blastall)
+    attr_accessor :db
+    
+    # Options for blastall
+    attr_accessor :options
+
+    # Server to submit the BLASTs to
+    attr_accessor :server
+
+    # Full path for blastall. (default: 'blastall').
+    attr_accessor :blastall
+
+    # Substitution matrix for blastall -M
+    attr_accessor :matrix
+
+    # Filter option for blastall -F (T or F).
+    attr_accessor :filter
+
+    # Returns a String containing blast execution output in as is the Bio::Blast#format.
+    attr_reader :output
+
+    # Output report format for blastall -m 
+    #
+    # 0, pairwise; 1; 2; 3; 4; 5; 6; 7, XML Blast outpu;, 8, tabular; 
+    # 9, tabular with comment lines; 10, ASN text; 11, ASN binery [intege].
+    attr_reader :format
+
+    #
+    attr_writer :parser  # to change :xmlparser, :rexml, :tab
+
+
+    # Creates a Bio::Blast factory object.
+    # 
+    # To run any BLAST searches, a factory has to be created that describes a
+    # certain BLAST pipeline: the program to use, the database to search, any
+    # options and the server to use. E.g.
+    # 
+    #   blast_factory = Bio::Blast.new('blastn','dbsts', '-e 0.0001 -r 4', 'genomenet')
+    # 
+    # ---
+    # *Arguments*:
+    # * _program_ (required): 'blastn', 'blastp', 'blastx', 'tblastn' or 'tblastx'
+    # * _db_ (required): name of the (local or remote) database
+    # * _options_: blastall options \
+    # (see http://www.genome.jp/dbget-bin/show_man?blast2)
+    # * _server_: server to use (e.g. 'genomenet'; DEFAULT = 'local')
+    # *Returns*:: Bio::Blast factory object
+    def initialize(program, db, opt = [], server = 'local')
+      @program  = program
+      @db       = db
+      @server   = server
+
+      @blastall = 'blastall'
+      @matrix   = nil
+      @filter   = nil
+
+      @output   = ''
+      @parser   = nil
+      @format   = 0
+
+      set_options(opt)
+    end      
+
+
+    # This method submits a sequence to a BLAST factory, which performs the
+    # actual BLAST.
+    # 
+    #   fasta_sequences = Bio::FlatFile.open(Bio::FastaFormat, 'my_sequences.fa')
+    #   report = blast_factory.query(fasta_sequences)
+    # 
+    # ---
+    # *Arguments*:
+    # * _query_ (required): single- or multiple-FASTA formatted sequence(s)
+    # *Returns*:: a Bio::Blast::Report object
+    def query(query)
+      return self.send("exec_#{@server}", query.to_s)
+    end
+
+    # Returns options of blastall
+    def option
+      # backward compatibility
+      Bio::Command.make_command_line(@options)
+    end
+
+    # Set options for blastall
+    def option=(str)
+      # backward compatibility
+      @options = Shellwords.shellwords(str)
+    end
+
+    private
+
+    def set_options(opt = nil)
+      opt = @options unless opt
+      begin
+        a = opt.to_ary
+      rescue NameError #NoMethodError
+        # backward compatibility
+        a = Shellwords.shellwords(opt)
+      end
+      unless a.find { |x| /\A\-m/ =~ x.to_s } then
+        if defined?(XMLParser) or defined?(REXML)
+          @format = 7
+        else
+          @format = 8
+        end
+      else
+        @format = a[a.index('-m') + 1].to_i
+      end
+      @options = [ *a ]
+    end
+
+
+    def parse_result(data)
+      Report.new(data, @parser)
+    end
+
+
+    def make_command_line
+      set_options
+      cmd = [ @blastall, '-p', @program, '-d', @db ]
+      if @matrix
+        cmd.concat([ '-M', @matrix ]) 
+        i = @options.index('-M')
+        @options.delete_at(i)
+        @options.delete_at(i)
+      end
+      if @filter
+        cmd.concat([ '-F', @filter ]) 
+        i = @options.index('-F')
+        @options.delete_at(i)
+        @options.delete_at(i)
+      end
+      if @format
+        cmd.concat([ '-m', @format.to_s ])
+        i = @options.index('-m')
+        @options.delete_at(i)
+        @options.delete_at(i)
+      end
+      cmd.concat(@options) if @options
+    end
+
+
+    def exec_local(query)
+      cmd = make_command_line
+      report = nil
+
+      @output = Bio::Command.query_command(cmd, query)
+      report = parse_result(@output)
+
+      return report
+    end
+
+
+    def exec_genomenet_tab(query)
+      @format = 8
+      exec_genomenet(query)
+    end
+
+
+    def exec_genomenet(query)
+      host = "blast.genome.jp"
+      #path = "/sit-bin/nph-blast"
+      path = "/sit-bin/blast" #2005.08.12
+
+      matrix = @matrix ? @matrix : 'blosum62'
+      filter = @filter ? @filter : 'T'
+
+      opt = []
+      opt.concat([ '-m', @format.to_s ]) if @format
+      opt.concat(@options) if @options
+
+      form = {
+        'style'          => 'raw',
+        'prog'           => @program,
+        'dbname'         => @db,
+        'sequence'       => CGI.escape(query),
+        'other_param'    => CGI.escape(Bio::Command.make_command_line_unix(opt)),
+        'matrix'         => matrix,
+        'filter'         => filter,
+        'V_value'        => 500, # default value for GenomeNet
+        'B_value'        => 250, # default value for GenomeNet
+        'alignment_view' => 0,
+      }
+
+      data = []
+
+      form.each do |k, v|
+        data.push("#{k}=#{v}") if v
+      end
+
+      report = nil
+
+      begin
+        http = Bio::Command.new_http(host)
+        http.open_timeout = 300
+        http.read_timeout = 600
+        result, = http.post(path, data.join('&'))
+        @output = result.body
+        # workaround 2005.08.12
+        if /\<A +HREF=\"(http\:\/\/blast\.genome\.jp(\/tmp\/[^\"]+))\"\>Show all result\<\/A\>/i =~ @output.to_s then
+          result, = http.get($2)
+          @output = result.body
+          txt = @output.to_s.split(/\<pre\>/)[1]
+          raise 'cannot understand response' unless txt
+          txt.sub!(/\<\/pre\>.*\z/m, '')
+          txt.sub!(/.*^ \-{20,}\s*/m, '')
+          @output = txt.gsub(/\&lt\;/, '<')
+          report = parse_result(@output)
+        else
+          raise 'cannot understand response'
+        end
+      end
+
+      return report
+    end
+
+    def exec_ncbi(query)
+      raise NotImplementedError
+    end
+
+  end # class Blast
+
+end # module Bio
+