end
add_species = false
- if cla.is_option_set? ADD_SPECIES
+ if cla.is_option_set? ADD_SPECIES
add_species = true
end
-
+
if ( add_domain_number_as_letter && add_domain_number_as_digit )
puts( "attempt to add domain number as letter and digit at the same time" )
print_help
add_domain_number_as_digit,
add_domain_number_as_letter,
trim,
- add_species,
+ add_species,
log )
rescue ArgumentError, IOError, StandardError => e
Util.fatal_error( PRG_NAME, "error: " + e.to_s, STDOUT )
add_domain_number_as_digit,
add_domain_number_as_letter,
trim_name,
- add_species,
+ add_species,
log )
Util.check_file_for_readability( hmmsearch_output )
if domain_id != target_name
next
end
-
-
+
+
sequence = $4
number = $10.to_i
out_of = $11.to_i
add_domain_number_as_digit,
add_domain_number_as_letter,
trim_name,
- add_species )
+ add_species )
if ( number < 1 || out_of < 1 || number > out_of )
error_msg = "impossible: number=" + number.to_s + ", out of=" + out_of.to_s
raise ArgumentError, error_msg
end
# hmmsearch is 1 based, wheres sequences are 0 bases in this package.
seq = in_msa.get_sequence( seqs[ 0 ] ).get_subsequence( seq_from - 1, seq_to - 1 )
-
+
orig_name = seq.get_name
-
+
seq.set_name( orig_name.split[ 0 ] )
-
+
if add_position
seq.set_name( seq.get_name + "_" + seq_from.to_s + "-" + seq_to.to_s )
end
elsif ( add_domain_number )
seq.set_name( seq.get_name + "~" + number.to_s + "-" + out_of.to_s )
end
-
+
if add_species
a = orig_name.rindex "["
b = orig_name.rindex "]"
unless a && b
- error_msg = "species not found in " + orig_name
- raise StandardError, error_msg
+ error_msg = "species not found in " + orig_name
+ raise StandardError, error_msg
end
species = orig_name[ a .. b ]
seq.set_name( seq.get_name + " [" + species + "]" )
- end
-
+ end
+
end
# if ( seq.get_name.length > 10 )