+ def process_hmmscan_datas( hmmscan_datas,
+ in_msa,
+ add_position,
+ add_domain_number,
+ add_species,
+ out_msa,
+ out_msa_singles,
+ out_msa_pairs,
+ out_msa_isolated,
+ min_linker )
+
+ actual_out_of = hmmscan_datas.size
+
+ hmmscan_datas.each_with_index do |hmmscan_data, index|
+ if hmmscan_data.number < ( index + 1 )
+ error_msg = "hmmscan_data.number < ( index + 1 ) (this should not have happened)"
+ raise StandardError, error_msg
+ end
+
+ extract_domain( hmmscan_data.seq_name,
+ index + 1,
+ actual_out_of,
+ hmmscan_data.env_from,
+ hmmscan_data.env_to,
+ in_msa,
+ out_msa,
+ add_position,
+ add_domain_number,
+ add_species )
+
+ if min_linker
+ if actual_out_of == 1
+ extract_domain( hmmscan_data.seq_name,
+ 1,
+ 1,
+ hmmscan_data.env_from,
+ hmmscan_data.env_to,
+ in_msa,
+ out_msa_singles,
+ add_position,
+ add_domain_number,
+ add_species )
+ else
+ first = index == 0
+ last = index == hmmscan_datas.length - 1
+
+ if ( ( first && ( ( hmmscan_datas[ index + 1 ].env_from - hmmscan_data.env_to ) > min_linker) ) ||
+ ( last && ( ( hmmscan_data.env_from - hmmscan_datas[ index - 1 ].env_to ) > min_linker ) ) ||
+ ( !first && !last && ( ( hmmscan_datas[ index + 1 ].env_from - hmmscan_data.env_to ) > min_linker ) &&
+ ( ( hmmscan_data.env_from - hmmscan_datas[ index - 1 ].env_to ) > min_linker ) ) )
+
+ extract_domain( hmmscan_data.seq_name,
+ index + 1,
+ actual_out_of,
+ hmmscan_data.env_from,
+ hmmscan_data.env_to,
+ in_msa,
+ out_msa_isolated,
+ add_position,
+ add_domain_number,
+ add_species )
+
+ elsif !first
+ extract_domain( hmmscan_data.seq_name,
+ index.to_s + "+" + ( index + 1 ).to_s,
+ actual_out_of,
+ hmmscan_datas[ index - 1 ].env_from,
+ hmmscan_data.env_to,
+ in_msa,
+ out_msa_pairs,
+ add_position,
+ add_domain_number,
+ add_species )
+ end
+ end
+ end
+ end # each
+ end # def process_hmmscan_data
+