# 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/MT-A70.hmm.txt # target sequence database: proteomes/Petromyzon_marinus.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: MT-A70 [M=172] Accession: PF05063.13 Description: MT-A70 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5.1e-59 198.6 0.0 7.8e-59 198.0 0.0 1.3 1 tr|S4RU24|S4RU24_PETMA Methyltransferase like 3 OS=Petromyzo 1.4e-24 86.6 0.1 1.9e-24 86.1 0.1 1.2 1 tr|S4RYP4|S4RYP4_PETMA Methyltransferase like 4 OS=Petromyzo Domain annotation for each sequence (and alignments): >> tr|S4RU24|S4RU24_PETMA Methyltransferase like 3 OS=Petromyzon marinus OX=7757 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 198.0 0.0 1.2e-62 7.8e-59 1 172 [] 410 571 .. 410 571 .. 0.98 Alignments for each domain: == domain 1 score: 198.0 bits; conditional E-value: 1.2e-62 ECEEEE---BHHCC..HTT--SS--B-HHHHHTS-GGGCECCEEEEEEEEETTCHHHHHHHHHHHCT-EEEEEEEEEEE-TTSCBS- CS MT-A70 1 fdviiidPpweeklvarkkelpyktlneeelkaLpipelldkrslvflWvtnrellelkrellkkWglelveeltwlkinklgepia 87 f+v+++dPpw++++ elpy+tl+++e+++L+ip+l+d++ +flWvt+r+++ + re+lk Wg+e+v+el+w+k+n+l+++i+ tr|S4RU24|S4RU24_PETMA 410 FSVVMADPPWDIHM-----ELPYGTLTDDEMRRLNIPVLQDDG-YIFLWVTGRAMELG-RECLKLWGYERVDELIWVKTNQLQRIIR 489 89************.....9*********************99.*************9.**************************** PP BS-TTSSSBB-EEEEEEEEES---TTT.TTB-TTSSBSEEEEE-CTTTSS-THHHHHHHHHCTCS-EEEET--G..--TTEEEEETT CS MT-A70 88 tlesgavlkkskEhllvgvkgnvkrstdkiinknvdtdliisperehsrKPpelfellehllesarkLelFaRs..lrpgwltvGne 172 t+++g++l++ kEh+lvg+kgn++ +n++ d+d++++++r++s+KP+e++ ++e+l++++rk+elF+R ++p+w+t+Gn+ tr|S4RU24|S4RU24_PETMA 490 TGRTGHWLNHGKEHCLVGMKGNTE-----GFNRGLDCDVLVAEVRSTSHKPDEIYGIIERLSPGTRKIELFGRPhnIQPNWVTLGNQ 571 ***********************8.....67899***************************************999*********97 PP >> tr|S4RYP4|S4RYP4_PETMA Methyltransferase like 4 OS=Petromyzon marinus OX=7757 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 86.1 0.1 2.8e-28 1.9e-24 1 139 [. 147 294 .. 147 295 .] 0.86 Alignments for each domain: == domain 1 score: 86.1 bits; conditional E-value: 2.8e-28 ECEEEE---BHHCC..HTT--SS--B-HHHHHTS-GGGCECCEEEEEEEEETTCHHHH..HHHHHHHCT-EEEEEEEEEEE-TTSCB CS MT-A70 1 fdviiidPpweeklvarkkelpyktlneeelkaLpipelldkrslvflWvtnrellel..krellkkWglelveeltwlkinklgep 85 fd+i++dPpwe+k+v+r+k +y l e ++ +Lp+p l+++++lv Wvtnr+ l + krel ++W + e+ w+k++ +ge tr|S4RYP4|S4RYP4_PETMA 147 FDLIVLDPPWENKSVKRSK--RYRALPEWQIARLPVPALASPGCLVLTWVTNRQRLLQfvKRELYPSWRVSDAVEWHWVKVTSHGEF 231 89****************9..9********************************999999*************************** PP S-BS-TTSSSBB-EEEEEEEEES.--.-TTT.TTB-TTS.........SBSEEEE.E-CTTTSS-T CS MT-A70 86 iatlesgavlkkskEhllvgvkg.nv.krstdkiinknv.........dtdliis.perehsrKPp 139 +++l+s +kk++E l++g ++ ++ +++ +++n +l++s p++ hs+KPp tr|S4RYP4|S4RYP4_PETMA 232 VFPLDSP--HKKPYETLVLGRVAaTPdMERYSQ-THRNPceitalklpTRKLLVSvPSILHSQKPP 294 **98888..99*******984331331345555.555544566777866677777999*******8 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (172 nodes) Target sequences: 13539 (5048892 residues searched) Passed MSV filter: 466 (0.0344191); expected 270.8 (0.02) Passed bias filter: 364 (0.0268853); expected 270.8 (0.02) Passed Vit filter: 10 (0.000738607); 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.08u 0.00s 00:00:00.08 Elapsed: 00:00:00.04 # Mc/sec: 21710.24 // [ok]