# 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/Pan_troglodytes.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 6.8e-137 456.1 0.0 8.9e-137 455.7 0.0 1.1 1 tr|A0A2I3RSR0|A0A2I3RSR0_PANTR U6 small nuclear RNA (adenine-(4 ------ inclusion threshold ------ 0.026 14.2 0.0 0.043 13.5 0.0 1.2 1 tr|A0A2I3SII7|A0A2I3SII7_PANTR SHC SH2-domain binding protein 1 0.034 13.8 0.0 0.075 12.7 0.0 1.5 1 tr|A0A2I3SZZ8|A0A2I3SZZ8_PANTR Transcription factor OS=Pan trog Domain annotation for each sequence (and alignments): >> tr|A0A2I3RSR0|A0A2I3RSR0_PANTR U6 small nuclear RNA (adenine-(43)-N(6))-methyltransferase OS=Pan troglodytes OX=9598 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 455.7 0.0 1.2e-140 8.9e-137 1 299 [] 1 291 [. 1 291 [. 0.99 Alignments for each domain: == domain 1 score: 455.7 bits; conditional E-value: 1.2e-140 XXXXXXXXXXXXXXXXX.XXXXXXXXXXXXXXXXXXXXXXXXXS-TTSHHHHHHHHHHHHHHHH----.---TTS---- CS Methyltransf_10 1 aalksglhprnrhkdky.dfaalikvvPelkqfvilnpkGrvsidfadPlavkalnkallrefygvsewdiPdgfliPp 78 +al++++h+rnr+kdk+ dfa+l++++P++kq+v++n++Grvs++f+dP+av+al+++llre++g+s+ diP ++liP+ tr|A0A2I3RSR0|A0A2I3RSR0_PANTR 1 MALSKSMHARNRYKDKPpDFAYLASKYPDFKQHVQINLNGRVSLNFKDPEAVRALTCTLLREDFGLSI-DIPLERLIPT 78 6899************99**************************************************.********** PP HHHHHHHHHHHHHHHT--TCGCS---EEEEES-TTTTHHHHHHHHHH--EEEEEES-HHHHHHHHHHHHHT.T-TTTEE CS Methyltransf_10 79 vPgradyihyvadllgkadsdisaqvraldiGvGanCiyPllGvteyGWrfvgsevddvslasakaivesnsklsdaik 157 vP+r++yih+v+dl+g++dsd+s+ +r++diG+Ga+CiyPllG+t++GW+f+++evdd+++++ak++ve+n +lsd+ik tr|A0A2I3RSR0|A0A2I3RSR0_PANTR 79 VPLRLNYIHWVEDLIGHQDSDKSTLRRGIDIGTGASCIYPLLGATLNGWYFLATEVDDMCFNYAKKNVEQN-NLSDLIK 156 ***********************************************************************.******* PP EEE--ST-SSTTTSTT--S.-EEEEEE-----SS------------...---.....---------TTTTHHHHTHHHH CS Methyltransf_10 158 lrkqkekkaifdglieene.rydvtlCnPPfhaslaeakagserklknlgkanktralpslnfgGqkaelyceGGevaf 235 ++k+++k+ ++d+l+ee+e +yd+++CnPPf+a+++eak++++r++ +++ +ps++++G+++e+++eGGe++f tr|A0A2I3RSR0|A0A2I3RSR0_PANTR 157 VVKVPQKTLLMDALKEESEiIYDFCMCNPPFFANQLEAKGVNSRNP---RRP-----PPSSVNTGGITEIMAEGGELEF 227 *******************99*************************...776.....899******************* PP HHHHHHHHHHHGGGEEEEEEEESSGGGHHHHHHHHHHTT-SEEEEEEEEETTEEEEEEEEES-- CS Methyltransf_10 236 ikklieeskqlakqvlWfttlvskkenlaplkeelkklGaskvtvvelaqGqklsrfiaWsfld 299 +k++i++s ql+k+++W++++++kk++laplkeel+++G++kvt++e++qG+++++++aWsf+d tr|A0A2I3RSR0|A0A2I3RSR0_PANTR 228 VKRIIHDSLQLKKRLRWYSCMLGKKCSLAPLKEELRIQGVPKVTYTEFCQGRTMRWALAWSFYD 291 **************************************************************97 PP >> tr|A0A2I3SII7|A0A2I3SII7_PANTR SHC SH2-domain binding protein 1 OS=Pan troglodytes OX=9598 GN=SHCBP1 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.5 0.0 5.6e-06 0.043 166 285 .. 530 648 .. 517 654 .. 0.80 Alignments for each domain: == domain 1 score: 13.5 bits; conditional E-value: 5.6e-06 SSTTTSTT--S-EEEEEE-----SS---.---------...---.....---------TTTTH.HHHTHHHHHHHHHHH CS Methyltransf_10 166 aifdglieenerydvtlCnPPfhaslae.akagserklknlgkanktralpslnfgGqkaely.ceGGevafikkliee 242 +++ +i++ne y v+l P +++l+e a+ g+e ++ l k +t+ +p + + +el+ c G +++ +++ + tr|A0A2I3SII7|A0A2I3SII7_PANTR 530 SMVNNIIHNNEGYGVVLVKPTIFSDLQEnAEDGTEENK-AL-KI-QTSGEPDVAERVDLEELIeCATGKMELCARI-DP 604 6899****************988776552777777766.78.44.78888999998888887559******99995.78 PP HHHHGGGEEEEEEEESSGGGHHHHH.HHHHHTT-SEEEEEEEEE CS Methyltransf_10 243 skqlakqvlWfttlvskkenlaplk.eelkklGaskvtvvelaq 285 s+q++ + + l++ +++ ++k ++l +lG+++ ++q tr|A0A2I3SII7|A0A2I3SII7_PANTR 605 SEQVEGNCEIVNELIAASTQKGQIKkKRLSELGITQADDNLMSQ 648 99999999999999998888877752579999988776666655 PP >> tr|A0A2I3SZZ8|A0A2I3SZZ8_PANTR Transcription factor OS=Pan troglodytes OX=9598 GN=TFDP1 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.7 0.0 9.8e-06 0.075 18 77 .. 4 66 .. 2 76 .. 0.89 Alignments for each domain: == domain 1 score: 12.7 bits; conditional E-value: 9.8e-06 XXXXXXXXXXXXXXXXXXXX...XXXXXS-TTSHHHHHHHHHHHHHHHH----.---TTS--- CS Methyltransf_10 18 dfaalikvvPelkqfvilnp...kGrvsidfadPlavkalnkallrefygvsewdiPdgfliP 77 d+a li+ + elk f+ n+ kG vs+ P +v l k ll + +g s+ +i ++ +i tr|A0A2I3SZZ8|A0A2I3SZZ8_PANTR 4 DLAGLIEANGELKVFIDQNLspgKGVVSLVAVHPSTVNPLGKQLLPKTFGQSNVNIAQQVVIG 66 99***************99744568999999*************************9988774 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (299 nodes) Target sequences: 23008 (11942255 residues searched) Passed MSV filter: 460 (0.019993); expected 460.2 (0.02) Passed bias filter: 375 (0.0162987); expected 460.2 (0.02) Passed Vit filter: 26 (0.00113004); expected 23.0 (0.001) Passed Fwd filter: 3 (0.000130389); expected 0.2 (1e-05) Initial search space (Z): 23008 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.23u 0.01s 00:00:00.24 Elapsed: 00:00:00.08 # Mc/sec: 44634.18 // [ok]