fs_e_value_threshold,
hmm_for_protein_output,
i_e_value_threshold,
- true )
+ false,
+ species )
end
hmmscan_results_per_protein.clear
end
end
end
end
+
if !hmm_for_protein_output.empty? && !hmmscan_results_per_protein.empty?
process_hmmscan_results_per_protein( hmmscan_results_per_protein,
fs_e_value_threshold,
hmm_for_protein_output,
i_e_value_threshold,
- true )
+ false,
+ species )
end
outfile.flush()
fs_e_value_threshold,
hmm_for_protein_output,
i_e_value_threshold,
- uniprotkb )
+ uniprotkb,
+ species )
+
+
dc = 0
# filter according to i-Evalue threshold
hmmscan_results_per_protein_filtered = []
hmmscan_results_per_protein.each do | r |
+
+ puts r.model
+ puts r.fs_e_value
+
if r.model == hmm_for_protein_output
if r.fs_e_value > fs_e_value_threshold
return
s << r.model + " "
end
s << "\t"
- e = UniprotKB::get_entry_by_id( process_id( own.query ) )
+ puts s
+ #e = UniprotKB::get_entry_by_id( process_id( own.query ) )
- if e != nil
- s << uniprot_annotation( e )
- s << "\uniprot_annotationt"
- end
+ #if e != nil
+ # s << uniprot_annotation( e )
+ # # s << "\uniprot_annotationt"
+ #end
overview = make_overview( hmmscan_results_per_protein_filtered, hmm_for_protein_output )