From: cmzmasek Date: Sat, 25 Feb 2017 01:48:39 +0000 (-0800) Subject: in progress... X-Git-Url: http://source.jalview.org/gitweb/?a=commitdiff_plain;h=5e7bf3deb2ae93f542e0b58e260f2afa728b7210;p=jalview.git in progress... --- diff --git a/forester/ruby/evoruby/lib/evo/io/parser/hmmscan_multi_domain_extractor.rb b/forester/ruby/evoruby/lib/evo/io/parser/hmmscan_multi_domain_extractor.rb index 56bf65f..7ce016d 100644 --- a/forester/ruby/evoruby/lib/evo/io/parser/hmmscan_multi_domain_extractor.rb +++ b/forester/ruby/evoruby/lib/evo/io/parser/hmmscan_multi_domain_extractor.rb @@ -280,6 +280,8 @@ module Evoruby query_seq = nil + concatenated_domains = '' + passed_domains.each do |domain| domain_name = domain.model @@ -301,7 +303,7 @@ module Evoruby raise StandardError, error_msg end - extract_domain( query_seq, + extracted_domain_seq = extract_domain( query_seq, domain_counter[domain_name], passed_domains_counts[domain_name], domain.env_from, @@ -309,6 +311,8 @@ module Evoruby in_msa, out_domain_msas[ domain_name ] ) + concatenated_domains += extracted_domain_seq + if domain.env_from < min_env_from min_env_from = domain.env_from end @@ -335,6 +339,7 @@ module Evoruby out_domain_architecture_msa ) puts passed_domains_counts + puts concatenated_domains true end @@ -427,6 +432,8 @@ module Evoruby end out_msa.add_sequence( seq ) + + seq.get_sequence_as_string end def is_ignorable?( line )