# 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/TRM13.hmm.txt # target sequence database: proteomes/Homo_sapiens.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: TRM13 [M=262] Accession: PF05206.13 Description: Methyltransferase TRM13 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.6e-91 308.7 0.2 2.7e-91 308.0 0.2 1.3 1 sp|Q9NUP7|TRM13_HUMAN tRNA:m(4)X modification enzyme TRM13 0.0001 24.2 1.3 0.68 11.7 1.1 2.5 2 sp|Q8NEC7-3|GSTCD_HUMAN Isoform 2 of Glutathione S-transfera 0.0025 19.7 5.9 0.84 11.4 1.1 2.8 3 sp|Q8NEC7|GSTCD_HUMAN Glutathione S-transferase C-terminal ------ inclusion threshold ------ 0.15 13.8 0.0 21 6.8 0.0 2.2 2 sp|Q96FB5|RRNAD_HUMAN Protein RRNAD1 OS=Homo sapiens GN=RR Domain annotation for each sequence (and alignments): >> sp|Q9NUP7|TRM13_HUMAN tRNA:m(4)X modification enzyme TRM13 homolog OS=Homo sapiens GN=TRMT13 PE=1 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 308.0 0.2 1.2e-95 2.7e-91 1 261 [. 165 469 .. 165 470 .. 0.96 Alignments for each domain: == domain 1 score: 308.0 bits; conditional E-value: 1.2e-95 TRM13 1 lqqsslignleelkllekssayvEfGaGkgelsryvnqaleedkkenskyvLidrasnrlkaDrklrkdeselaikRlkidIkdLnls 88 +qq+s++gn+e+lkll +++++vEfGaGkg+ls++v+ al+++ e+++++L+++ ++r+k+D k+rk++s +++Rl+idI++L+l+ sp|Q9NUP7|TRM13_HUMAN 165 KQQASILGNIENLKLLGPRRCFVEFGAGKGKLSHWVDIALKDA--EKVHFILVEKVTTRFKVDGKHRKKNS--VFERLQIDIQHLCLN 248 69**************************************777..9************************6..*************** PP TRM13 89 aleele.kkkkvvavsKHLCGaATDLtLrcllnsekask....................................kakleglliAlCC 139 ++++l+ +k +vv+++KHLCG+ATDL+Lrcl+++++as ++ ++g++iAlCC sp|Q9NUP7|TRM13_HUMAN 249 KIPVLReEKLPVVGIGKHLCGMATDLALRCLVETYAASFeerneeplakrikndktekeiytlakegneknvpekWNPVAGIVIALCC 336 ******999***************************9999************************************************ PP TRM13 140 hhvCswkeyvnkeyleelgitkeefqiltkmvSWavsgkrkse...........aeeedekeekekeekkeeeeseeekksklsseer 216 hh+C+w++yv+key+++lg+ + ef+++++m+SWa++g+rk++ ++++d +ee+++ + ++++ + + ls+ee+ sp|Q9NUP7|TRM13_HUMAN 337 HHRCDWRHYVGKEYFRALGLGAVEFHYFQRMSSWATCGMRKTSletsnsttkrqDNQNDDSEEHDDGGYRITDDGADCLPGLLSVEEK 424 ****************************************8889999999987656666677777777788999999999******** PP TRM13 217 eeiGlkakrlidegRllalkekgfeaelvkYvekevslEnvlLla 261 ++iG+ +k lid+gR+++l++kgf+++l++Y+++ vslEnvlL+a sp|Q9NUP7|TRM13_HUMAN 425 KKIGHLCKLLIDQGRIQYLQQKGFSPALQYYTDPLVSLENVLLTA 469 *******************************************87 PP >> sp|Q8NEC7-3|GSTCD_HUMAN Isoform 2 of Glutathione S-transferase C-terminal domain-containing protein OS=Homo sapiens # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 11.7 1.1 3.1e-05 0.68 75 169 .. 391 485 .. 347 493 .. 0.71 2 ! 10.4 0.0 7.5e-05 1.7 215 260 .. 497 542 .. 487 544 .. 0.92 Alignments for each domain: == domain 1 score: 11.7 bits; conditional E-value: 3.1e-05 TRM13 75 ikRlkidIkdLnlsaleelekkkkv......vavsKHLCGaATDLtLrcllnsekaskkakleglliAlCChhvCswkeyvn..ke 152 + R k +L ls++ ++ + ++ ++v+ H CG+ATD+ ++ +++++++ ++ CC + + n k+ sp|Q8NEC7-3|GSTCD_HUMAN 391 LIRAKKRSDELGLSNIWFIQANMEYftgmfnIGVALHACGVATDMVIEHCIKTRAS--------FVTCPCCYGFIQNTSKFNfpKS 468 5555555566666666666655555566666999*************999988774........5666677655555444444422 PP TRM13 153 yleelgitkeefqiltk 169 ++ ++++e il++ sp|Q8NEC7-3|GSTCD_HUMAN 469 EQFKKTLSYKEHMILCR 485 23345555555555555 PP == domain 2 score: 10.4 bits; conditional E-value: 7.5e-05 TRM13 215 ereeiGlkakrlidegRllalkekgfeaelvkYvekevslEnvlLl 260 +r+ iG+++ l+d R+ a +e g++++++ +++ s+ n +++ sp|Q8NEC7-3|GSTCD_HUMAN 497 QRRLIGKQCMCLVDLDRARAAEECGYSVQVISMEPESCSPKNNMIV 542 7999*************************************99886 PP >> sp|Q8NEC7|GSTCD_HUMAN Glutathione S-transferase C-terminal domain-containing protein OS=Homo sapiens GN=GSTCD PE=1 S # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -3.1 0.2 1 2.2e+04 127 174 .. 97 147 .. 87 156 .. 0.70 2 ! 11.4 1.1 3.8e-05 0.84 75 169 .. 478 572 .. 436 581 .. 0.71 3 ! 10.3 0.0 8.1e-05 1.8 215 260 .. 584 629 .. 563 631 .. 0.92 Alignments for each domain: == domain 1 score: -3.1 bits; conditional E-value: 1 TRM13 127 kakleglliAlCChhvCswkeyvn....keyleelgitkeefqiltkmvSWa 174 +++++ + +A+ +h+ + + y ke+le lg++k+ +++ +++ W sp|Q8NEC7|GSTCD_HUMAN 97 SDNFCRAGLAVVLRHIIQKS-YEAdplkKELLELLGFKKTCLKACAEVSQWT 147 46777777788888887654.4333666888888999999999988888885 PP == domain 2 score: 11.4 bits; conditional E-value: 3.8e-05 TRM13 75 ikRlkidIkdLnlsaleelekkkkv......vavsKHLCGaATDLtLrcllnsekaskkakleglliAlCChhvCswkeyvn..keyl 154 + R k +L ls++ ++ + ++ ++v+ H CG+ATD+ ++ +++++++ ++ CC + + n k+ sp|Q8NEC7|GSTCD_HUMAN 478 LIRAKKRSDELGLSNIWFIQANMEYftgmfnIGVALHACGVATDMVIEHCIKTRAS--------FVTCPCCYGFIQNTSKFNfpKSEQ 557 5555555566666666666645555566666999*************999988774........566667765555544444442223 PP TRM13 155 eelgitkeefqiltk 169 ++ ++++e il++ sp|Q8NEC7|GSTCD_HUMAN 558 FKKTLSYKEHMILCR 572 345555555555555 PP == domain 3 score: 10.3 bits; conditional E-value: 8.1e-05 TRM13 215 ereeiGlkakrlidegRllalkekgfeaelvkYvekevslEnvlLl 260 +r+ iG+++ l+d R+ a +e g++++++ +++ s+ n +++ sp|Q8NEC7|GSTCD_HUMAN 584 QRRLIGKQCMCLVDLDRARAAEECGYSVQVISMEPESCSPKNNMIV 629 7999*************************************99886 PP >> sp|Q96FB5|RRNAD_HUMAN Protein RRNAD1 OS=Homo sapiens GN=RRNAD1 PE=2 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 4.7 0.0 0.0041 90 9 55 .. 144 191 .. 138 214 .. 0.76 2 ? 6.8 0.0 0.00097 21 226 261 .. 423 458 .. 417 459 .. 0.90 Alignments for each domain: == domain 1 score: 4.7 bits; conditional E-value: 0.0041 TRM13 9 nleelkllekssayvEfGaGkgelsryvnqaleed.kkenskyvLidr 55 +++l +++ +++v+ G+G+g+lsr+++ l k+ + L++r sp|Q96FB5|RRNAD_HUMAN 144 LVKKLSDFTGCTQVVDVGSGQGHLSRFMALGLGLMvKSIEGDQRLVER 191 6778888899999**************987554332344445555655 PP == domain 2 score: 6.8 bits; conditional E-value: 0.00097 TRM13 226 lidegRllalkekgfeaelvkYvekevslEnvlLla 261 li Rll+l+e+gf+ael +e+s+ n++L+a sp|Q96FB5|RRNAD_HUMAN 423 LILLDRLLYLQEQGFHAELLPIFSPELSPRNLVLVA 458 66667*****************************98 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (262 nodes) Target sequences: 88473 (35070517 residues searched) Passed MSV filter: 5329 (0.0602331); expected 1769.5 (0.02) Passed bias filter: 2331 (0.026347); expected 1769.5 (0.02) Passed Vit filter: 209 (0.0023623); expected 88.5 (0.001) Passed Fwd filter: 6 (6.78173e-05); expected 0.9 (1e-05) Initial search space (Z): 88473 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.89u 0.03s 00:00:00.92 Elapsed: 00:00:00.27 # Mc/sec: 34031.39 // [ok]