# hmmsearch :: search profile(s) against a sequence database # HMMER 3.1b2 (February 2015); http://hmmer.org/ # Copyright (C) 2015 Howard Hughes Medical Institute. # Freely distributed under the GNU General Public License (GPLv3). # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # query HMM file: pfam_list/Methyltransf_10.hmm.txt # target sequence database: proteomes/Trichomonas_vaginalis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Methyltransf_10 [M=299] Accession: PF05971.11 Description: Protein of unknown function (DUF890) Scores for complete sequences (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Sequence Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- ------ inclusion threshold ------ 0.37 11.5 1.2 2.3 8.9 0.1 2.0 2 tr|A2F938|A2F938_TRIVA Uncharacterized protein OS=Trichomona Domain annotation for each sequence (and alignments): >> tr|A2F938|A2F938_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_429190 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 1.7 0.1 0.0073 3.7e+02 148 189 .. 32 76 .. 18 132 .. 0.79 2 ? 8.9 0.1 4.6e-05 2.3 143 211 .. 156 230 .. 143 288 .. 0.74 Alignments for each domain: == domain 1 score: 1.7 bits; conditional E-value: 0.0073 HT.T-TTTEEEEE--ST-SSTTTSTT--S..-EEE.EEE-----S CS Methyltransf_10 148 snsklsdaiklrkqkekkaifdglieene..rydv.tlCnPPfha 189 + sk +i+++ k+ i+++++++ + +++ + C PP + tr|A2F938|A2F938_TRIVA 32 QSSKNILTIRCKELKDLIDIYNSIVDKYQelMFKMpSNCFPPSFP 76 4445556799999***********998652667764569999764 PP == domain 2 score: 8.9 bits; conditional E-value: 4.6e-05 HHHHHHT.T-TTTEEEEE--ST-SSTTTSTT--S-EE.EEEE-----SS------------...---......... CS Methyltransf_10 143 kaivesnsklsdaiklrkqkekkaifdglieeneryd.vtlCnPPfhaslaeakagserklknlgkank......t 211 ++ ++n+ +++aik k +++++ +gl + y+ + lC f +++++ ++ s r+ n+ ++ t tr|A2F938|A2F938_TRIVA 156 QKLFKTNTSINTAIKNEKIVNQEMVSKGLDLYRNTYKaFKLCKDFFSSQKKDNNENSAREF-NIEDLKIalslvnT 230 678899************************9999997256787666666666666666766.77554223333330 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (299 nodes) Target sequences: 50190 (17055200 residues searched) Passed MSV filter: 1005 (0.0200239); expected 1003.8 (0.02) Passed bias filter: 637 (0.0126918); expected 1003.8 (0.02) Passed Vit filter: 46 (0.000916517); expected 50.2 (0.001) Passed Fwd filter: 1 (1.99243e-05); expected 0.5 (1e-05) Initial search space (Z): 50190 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.32u 0.01s 00:00:00.33 Elapsed: 00:00:00.12 # Mc/sec: 42495.87 // [ok]