2.2g hmmpfam - search sequences against an HMM database S. Eddy hmm hmmpfam hmmpfam perl "hmmpfam" 0 description Description of hmmpfam hmmpfam reads a single sequence from seqfile and compares it against all the HMMs in hmmfile looking for significantly similar sequence matches. The output consists of three sections: a ranked list of the best scoring HMMs, a list of the best scoring domains in order of their occurrence in the sequence, and alignments for all the best scoring domains. A sequence score may be higher than a domain score for the same sequence if there is more than one domain in the sequence; the sequence score takes into account all the domains. All sequences scoring above the -E and -T cutoffs are shown in the first list, then every domain found in this list is shown in the second list of domain hits. If desired, E-value and score thresholds may also be applied to the domain list using the --domE and --domT options. toto perl "" 1000 seqfile Sequence file perl " $value" 3 8 seqfile perl 1 HMMDB HMM database perl " $value" Pfam 2 Pfam Pfam PfamFrag PfamFrag nucleic_acid Force nucleic acid models and sequences (-n) perl ($value) ? " -n" : "" 1 n_best number of reported alignments (-A n) perl ($value) ? " -A $value" : "" 1 Limits the alignment output to the n best scoring domains. -A0 shuts off the alignment output and can be used to reduce the size of output files. E_value_cutoff E_value cutoff (-E x) perl (defined $value && $value != $vdef) ? " -E $value" : "" 10.0 1 Set the E-value cutoff for the per-sequence ranked hit list to x, where x is a positive real number. The default is 10.0. Hits with E-values better than (less than) this threshold will be shown. Bit_cutoff Bit score cutoff (-T x) perl ($value) ? " -T $value" : "" 1 Set the bit score cutoff for the per-sequence ranked hit list to x, where x is a real number. The default is negative infinity; by default, the threshold is controlled by E-value and not by bit score. Hits with bit scores better than (greater than) this threshold will be shown. E_value_calculation Control of E_value calculation (-Z n) perl ($value) ? " -Z $value" : "" 1 Calculate the E-value scores as if we had seen a sequence database of n sequences. The default is arbitrarily set to 59021, the size of Swissprot 34. expert_options Expert Options acc Report accessions instead of names in the output reports (--acc) perl ($value) ? " --acc" : "" compat Use the output format of HMMER 2.1.1 (--compat) perl ($value) ? " --compat" : "" 1 domE E-value cutoff for the per-domain ranked hit list (--domE x) perl ($value) ? " --domE $value" : "" 1 Set the E-value cutoff for the per-domain ranked hit list to x, where x is a positive real number. The default is infinity; by default, all domains in the sequences that passed the first threshold will be reported in the second list, so that the number of domains reported in the per-sequence list is consistent with the number that appear in the per-domain list. domT bit score cutoff for the per-domain ranked hit list (--domT x) perl ($value) ? " --domT $value" : "" 1 Set the bit score cutoff for the per-domain ranked hit list to x, where x is a real number. The default is negative infinity; by default, all domains in the sequences that passed the first threshold will be reported in the second list, so that the number of domains reported in the per-sequence list is consistent with the number that appear in the per-domain list. Important note: only one domain in a sequence is absolutely controlled by this parameter, or by --domT. The second and subsequent domains in a sequence have a de facto bit score threshold of 0 because of the details of how HMMER works. HMMER requires at east one pass through the main model per sequence; to do more than one pass (more than one domain) the multidomain alignment must have a better score than the single domain alignment, and hence the extra domains must contribute positive score. See the Users' Guide for more detail. forward forward algorithm (--forward) perl ($value) ? " --forward" : "" 1 Use the Forward algorithm instead of the Viterbi algorithm to determine the per-sequence scores. Per-domain scores are still determined by the Viterbi algorithm. Some have argued that Forward is a more sensitive algorithm for detecting remote sequence homologues; my experiments with HMMER have not confirmed this, however. null2 turns off the second post processing step (--null2) perl ($value)? " --null2" : "" 1 Turn off the post hoc second null model. By default, each alignment is rescored by a postprocessing step that takes into account possible biased composition in either the HMM or the target sequence. This is almost essential in database searches, especially with local alignment models. There is a very small chance that this postprocessing might remove real matches, and in these cases --null2 may improve sensitivity at the expense of reducing specificity by letting biased composition hits through. xnu turns on XNU filtering (--xnu) perl ($value) ? " --xnu" : "" 1 Turn on XNU filtering of target protein sequences. Has no effect on nucleic acid sequences. In trial experiments, --xnu appears to perform less well than the default post hoc null2 model.