# 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/GCD14.hmm.txt # target sequence database: proteomes/Trichomonas_vaginalis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: GCD14 [M=247] Accession: PF08704.9 Description: tRNA methyltransferase complex GCD14 subunit 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.8e-06 27.9 0.0 7.2e-06 27.3 0.0 1.2 1 tr|A2E9S9|A2E9S9_TRIVA Fibrillarin, putative OS=Trichomonas 2.2e-05 25.8 0.0 4.2e-05 24.8 0.0 1.3 1 tr|A2DXE8|A2DXE8_TRIVA Fibrillarin, putative OS=Trichomonas ------ inclusion threshold ------ 1.5 9.9 5.1 1.8 9.6 5.1 1.2 1 tr|A2DP68|A2DP68_TRIVA Uncharacterized protein OS=Trichomona 3.3 8.8 4.5 48 5.0 0.0 2.5 2 tr|A2EJG2|A2EJG2_TRIVA Uncharacterized protein OS=Trichomona Domain annotation for each sequence (and alignments): >> tr|A2E9S9|A2E9S9_TRIVA Fibrillarin, putative OS=Trichomonas vaginalis OX=5722 GN=TVAG_282140 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 27.3 0.0 5.8e-10 7.2e-06 37 164 .. 157 292 .. 146 307 .. 0.76 Alignments for each domain: == domain 1 score: 27.3 bits; conditional E-value: 5.8e-10 GCD14 37 elkpGsvvvesGtGsgslshaiartvaPtGhlftfefheqraekareefrehelsqlvtvthrDvc.keGfdeevsgkvDavfLDlP 122 + pG++v+ G+ sg+ + ++ v+PtG+++++ef + r ar+ +l ++ D+ + + + + vD +f D+ tr|A2E9S9|A2E9S9_TRIVA 157 GICPGAKVLYLGAASGTTVSHVSDIVGPTGRVYAVEFSH-R--VARDLVSLAKLRPNIVPIINDARyPNKY-RMLVPMVDYMFADVA 239 6889**************9******************75.3..5777888889988777777888735555.66779********** PP GCD14 123 aPwea...iphaakalkkegg.rlCsfsPCieqv.......qrtceaLrelgf 164 +P++a +a + lk+ gg +lC + C++ q++ e+L++ gf tr|A2E9S9|A2E9S9_TRIVA 240 QPDQArifAINANNFLKSGGGfMLCVKASCVNSSedpetiyQQAQEELTKEGF 292 ****9432245667788777657899999*97542222223566666666666 PP >> tr|A2DXE8|A2DXE8_TRIVA Fibrillarin, putative OS=Trichomonas vaginalis OX=5722 GN=TVAG_380100 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 24.8 0.0 3.3e-09 4.2e-05 35 168 .. 143 284 .. 138 293 .. 0.70 Alignments for each domain: == domain 1 score: 24.8 bits; conditional E-value: 3.3e-09 GCD14 35 kLelkpGsvvvesGtGsgslshaiartvaPtGhlftfefheqraekareefrehelsqlvtvthrDvc.keGfdeevsgkvDavfLD 120 k+ +kpG+ v+ G+ sg+ + ++ v+PtG+++++ef +r+ + + ++ + + ++v + D+ + + + + vD +f D tr|A2DXE8|A2DXE8_TRIVA 143 KFGIKPGASVLYLGAASGTTVSHVSDIVGPTGRVYAVEFS-HRVGRELCNLAKSRP-NIVPI-IEDARhPQKY-RMLVPMVDYMFAD 225 6779****************9******************6.577777777766553.33333.3455514445.556799******* PP GCD14 121 lPaPwea...iphaakalkkegg.rlCsfsPCieqv.......qrtceaLrelgfteie 168 + +P++a + + lk++gg +lC + C++ + ++L++ gf++ e tr|A2DXE8|A2DXE8_TRIVA 226 VAQPDQArifVLNGQFFLKNNGGyMLCVKASCVDSTldphkvyAMAQQELKDAGFQNEE 284 ******96444555556788776478999999986411111112234556666666555 PP >> tr|A2DP68|A2DP68_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_015980 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 9.6 5.1 0.00014 1.8 153 232 .. 14 95 .. 6 143 .. 0.77 Alignments for each domain: == domain 1 score: 9.6 bits; conditional E-value: 0.00014 GCD14 153 qrtceaLrelgfteietlevllrtyevrks..kleeidlaekrleenerkriekakrrkaskdakvedkssks.stevllkPv 232 r +e+L e + +ei tl +++ + r +++++ +++kr e+ e+++ ++a ++k+sk++k++ kss+ t v kP tr|A2DP68|A2DP68_TRIVA 14 ERVVEKLIE-KNNEIATLKDEIKELKNRLLpyEMNTMTETQKRAEQIENDKYKQAIKQKTSKETKTSSKSSQNtATPVPEKPS 95 678889988.779*****99999999998711455556678888888888888888888888888777776642344444443 PP >> tr|A2EJG2|A2EJG2_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_050340 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 1.8 3.8 0.035 4.4e+02 175 234 .. 36 97 .. 21 103 .. 0.63 2 ? 5.0 0.0 0.0038 48 51 91 .. 258 298 .. 237 302 .. 0.83 Alignments for each domain: == domain 1 score: 1.8 bits; conditional E-value: 0.035 GCD14 175 rtyevrkskleeidlaekrleenerkriekakrrkaskdakvedkssks..stevllkPvse 234 ++v ++ + ++ + k+ ++n+++r+ek+ ++++++++++++s k ++ kPv++ tr|A2EJG2|A2EJG2_TRIVA 36 SDFDVPEMYFYFTSNGSKANKKNKKTRAEKRLAKSQKEQEAKSEESPKEeiPAKEEEKPVED 97 56777777777777788888888888887776666544444444333322245555666655 PP == domain 2 score: 5.0 bits; conditional E-value: 0.0038 GCD14 51 sgslshaiartvaPtGhlftfefheqraekareefrehels 91 g s ++ r + P G + fe ++ ++ +eef+ + ++ tr|A2EJG2|A2EJG2_TRIVA 258 IGDKSMSLKRVMYPDGEVEEFESKKNLTDADLEEFKANWIK 298 56668889****************************98776 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (247 nodes) Target sequences: 50190 (17055200 residues searched) Passed MSV filter: 4050 (0.0806934); expected 1003.8 (0.02) Passed bias filter: 1839 (0.0366408); expected 1003.8 (0.02) Passed Vit filter: 236 (0.00470213); expected 50.2 (0.001) Passed Fwd filter: 47 (0.000936442); 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: 1.13u 0.03s 00:00:01.15 Elapsed: 00:00:00.32 # Mc/sec: 13164.48 // [ok]