What criteria are used to determine significant HMMER hits?
0
0
Entering edit mode
9 months ago
O.rka ▴ 710

I'm looking at the code here but not very familiar w/ Ruby: https://github.com/takaram/kofam_scan/blob/master/lib/kofam_scan/result.rb

module KofamScan
  class Result
    extend Autoload

    autoload :WithEvalueThreshold
    autoload :WithThresholdScale
    autoload :WithThresholdScaleAndEvalueThreshold

    def self.create(query_list, threshold_scale: nil, e_value_threshold: nil)
      if threshold_scale && e_value_threshold
        WithThresholdScaleAndEvalueThreshold.new(query_list, threshold_scale, e_value_threshold)
      elsif e_value_threshold
        WithEvalueThreshold.new(query_list, e_value_threshold)
      elsif threshold_scale
        WithThresholdScale.new(query_list, threshold_scale)
      else
        Result.new(query_list)
      end
    end

Here's what ko_list looks like:

knum    threshold   score_type  profile_type    F-measure   nseq    nseq_used   alen    mlen    eff_nseq    re/pos  definition
K00001  357.90  domain  all 0.256163    2367    1915    1975    464 14.05   0.590   alcohol dehydrogenase [EC:1.1.1.1]
K00002  443.17  full    all 0.430391    2376    2273    5878    503 7.51    0.590   alcohol dehydrogenase (NADP+) [EC:1.1.1.2]
K00003  286.37  domain  all 0.945500    6268    5369    3257    782 7.03    0.590   homoserine dehydrogenase [EC:1.1.1.3]
K00004  369.60  domain  trim    0.809403    1572    1320    1364    436 5.41    0.590   (R,R)-butanediol dehydrogenase / meso-butanediol dehydrogenase / diacetyl reductase [EC:1.1.1.4 1.1.1.- 1.1.1.303]
K00005  320.93  full    all 0.984022    1449    1051    682 366 1.99    0.590   glycerol dehydrogenase [EC:1.1.1.6]
K00006  316.47  full    all 0.899202    2118    1971    3274    549 4.64    0.590   glycerol-3-phosphate dehydrogenase (NAD+) [EC:1.1.1.8]
K00007  520.80  full    all 0.998236    358 283 834 469 1.01    0.589   D-arabinitol 4-dehydrogenase [EC:1.1.1.11]
K00008  420.17  full    all 0.500025    3737    3281    3597    524 8.29    0.590   L-iditol 2-dehydrogenase [EC:1.1.1.14]
K00009  147.27  full    trim    0.997067    1556    1192    1443    459 3.28    0.590   mannitol-1-phosphate 5-dehydrogenase [EC:1.1.1.17]

I'm seeing a score threshold but not e-value threshold.

Can anyone elaborate on the scheme used for determining significant hits?

kegg protein kofamscan annotation metagenomics • 354 views
ADD COMMENT

Login before adding your answer.

Traffic: 1599 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6