# 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/Branchiostoma_floridae.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.9e-87 293.8 1.6 2.8e-87 293.2 1.6 1.3 1 tr|C3YCH1|C3YCH1_BRAFL Uncharacterized protein OS=Branchiost 0.0027 18.0 0.3 3.4 7.8 0.0 2.4 3 tr|C3Z7E1|C3Z7E1_BRAFL Uncharacterized protein OS=Branchiost ------ inclusion threshold ------ 0.14 12.3 0.0 0.18 11.9 0.0 1.2 1 tr|C3ZHD7|C3ZHD7_BRAFL Uncharacterized protein OS=Branchiost 0.23 11.6 0.0 0.34 11.0 0.0 1.3 1 tr|C3Y0I3|C3Y0I3_BRAFL Uncharacterized protein OS=Branchiost Domain annotation for each sequence (and alignments): >> tr|C3YCH1|C3YCH1_BRAFL Uncharacterized protein OS=Branchiostoma floridae OX=7739 GN=BRAFLDRAFT_285163 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 293.2 1.6 3.9e-91 2.8e-87 1 261 [. 160 455 .. 160 456 .. 0.94 Alignments for each domain: == domain 1 score: 293.2 bits; conditional E-value: 3.9e-91 TRM13 1 lqqsslignleelkllekssayvEfGaGkgelsryvnqaleedkkenskyvLidrasnrlkaDrklrkdeselaikRlkidIkdLnl 87 +qq+slig+l++l+ll+++ +++EfGaG+g+ls+yv+++l e+ ++s+++ idra +++k+D+ +r+ e+ + ++Rlk+dI+dLnl tr|C3YCH1|C3YCH1_BRAFL 160 KQQASLIGHLDHLDLLKSDVCFIEFGAGRGRLSHYVQMCLPEQ--AASHFIVIDRANTKHKFDKNIRRGEKGPDFERLKVDIEDLNL 244 69**************************************999..99**************************************** PP TRM13 88 saleele.kkkkvvavsKHLCGaATDLtLrcllnsekask.....................................kakleglliA 136 ++++ +k +v a+sKHLCGaATDLtLrcl+++ ++++ + e +liA tr|C3YCH1|C3YCH1_BRAFL 245 GLVPSMAeHKGPVAAISKHLCGAATDLTLRCLIETLHCKSsrnsskevftssstsetasrqhpvsdniaqesnttekQLVDECILIA 331 *****998889999**********************998899********************************9744445789*** PP TRM13 137 lCChhvCswkeyvnkeyleelgitkeefqiltkmvSWavsgkrkseaeeedekeekekeekkeeeeseeekksklsseereeiGlka 223 lCChh+C w +yv+k ++++ g+t+e+f++l +m+SWa++g + ++++ +e+e +e++eee++ee+ ++eereeiG+++ tr|C3YCH1|C3YCH1_BRAFL 332 LCCHHRCTWTSYVGKGFFTSIGFTEEDFHLLRSMSSWATCGFARGQTKNAG-GNEEEVKENSEEEGEEEKAAAIWTVEEREEIGRQC 417 *****************************************9977777766.44555556667888888999999************ PP TRM13 224 krlidegRllalkekgfeaelvkYvekevslEnvlLla 261 krl+d+gRl++l++ +f+++lv Yve++vslEn++++a tr|C3YCH1|C3YCH1_BRAFL 418 KRLLDQGRLWYLQQCDFNTQLVHYVEPDVSLENIAIVA 455 ***********************************997 PP >> tr|C3Z7E1|C3Z7E1_BRAFL Uncharacterized protein OS=Branchiostoma floridae OX=7739 GN=BRAFLDRAFT_200329 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 7.6 0.2 0.00054 3.9 20 109 .. 156 263 .. 134 288 .. 0.74 2 ? -2.9 0.0 0.86 6.1e+03 31 64 .. 300 332 .. 269 339 .. 0.73 3 ! 7.8 0.0 0.00048 3.4 229 261 .. 428 460 .. 415 461 .. 0.89 Alignments for each domain: == domain 1 score: 7.6 bits; conditional E-value: 0.00054 TRM13 20 sayvEfGaGkgelsryvnqaleed.............kkenskyvLidrasnrlkaDrklrkdese.....laikRlkidIkdLnls 88 s++v++G+G+g+lsr++ a++ + + k v d +++++k+ + de +++ lk d++++ tr|C3Z7E1|C3Z7E1_BRAFL 156 SNVVDIGSGQGHLSRVLSFAHRLNvtsveaegchvtgAARFDKQVVTDITKEKMKKQQEADIDEVLhaqpsHVTCFLKP---DISVD 239 5799***************999889999999999998888888899999999999877777666434333333333333...55566 PP TRM13 89 aleele...kkkkvvavsKHLCGa 109 +++++ ++k + +v+ H CG tr|C3Z7E1|C3Z7E1_BRAFL 240 DFLSVDapcDRKAFSMVGLHTCGD 263 666666666888889999999995 PP == domain 2 score: -2.9 bits; conditional E-value: 0.86 TRM13 31 elsryvnqaleedkkenskyvLidrasnrlkaDr 64 ++++ +n ++ + + + + id ++rlk D tr|C3Z7E1|C3Z7E1_BRAFL 300 QMAQNLNAEMSYQ-ARELACHAIDMYTERLKSDD 332 5555555555444.44455556888888888886 PP == domain 3 score: 7.8 bits; conditional E-value: 0.00048 TRM13 229 egRllalkekgfeaelvkYvekevslEnvlLla 261 R+++lke+g + +lv ++ +s+ n++Lla tr|C3Z7E1|C3Z7E1_BRAFL 428 LDRMQYLKEHGTNCQLVPLFNPRLSPRNFVLLA 460 56*****************************98 PP >> tr|C3ZHD7|C3ZHD7_BRAFL Uncharacterized protein OS=Branchiostoma floridae OX=7739 GN=BRAFLDRAFT_99360 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.9 0.0 2.6e-05 0.18 2 41 .. 55 94 .. 54 126 .. 0.83 Alignments for each domain: == domain 1 score: 11.9 bits; conditional E-value: 2.6e-05 TRM13 2 qqsslignleelkllekssayvEfGaGkgelsryvnqale 41 qq+ + ++++ +++ ++++ G+G ge+s+yv+q + tr|C3ZHD7|C3ZHD7_BRAFL 55 QQEIAVMLMHQNFEWQEGDTVLDVGCGTGEISKYVAQQHG 94 666667778888889*********************9874 PP >> tr|C3Y0I3|C3Y0I3_BRAFL Uncharacterized protein OS=Branchiostoma floridae OX=7739 GN=BRAFLDRAFT_92966 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.0 0.0 4.8e-05 0.34 21 44 .. 49 73 .. 23 94 .. 0.78 Alignments for each domain: == domain 1 score: 11.0 bits; conditional E-value: 4.8e-05 TRM13 21 ayvEfGaGkgelsryvnqaleed.k 44 + +++G+Gkg++s+y+n+ + e tr|C3Y0I3|C3Y0I3_BRAFL 49 VALDLGCGKGYISKYLNKDIVETlY 73 6789**************9987752 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (262 nodes) Target sequences: 28542 (13248087 residues searched) Passed MSV filter: 1724 (0.0604022); expected 570.8 (0.02) Passed bias filter: 769 (0.0269428); expected 570.8 (0.02) Passed Vit filter: 72 (0.0025226); expected 28.5 (0.001) Passed Fwd filter: 6 (0.000210217); expected 0.3 (1e-05) Initial search space (Z): 28542 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.34u 0.01s 00:00:00.35 Elapsed: 00:00:00.10 # Mc/sec: 34709.99 // [ok]