# 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/Mus_musculus.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 6.9e-94 314.5 0.0 9.3e-94 314.1 0.0 1.2 1 sp|Q8BYH3|TRM13_MOUSE tRNA:m(4)X modification enzyme TRM13 h 1.6e-05 24.9 0.4 0.27 11.1 0.9 2.7 2 sp|Q5RL51|GSTCD_MOUSE Glutathione S-transferase C-terminal d ------ inclusion threshold ------ 0.14 11.9 0.0 15 5.3 0.0 2.3 2 sp|Q8BZG5|RRNAD_MOUSE Protein RRNAD1 OS=Mus musculus OX=1009 Domain annotation for each sequence (and alignments): >> sp|Q8BYH3|TRM13_MOUSE tRNA:m(4)X modification enzyme TRM13 homolog OS=Mus musculus OX=10090 GN=Trmt13 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 314.1 0.0 1.3e-97 9.3e-94 1 261 [. 165 468 .. 165 469 .. 0.98 Alignments for each domain: == domain 1 score: 314.1 bits; conditional E-value: 1.3e-97 TRM13 1 lqqsslignleelkllekssayvEfGaGkgelsryvnqaleedkkenskyvLidrasnrlkaDrklrkdeselaikRlkidIkdLnls 88 +qq+s++gn+e+lkll +++++vEfGaGkg+ls++v+ al+++ en++++L++r ++r+k+D k+rk++s +++Rl+idI++L+l+ sp|Q8BYH3|TRM13_MOUSE 165 KQQASILGNIEKLKLLGPRRCFVEFGAGKGKLSHWVDIALKDA--ENVHFILVERVTTRFKVDGKHRKKDS--VFERLQIDIQHLCLN 248 69**************************************777..*************************7..*************** PP TRM13 89 aleele.kkkkvvavsKHLCGaATDLtLrcllnsekask....................................kakleglliAlCC 139 ++++l+ + +vv+++KHLCG+ATD++Lrcl+++++as ++ ++g++iAlCC sp|Q8BYH3|TRM13_MOUSE 249 RVPVLReGRLPVVGIGKHLCGVATDVALRCLVETYAASFeekdeeplakrikndktekesntlakegsekdvpetWTPVAGIVIALCC 336 ******999***************************99999*********************************************** PP TRM13 140 hhvCswkeyvnkeyleelgitkeefqiltkmvSWavsgkrkse..........aeeedekeekekeekkeeeeseeekksklsseere 217 hh+C+w++yv+key+++lg+ + ef+++++m+SWa++g+r+s ++++de+ee+++ ++ +++++++ + l++ee++ sp|Q8BYH3|TRM13_MOUSE 337 HHRCDWRHYVGKEYFKALGLGAVEFYYFQRMSSWATCGMRTSLegsdvtperkDAQRDENEEHDDGGDRLTDGNTDSLPGILTVEEKK 424 ****************************************8889*******998889999999999********************** PP TRM13 218 eiGlkakrlidegRllalkekgfeaelvkYvekevslEnvlLla 261 +iG+ +k lid+gRl++l++kgf+++l++Y+++ vslEnvlL+a sp|Q8BYH3|TRM13_MOUSE 425 KIGHLCKLLIDQGRLQYLQQKGFSPALQYYTDPLVSLENVLLTA 468 ******************************************97 PP >> sp|Q5RL51|GSTCD_MOUSE Glutathione S-transferase C-terminal domain-containing protein OS=Mus musculus OX=10090 GN=Gst # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 11.1 0.9 3.6e-05 0.27 75 140 .. 479 542 .. 438 581 .. 0.72 2 ! 11.0 0.0 3.7e-05 0.28 214 260 .. 584 630 .. 555 632 .. 0.89 Alignments for each domain: == domain 1 score: 11.1 bits; conditional E-value: 3.6e-05 TRM13 75 ikRlkidIkdLnlsaleelekkkkv......vavsKHLCGaATDLtLrcllnsekaskkakleglliAlCCh 140 + R k +L ls++ ++ + ++ ++v+ H CG+ATD+ ++ +++++++ ++ CC sp|Q5RL51|GSTCD_MOUSE 479 LIRAKKRSDELGLSNIWFIQANMEYftgmfnIGVALHACGVATDMVIEHCIQTRAS--------FITCPCCY 542 5555555566666666666645555566666999*************999988764........34445555 PP == domain 2 score: 11.0 bits; conditional E-value: 3.7e-05 TRM13 214 eereeiGlkakrlidegRllalkekgfeaelvkYvekevslEnvlLl 260 er+ iG+++ l+d R+ a e+g++++++ +++ s+ n +++ sp|Q5RL51|GSTCD_MOUSE 584 PERRLIGKQCMGLVDLDRAAAAGEHGYSVQVISMEPESCSPKNNMIV 630 47999*************************************99886 PP >> sp|Q8BZG5|RRNAD_MOUSE Protein RRNAD1 OS=Mus musculus OX=10090 GN=Rrnad1 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 4.2 0.0 0.0043 32 8 38 .. 143 173 .. 138 215 .. 0.81 2 ? 5.3 0.0 0.002 15 226 261 .. 423 458 .. 417 459 .. 0.90 Alignments for each domain: == domain 1 score: 4.2 bits; conditional E-value: 0.0043 TRM13 8 gnleelkllekssayvEfGaGkgelsryvnq 38 + +++l +l+ +++v+ G+G+g+lsr++ sp|Q8BZG5|RRNAD_MOUSE 143 ELVKKLSDLTGCTQVVDVGSGQGHLSRFMSL 173 5678888999999**************9876 PP == domain 2 score: 5.3 bits; conditional E-value: 0.002 TRM13 226 lidegRllalkekgfeaelvkYvekevslEnvlLla 261 li R+l+l+e+gf ael +e+s+ n++L+a sp|Q8BZG5|RRNAD_MOUSE 423 LILLDRMLYLQEQGFYAELLPIFSPELSPRNLVLVA 458 66667*****************************98 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (262 nodes) Target sequences: 22289 (11757923 residues searched) Passed MSV filter: 1529 (0.0685989); expected 445.8 (0.02) Passed bias filter: 648 (0.0290726); expected 445.8 (0.02) Passed Vit filter: 58 (0.00260218); expected 22.3 (0.001) Passed Fwd filter: 5 (0.000224326); expected 0.2 (1e-05) Initial search space (Z): 22289 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.31u 0.01s 00:00:00.32 Elapsed: 00:00:00.10 # Mc/sec: 30805.76 // [ok]