This is a test version of Biostars. For the public version, visit https://www.biostars.org.
What criteria are used to determine significant HMMER hits?

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

0 answers

No answers yet.

Log in to answer this question.