# 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_8.hmm.txt # target sequence database: proteomes/Monosiga_brevicollis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Methyltransf_8 [M=219] Accession: PF05148.15 Description: Hypothetical methyltransferase 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.2e-74 248.1 0.1 3.9e-72 241.7 0.1 2.1 1 tr|A9VA61|A9VA61_MONBE Predicted protein OS=Monosiga brevico 0.00063 18.7 0.0 0.13 11.1 0.0 2.2 2 tr|A9UVH5|A9UVH5_MONBE Predicted protein OS=Monosiga brevico ------ inclusion threshold ------ 0.028 13.3 0.0 0.054 12.4 0.0 1.4 1 tr|A9UU62|A9UU62_MONBE Predicted protein OS=Monosiga brevico Domain annotation for each sequence (and alignments): >> tr|A9VA61|A9VA61_MONBE Predicted protein OS=Monosiga brevicollis OX=81824 GN=29115 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 241.7 0.1 1.3e-75 3.9e-72 7 219 .] 216 449 .] 212 449 .] 0.90 Alignments for each domain: == domain 1 score: 241.7 bits; conditional E-value: 1.3e-75 -HHHHHHHHHHHHHHHHHS-HHHHHHHHHH-HHHHHHHHHHHHHHHCTSSS-HHHHHHHHHCTS-TTS-EEEES-TT-HHHHH--S- CS Methyltransf_8 7 kkrldgarfrflneklytgesskavdlfkedpeafdlyhkgfrsqvkkWPenPldviirklkkkkdnlvvadlGcGeakiaeavrek 93 + rl+garfr+lneklyt++++ a++ fke pe+fd+yhkgf++qv++WP+nP+d++i+ + +k+++lvvad+GcGeak+ +v +k tr|A9VA61|A9VA61_MONBE 216 SARLQGARFRMLNEKLYTTTGDDAFRWFKESPELFDVYHKGFATQVQRWPVNPVDRMIEFVLQKPAKLVVADMGCGEAKLGASVPNK 302 679******************************************************************************777666 PP ..EEEEESS-SSTTEEES-TTS-S--TT-EEEEEEES---SS-HHHHHHHHHHHEEEEEEEEEEE.................EGGG- CS Methyltransf_8 94 ekvrsfdlvavnervvacdlakvPledesvdvavfclslmgtnladflkeanrvlkvGGllkvae.................vksrf 163 v+sfdlva n+ v+acd+a+vPl+de +d+ +fcl+lmgtn+ d+l ea r+lk+ G+lk+ae vksr+ tr|A9VA61|A9VA61_MONBE 303 --VHSFDLVAANPSVTACDIAHVPLADEKIDICIFCLALMGTNYVDYLLEAFRILKPRGILKIAElqgdakahqnawpsceqVKSRI 387 ..**************************************************************855555555666666666789** PP S-HHHHHHHHHCTTEEEEEEE--STTEEEEEEEE-S......SS-TT---GG--------XX CS Methyltransf_8 164 esvklfaraltklGfevkskdlsnkqfvlfefkkve......kkgkkrelglklkpclykkr 219 + f+r ++ +Gfe +s d n++fv +ef ++ + ++ + l+pc+yk+r tr|A9VA61|A9VA61_MONBE 388 GESSAFVRMVQDVGFELTSEDARNTHFVDYEFVRMRtpgrrqLIRNSQQTQQVLQPCIYKRR 449 9*******************************988733432233445566888********9 PP >> tr|A9UVH5|A9UVH5_MONBE Predicted protein OS=Monosiga brevicollis OX=81824 GN=6966 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 5.5 0.0 0.0023 7.1 40 85 .. 32 81 .. 20 92 .. 0.77 2 ! 11.1 0.0 4.3e-05 0.13 110 157 .. 130 177 .. 112 195 .. 0.85 Alignments for each domain: == domain 1 score: 5.5 bits; conditional E-value: 0.0023 HHHHHHHHHHHHHCTSSS-HHHHHHHHHCTS-T....TS-EEEES-TT-H CS Methyltransf_8 40 afdlyhkgfrsqvkkWPenPldviirklkkkkd....nlvvadlGcGeak 85 +f y + +++ ++P d i++++ + d + + ad+GcG + tr|A9UVH5|A9UVH5_MONBE 32 VFGQYADEYDQLRPHYPASLWDLAIKHVQMRTDtepcSWTAADIGCGTGR 81 577788888888889999999********665511116789******877 PP == domain 2 score: 11.1 bits; conditional E-value: 4.3e-05 ES-TTS-S--TT-EEEEEEES---SS-HHHHHHHHHHHEEEEEEEEEE CS Methyltransf_8 110 acdlakvPledesvdvavfclslmgtnladflkeanrvlkvGGllkva 157 c + l +s d++ + + +l l+e++r+lk+GG+l +a tr|A9UVH5|A9UVH5_MONBE 130 QCTAEQTGLPAHSCDLVTCLQAWHWVDLVPGLEEMRRILKPGGVLALA 177 677788888889999977666778899999***************987 PP >> tr|A9UU62|A9UU62_MONBE Predicted protein OS=Monosiga brevicollis OX=81824 GN=36263 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.4 0.0 1.8e-05 0.054 67 163 .. 579 677 .. 569 695 .. 0.67 Alignments for each domain: == domain 1 score: 12.4 bits; conditional E-value: 1.8e-05 HCTS-TTS-EEEES-TT-HHHHH--S-..EE.EEESS-SSTTEEES-T..TS-S--TT-EEEEEEES--.-SS-HHHHHHHHHHHEE CS Methyltransf_8 67 lkkkkdnlvvadlGcGeakiaeavrekekvr.sfdlvavnervvacdl..akvPledesvdvavfclsl.mgtnladflkeanrvlk 149 lk++ +++ + d+ +a + e + ++lv + d + +P+ +svd++ f sl m + ++ l++a r+l tr|A9UU62|A9UU62_MONBE 579 LKRRFKKIYILDISASLLDVARRRVEAAGLSdVVELVLAD-FTDLSDKgkSALPA-PGSVDLVTFSYSLsMIPDKRSALQQATRLLA 663 6666666666666666666663333332222022222222.22222322144554.589*******9998899************** PP EEEEEEEEEEGGG- CS Methyltransf_8 150 vGGllkvaevksrf 163 +GG l +a+ r tr|A9UU62|A9UU62_MONBE 664 PGGSLGIADFFYRN 677 ********987765 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (219 nodes) Target sequences: 9188 (5511963 residues searched) Passed MSV filter: 165 (0.0179582); expected 183.8 (0.02) Passed bias filter: 147 (0.0159991); expected 183.8 (0.02) Passed Vit filter: 12 (0.00130605); expected 9.2 (0.001) Passed Fwd filter: 3 (0.000326513); expected 0.1 (1e-05) Initial search space (Z): 9188 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.09u 0.00s 00:00:00.09 Elapsed: 00:00:00.04 # Mc/sec: 30178.00 // [ok]