lib/genevalidator/hsp.rb in genevalidator-1.6.1 vs lib/genevalidator/hsp.rb in genevalidator-1.6.2
- old
+ new
@@ -31,12 +31,11 @@
# Initializes the corresponding attribute of the hsp
# with respect to the column name of the tabular blast output
# Params:
# +column+: String with column name.
# +value+: Value of the column
- # +type+: type of the sequences: :nucleotide or :protein
- def init_tabular_attribute(hash, _type = :protein)
+ def init_tabular_attribute(hash)
@match_query_from = hash['qstart'].to_i if hash['qstart']
@match_query_to = hash['qend'].to_i if hash['qend']
@query_reading_frame = hash['qframe'].to_i if hash['qframe']
@hit_from = hash['sstart'].to_i if hash['sstart']
@hit_to = hash['send'].to_i if hash['send']
@@ -45,15 +44,13 @@
@align_len = hash['length'].to_i if hash['length']
@pidentity = hash['pident'].to_f if hash['pident']
@identity = hash['nident'].to_f if hash['nident']
@hsp_evalue = hash['evalue'].to_f if hash['evalue']
if hash['qseq']
- puts @query_alignment
query_seq_type = BlastUtils.guess_sequence_type(@query_alignment)
fail SequenceTypeError if query_seq_type != :protein
end
if hash['sseq']
- puts @hit_alignment
hit_seq_type = BlastUtils.guess_sequence_type(@hit_alignment)
fail SequenceTypeError if hit_seq_type != :protein
end
end
end