# 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/Methyltransf_10.hmm.txt # target sequence database: proteomes/Petromyzon_marinus.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Methyltransf_10 [M=299] Accession: PF05971.11 Description: Protein of unknown function (DUF890) 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.3e-83 280.3 0.0 1.5e-83 280.2 0.0 1.0 1 tr|S4RS03|S4RS03_PETMA Methyltransferase like 16 OS=Petromyz 9.7e-15 54.3 0.0 1.4e-14 53.8 0.0 1.2 1 tr|S4RH34|S4RH34_PETMA Uncharacterized protein OS=Petromyzon Domain annotation for each sequence (and alignments): >> tr|S4RS03|S4RS03_PETMA Methyltransferase like 16 OS=Petromyzon marinus OX=7757 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 280.2 0.0 2.2e-87 1.5e-83 27 255 .. 2 222 .] 1 222 [] 0.97 Alignments for each domain: == domain 1 score: 280.2 bits; conditional E-value: 2.2e-87 XXXXXXXXXXXXXXXXS-TTSHHHHHHHHHHHHHHHH----.---TTS----HHHHHHHHHHHHHHHT--TCGCS-.--EEEEES-T CS Methyltransf_10 27 PelkqfvilnpkGrvsidfadPlavkalnkallrefygvsewdiPdgfliPpvPgradyihyvadllgkadsdisa.qvraldiGvG 112 Pe+k++v+ +++Grv ++f++P+av+al++all+e++g+s+ +iP ++liP++P+r++yih+v+dllg++ s ++ +v+++diG+G tr|S4RS03|S4RS03_PETMA 2 PEFKKHVQTDLSGRVRLNFKEPEAVRALTCALLQEDFGLSI-EIPLERLIPTLPLRLNYIHWVEDLLGSSLSRDPEkSVHGIDIGTG 87 9****************************************.****************************999877699******** PP TTTHHHHHHHHHH--EEEEEES-HHHHHHHHHHHHHT.T-TTTEEEEE--ST-SSTTTSTT--S.-EEEEEE-----SS-------- CS Methyltransf_10 113 anCiyPllGvteyGWrfvgsevddvslasakaivesnsklsdaiklrkqkekkaifdglieene.rydvtlCnPPfhaslaeakags 198 a+CiyPllG+t +GW+f+++evd+v++++a+++v+ n l++++++ k+++k+ ++++l e e yd+++CnPPf+a+++ea +++ tr|S4RS03|S4RS03_PETMA 88 ASCIYPLLGATMNGWHFLATEVDNVCFTYATQNVARN-GLQGLVRVMKVPQKTLLVEALNEDPElAYDFCVCNPPFFANQQEAYGVN 173 *************************************.**************************9********************** PP ----...---.....---------TTTTHHHHTHHHHHHHHHHHHHHHGGGEEEEEE CS Methyltransf_10 199 erklknlgkanktralpslnfgGqkaelyceGGevafikklieeskqlakqvlWftt 255 +r+ ++ +ps++++G+ +e+++eGGe++f++k+i++s ql+k++++ ++ tr|S4RS03|S4RS03_PETMA 174 SRNS---HRP-----PPSSVNTGGATEIMAEGGELEFVRKIIHDSLQLQKRIRYTSC 222 **99...565.....899**********************************98765 PP >> tr|S4RH34|S4RH34_PETMA Uncharacterized protein OS=Petromyzon marinus OX=7757 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 53.8 0.0 2e-18 1.4e-14 251 298 .. 1 48 [. 1 49 [. 0.98 Alignments for each domain: == domain 1 score: 53.8 bits; conditional E-value: 2e-18 EEEEEEESSGGGHHHHHHHHHHTT-SEEEEEEEEETTEEEEEEEEES- CS Methyltransf_10 251 lWfttlvskkenlaplkeelkklGaskvtvvelaqGqklsrfiaWsfl 298 +W+t++++kk++laplkeel++ G+++ t++ ++qG+++++++aWsf tr|S4RH34|S4RH34_PETMA 1 RWYTCMLGKKCSLAPLKEELRRCGVPRTTTTTFCQGRTMRWALAWSFH 48 6*********************************************96 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (299 nodes) Target sequences: 13539 (5048892 residues searched) Passed MSV filter: 253 (0.0186868); expected 270.8 (0.02) Passed bias filter: 217 (0.0160278); expected 270.8 (0.02) Passed Vit filter: 9 (0.000664746); expected 13.5 (0.001) Passed Fwd filter: 2 (0.000147721); expected 0.1 (1e-05) Initial search space (Z): 13539 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.09u 0.00s 00:00:00.09 Elapsed: 00:00:00.05 # Mc/sec: 30192.37 // [ok]