in progress...
authorcmzmasek <chris.zma@outlook.com>
Sat, 25 Feb 2017 01:48:39 +0000 (17:48 -0800)
committercmzmasek <chris.zma@outlook.com>
Sat, 25 Feb 2017 01:48:39 +0000 (17:48 -0800)
forester/ruby/evoruby/lib/evo/io/parser/hmmscan_multi_domain_extractor.rb

index 56bf65f..7ce016d 100644 (file)
@@ -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 )