+ if ( domain_not_present_counter + passed_seqs.get_number_of_seqs + proteins_with_failing_domains ) != protein_counter
+ error_msg = "not present + passing + not passing != proteins in sequence (fasta) file (this should not have happened)"
+ raise StandardError, error_msg
+ end
+
+ puts
+ log << ld
+
+ log << ld
+ avg_pass = ( passing_target_length_sum / domain_pass_counter )
+ puts( "Passing target domain lengths: average: " + avg_pass.to_s )
+ log << "Passing target domain lengths: average: " + avg_pass.to_s
+ log << ld
+ puts( "Passing target domain lengths: min-max: " + passing_target_length_min.to_s + " - " + passing_target_length_max.to_s)
+ log << "Passing target domain lengths: min-max: " + passing_target_length_min.to_s + " - " + passing_target_length_max.to_s
+ log << ld
+ puts( "Passing target domain iE: min-max: " + passing_target_ie_min.to_s + " - " + passing_target_ie_max.to_s)
+ log << "Passing target domain iE: min-max: " + passing_target_ie_min.to_s + " - " + passing_target_ie_max.to_s
+ log << ld
+ puts( "Passing target domains: sum: " + domain_pass_counter.to_s )
+ log << "Passing target domains: sum: " + domain_pass_counter.to_s
+ log << ld
+ log << ld
+ puts
+ sum = domain_pass_counter + domain_fail_counter
+ avg_all = overall_target_length_sum / sum
+ puts( "All target domain lengths: average: " + avg_all.to_s )
+ log << "All target domain lengths: average: " + avg_all.to_s
+ log << ld
+ puts( "All target domain lengths: min-max: " + overall_target_length_min.to_s + " - " + overall_target_length_max.to_s)
+ log << "All target domain lengths: min-max: " + overall_target_length_min.to_s + " - " + overall_target_length_max.to_s
+ log << ld
+ puts( "All target target domain iE: min-max: " + overall_target_ie_min.to_s + " - " + overall_target_ie_max.to_s)
+ log << "All target target domain iE: min-max: " + overall_target_ie_min.to_s + " - " + overall_target_ie_max.to_s
+ log << ld
+ puts( "All target domains: sum: " + sum.to_s )
+ log << "All target domains: sum: " + sum.to_s
+
+ puts
+ puts( "Proteins with passing target domain(s): " + passed_seqs.get_number_of_seqs.to_s )
+ puts( "Proteins with no passing target domain: " + proteins_with_failing_domains.to_s )
+ puts( "Proteins with no target domain : " + domain_not_present_counter.to_s )
+
+ log << ld