X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;ds=sidebyside;f=forester%2Fruby%2Fevoruby%2Flib%2Fevo%2Ftool%2Fhmmscan_summary.rb;h=92bceb2a045c0d9cd212ad0fd99521e3cfaf858a;hb=bb5057ba057d5b76362e1c27af04d4ec04feb753;hp=e08ec29264ead36af351f2d827402fd20160d34c;hpb=31a4f995796b7e872e89ad8f07236aef4576f266;p=jalview.git diff --git a/forester/ruby/evoruby/lib/evo/tool/hmmscan_summary.rb b/forester/ruby/evoruby/lib/evo/tool/hmmscan_summary.rb index e08ec29..92bceb2 100644 --- a/forester/ruby/evoruby/lib/evo/tool/hmmscan_summary.rb +++ b/forester/ruby/evoruby/lib/evo/tool/hmmscan_summary.rb @@ -320,6 +320,7 @@ module Evoruby 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, @@ -359,6 +360,8 @@ module Evoruby i_e_value_threshold, uniprotkb, species ) + + dc = 0 # filter according to i-Evalue threshold @@ -366,6 +369,9 @@ module Evoruby hmmscan_results_per_protein_filtered = [] hmmscan_results_per_protein.each do | r | + + puts r + if r.model == hmm_for_protein_output if r.fs_e_value > fs_e_value_threshold return @@ -404,12 +410,13 @@ module Evoruby 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 )