# 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/DREV.hmm.txt # target sequence database: proteomes/Trichomonas_vaginalis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: DREV [M=265] Accession: PF05219.11 Description: DREV methyltransferase 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.063 13.8 0.0 0.12 13.0 0.0 1.4 1 tr|A2EVN6|A2EVN6_TRIVA rRNA adenine N(6)-methyltransferase O 0.065 13.8 0.1 0.084 13.4 0.1 1.3 1 tr|A2EXH2|A2EXH2_TRIVA Uncharacterized protein OS=Trichomona 0.081 13.5 0.0 0.1 13.1 0.0 1.1 1 tr|A2FAI4|A2FAI4_TRIVA Uncharacterized protein OS=Trichomona 0.34 11.4 0.1 0.51 10.9 0.1 1.2 1 tr|A2FAE3|A2FAE3_TRIVA Uncharacterized protein OS=Trichomona Domain annotation for each sequence (and alignments): >> tr|A2EVN6|A2EVN6_TRIVA rRNA adenine N(6)-methyltransferase OS=Trichomonas vaginalis OX=5722 GN=TVAG_299080 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.0 0.0 9.3e-06 0.12 95 135 .. 44 85 .. 21 110 .. 0.82 Alignments for each domain: == domain 1 score: 13.0 bits; conditional E-value: 9.3e-06 DREV 95 dkklldlGaGdGeitkrv.aptfeevyatelsetmrdrlkkk 135 k+l++G G+G++t+ + + ev a e ++m lkkk tr|A2EVN6|A2EVN6_TRIVA 44 GDKILEIGPGNGNMTEEMlSREGIEVIAIEKDQRMCVELKKK 85 579************9983455679**************997 PP >> tr|A2EXH2|A2EXH2_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_253130 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.4 0.1 6.7e-06 0.084 62 154 .. 44 142 .. 33 220 .. 0.70 Alignments for each domain: == domain 1 score: 13.4 bits; conditional E-value: 6.7e-06 DREV 62 nGllgrGsmfvlseeqfrkllaiekespasdqedkklldlGaGdGeitk.rvaptfee........vyatelsetmrdrlkkkkynv 139 nG+lg G v s + + + k++ + ++ +k d+GaG G +++ +++++f+e + + e++++++ k y tr|A2EXH2|A2EXH2_TRIVA 44 NGMLG-GLPQVSSTDAIQSEQFVSKYQLNHGMGCQKCADIGAGIGRVSElILSKYFKEidlvepvqKFVDVAKEKLKNKVILKTYT- 128 67776.56667788888888889999*********************983578999994444444434444555555555555554. PP DREV 140 lsalelqetdvkldl 154 +a +++ d +d tr|A2EXH2|A2EXH2_TRIVA 129 CGAQDWK-IDGVFDC 142 3555665.3434444 PP >> tr|A2FAI4|A2FAI4_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_415010 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.1 0.0 8.1e-06 0.1 79 187 .. 22 136 .. 8 156 .. 0.82 Alignments for each domain: == domain 1 score: 13.1 bits; conditional E-value: 8.1e-06 DREV 79 rkllaiekespasdqedkklldlGaGdGeitkrvaptfee..........vyatelsetmrdrlkkkkynvlsalelqetdvkldli 155 +++ i k + + ++ ++ld+G +G + ++ + ++ + +e s +m + k +y+++ ++ + + k+dl+ tr|A2FAI4|A2FAI4_TRIVA 22 DHFVDIIKSQSNL-NKSMTVLDFGCEKGFLGLKLIDDAKKiafldpdpnaISIAEESLKM---MGKTNYEIINKVVEEYNGSKFDLV 104 5566777667666.7899************999999988888888886666666666666...67899****998777789****** PP DREV 156 lclnlldrcfdplklledirsalepngrviva 187 + l d l i+s l+p+g+ ++ tr|A2FAI4|A2FAI4_TRIVA 105 YASVSLHYVPDVKSALTKIKSLLNPDGKLVII 136 ***************************98874 PP >> tr|A2FAE3|A2FAE3_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_299880 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.9 0.1 4e-05 0.51 90 186 .. 33 136 .. 18 143 .. 0.79 Alignments for each domain: == domain 1 score: 10.9 bits; conditional E-value: 4e-05 DREV 90 asdqedkklldlGaGdGeitkrvaptfee...vyatelset....mrdrlkkkkynvlsalelqetdvkldlilclnlldrcfdplk 169 a +d ++ldlG G G + + ++ + te s++ ++ lk +y+ + + q + k+dl + + +d l tr|A2FAE3|A2FAE3_TRIVA 33 AVLSKDIHVLDLGCGLGCTGSSIINEVKTvtfLDVTEKSKKsvcdIMTELKTTNYEYFDQPIEQYHGEKFDLAIASLSFHHIDDYLS 119 3445899***********99999999999333347999987444456899*******99877777899***9888777888899999 PP DREV 170 lledirsalepngrviv 186 +++ + s l+ n++viv tr|A2FAE3|A2FAE3_TRIVA 120 VIKILPSHLNENANVIV 136 99999999999999998 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (265 nodes) Target sequences: 50190 (17055200 residues searched) Passed MSV filter: 1166 (0.0232317); expected 1003.8 (0.02) Passed bias filter: 478 (0.00952381); expected 1003.8 (0.02) Passed Vit filter: 33 (0.000657501); expected 50.2 (0.001) Passed Fwd filter: 4 (7.96972e-05); expected 0.5 (1e-05) Initial search space (Z): 50190 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.32u 0.02s 00:00:00.34 Elapsed: 00:00:00.15 # Mc/sec: 30130.85 // [ok]