Skip to content

Instantly share code, notes, and snippets.

@chasemc
Last active May 16, 2023 12:54
Show Gist options
  • Save chasemc/270ec6eaaa864269b9bed2f12816d55f to your computer and use it in GitHub Desktop.
Save chasemc/270ec6eaaa864269b9bed2f12816d55f to your computer and use it in GitHub Desktop.
hmmsearch a model with --cut_ga when the model has no defined cutoffs
hmmsearch --domtblout "a.dom" -Z 100 --cut_ga '/socialgene_nr_hmms_file_without_cutoffs_1_of_1.hmm' /example.fa
# hmmsearch :: search profile(s) against a sequence database
# HMMER 3.3.2 (Nov 2020); http://hmmer.org/
# Copyright (C) 2020 Howard Hughes Medical Institute.
# Freely distributed under the BSD open source license.
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# query HMM file: /socialgene_nr_hmms_file_without_cutoffs_1_of_1.hmm
# target sequence database: /example.fa
# per-dom hits tabular output: a.dom
# model-specific thresholding: GA cutoffs
# sequence search space set to: 100
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: DLq3VAeNtnYBFERw4zkJh1S_KMGoVfs5 [M=97]
Error: GA bit thresholds unavailable on model DLq3VAeNtnYBFERw4zkJh1S_KMGoVfs5
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment