# 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/Macaca_mulatta.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 8.6e-114 380.1 0.2 2.6e-93 312.9 0.0 2.1 2 tr|F7DQI0|F7DQI0_MACMU U6 small nuclear RNA (adenine-(43)-N( 0.0075 15.9 0.1 0.012 15.2 0.1 1.2 1 tr|F7EET4|F7EET4_MACMU SHC binding and spindle associated 1 Domain annotation for each sequence (and alignments): >> tr|F7DQI0|F7DQI0_MACMU U6 small nuclear RNA (adenine-(43)-N(6))-methyltransferase OS=Macaca mulatta OX=9544 GN=METTL # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 312.9 0.0 2.5e-97 2.6e-93 1 195 [. 1 195 [. 1 198 [. 0.99 2 ! 66.0 0.0 3.9e-22 4.1e-18 248 299 .] 197 248 .. 195 248 .. 0.97 Alignments for each domain: == domain 1 score: 312.9 bits; conditional E-value: 2.5e-97 XXXXXXXXXXXXXXXXX.XXXXXXXXXXXXXXXXXXXXXXXXXS-TTSHHHHHHHHHHHHHHHH----.---TTS----HHHHHHHH CS Methyltransf_10 1 aalksglhprnrhkdky.dfaalikvvPelkqfvilnpkGrvsidfadPlavkalnkallrefygvsewdiPdgfliPpvPgradyi 86 +al++++h+rnr+kdk+ dfa+l++++P++kq+v++n++Grvs++f+dP+av+al+++llre++g+s+ diP ++liP+vP+r++yi tr|F7DQI0|F7DQI0_MACMU 1 MALSKSMHARNRYKDKPpDFAYLASKYPDFKQHVQINLNGRVSLNFKDPEAVRALTCTLLREDFGLSI-DIPLERLIPTVPLRLNYI 86 6899************99**************************************************.****************** PP HHHHHHHT--TCGCS---EEEEES-TTTTHHHHHHHHHH--EEEEEES-HHHHHHHHHHHHHT.T-TTTEEEEE--ST-SSTTTSTT CS Methyltransf_10 87 hyvadllgkadsdisaqvraldiGvGanCiyPllGvteyGWrfvgsevddvslasakaivesnsklsdaiklrkqkekkaifdglie 173 h+v+dl+g++dsd+s+ +r++diG+Ga+CiyPllG+t++GW+f+++evdd+++++ak++ve+n +lsd+ik++k+++k+ ++d+l+e tr|F7DQI0|F7DQI0_MACMU 87 HWVEDLIGHQDSDKSTLRRGIDIGTGASCIYPLLGATLNGWYFLATEVDDMCFNYAKKNVEQN-NLSDLIKVVKVPQKTLLMDALKE 172 ***************************************************************.*********************** PP --S.-EEEEEE-----SS----- CS Methyltransf_10 174 ene.rydvtlCnPPfhaslaeak 195 e+e +yd+++CnPPf+a+++eak tr|F7DQI0|F7DQI0_MACMU 173 ESEiIYDFCMCNPPFFANQLEAK 195 ***99*****************8 PP == domain 2 score: 66.0 bits; conditional E-value: 3.9e-22 GGEEEEEEEESSGGGHHHHHHHHHHTT-SEEEEEEEEETTEEEEEEEEES-- CS Methyltransf_10 248 kqvlWfttlvskkenlaplkeelkklGaskvtvvelaqGqklsrfiaWsfld 299 k +lW++++++kk++laplkeel+++G++kvt++e++qG+++++++aWsf+d tr|F7DQI0|F7DQI0_MACMU 197 KCLLWYSCMLGKKCSLAPLKEELRIQGVPKVTYTEFCQGRTMRWALAWSFYD 248 5689**********************************************97 PP >> tr|F7EET4|F7EET4_MACMU SHC binding and spindle associated 1 OS=Macaca mulatta OX=9544 GN=SHCBP1 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.2 0.1 1.1e-06 0.012 166 285 .. 515 633 .. 502 639 .. 0.81 Alignments for each domain: == domain 1 score: 15.2 bits; conditional E-value: 1.1e-06 SSTTTSTT--S-EEEEEE-----SS---.---------...---.....---------TTTTH.HHHTHHHHHHHHHHHHHHHGGGE CS Methyltransf_10 166 aifdglieenerydvtlCnPPfhaslae.akagserklknlgkanktralpslnfgGqkaely.ceGGevafikklieeskqlakqv 250 +++ +i++ne y v+l P +++l+e a+ g+e ++ l k +t+ +p + g + +el+ c G +++ ++ + s+q++ + tr|F7EET4|F7EET4_MACMU 515 SMVNNIIHNNEGYGVVLVKPTIFSDLQEnAEDGTEENK-AL-KI-QTSGEPDVTEGVDLEELMqCATGKMELYART-DLSEQVEGNC 597 6899****************988776552777777766.78.43.89999***********987999**9999885.6688999999 PP EEEEEEESSGGGHHHHH.HHHHHTT-SEEEEEEEEE CS Methyltransf_10 251 lWfttlvskkenlaplk.eelkklGaskvtvvelaq 285 + l++ +++ ++k ++l +lG+++ ++q tr|F7EET4|F7EET4_MACMU 598 EIVNELIAASTQKGQMKkKRLSELGITQADDNLMSQ 633 999999988888888762578999988766665555 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (299 nodes) Target sequences: 21211 (10888574 residues searched) Passed MSV filter: 440 (0.020744); expected 424.2 (0.02) Passed bias filter: 360 (0.0169723); expected 424.2 (0.02) Passed Vit filter: 22 (0.0010372); expected 21.2 (0.001) Passed Fwd filter: 2 (9.42907e-05); expected 0.2 (1e-05) Initial search space (Z): 21211 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.22u 0.01s 00:00:00.23 Elapsed: 00:00:00.09 # Mc/sec: 36174.26 // [ok]