# 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_5.hmm.txt # target sequence database: proteomes/Macaca_mulatta.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Methyltransf_5 [M=309] Accession: PF01795.18 Description: MraW methylase family 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.6e-79 267.2 0.0 4.2e-79 266.5 0.0 1.3 1 tr|F7GN12|F7GN12_MACMU Uncharacterized protein OS=Macac 2.1e-43 149.3 0.0 2.8e-43 148.8 0.0 1.2 1 tr|A0A1D5Q2S8|A0A1D5Q2S8_MACMU Uncharacterized protein OS=Macac Domain annotation for each sequence (and alignments): >> tr|F7GN12|F7GN12_MACMU Uncharacterized protein OS=Macaca mulatta OX=9544 GN=METTL15 PE=3 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 266.5 0.0 3.9e-83 4.2e-79 3 258 .. 71 371 .. 69 400 .. 0.93 Alignments for each domain: == domain 1 score: 266.5 bits; conditional E-value: 3.9e-83 XX--TTTHHHHHHT--.................................-TT-EEEES--TTSHHHHHHHTT--TTEEEEEES-HHH CS Methyltransf_5 3 eHvpVLlqevvellep.................................kedgvyiDaTlGaggHsealLeqlselrligidrDksa 56 H+pV+++evv+ l+p + +++D+T+G+ggH++a+L++ s+ l+++drD++a tr|F7GN12|F7GN12_MACMU 71 LHIPVMVDEVVHCLSPqkgqisfmcivlfaypnnlvmwtrmvhlhtrklRLREIFLDMTFGSGGHTKAILQKESDIVLYALDRDPTA 157 6**********************************************99999*********************************** PP HHHHTT.....---EEEEES-GGGHHHHHHHTTTT..-EEEEEEE-S--HHHHH-GGG---SSS--BS---S-TTS..S-.-HHHHH CS Methyltransf_5 57 larakerLkefekrltlvkgnfkelkealaeanvt..kvdgillDLGvSslqldkaerGFsfaedapLDmRmdqsa..el.tAaevl 138 a a++ + + k+++ v g+f++ ++ l++a+v+ ++dg+l+DLG Ss+qld++erGFs ++d+pLDmRmd + ++ tAa+v+ tr|F7GN12|F7GN12_MACMU 158 YALAEHLSELYPKQIRAVLGQFSQAEALLTKAGVQpgTFDGVLMDLGCSSMQLDTPERGFSLRKDGPLDMRMDGGRypDMpTAADVV 244 ******77788899*******************76669***********************************87622346****** PP HH--HHHHHHHHHHHH--TTHHHHHHHHHHHHHHS---BHHHHHHHHHHHS-------.........-TTHHHHHHHHHHHHTHHHH CS Methyltransf_5 139 nrlseeeLakilkeygeekfakriakaivekrkkkpiettkelaevikkvvpakekrk.......kihpatrvFQAlRIavNdElee 218 n l++ +La+il++ygeek+a++ia aiv++r+ pi++t++la++++ ++p+++ + ++h at++FQAlRI+vN+El+e tr|F7GN12|F7GN12_MACMU 245 NALDQPALASILRTYGEEKHARKIASAIVQARSIYPITRTQQLASIVAGAFPPSAIYArkdllqrSTHIATKTFQALRIFVNNELNE 331 ****************************************************98876666888889********************* PP HHHHHHHHHHHEEEEEEEEEEESSHHHHHHHHHHHHHTT- CS Methyltransf_5 219 leelleaalellakggrlavisFhSledrivkkvfkeask 258 l + l++a + l++ggrl+ +sFhSledri+k+++ s tr|F7GN12|F7GN12_MACMU 332 LYTGLKTAQKFLRPGGRLVALSFHSLEDRIIKRFLLGISM 371 *********************************9977666 PP >> tr|A0A1D5Q2S8|A0A1D5Q2S8_MACMU Uncharacterized protein OS=Macaca mulatta OX=9544 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 148.8 0.0 2.7e-47 2.8e-43 3 255 .. 57 291 .. 55 331 .. 0.87 Alignments for each domain: == domain 1 score: 148.8 bits; conditional E-value: 2.7e-47 XX--TTTHHHHHHT---TT-EEEES--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT.....---EEEEES-GGGH CS Methyltransf_5 3 eHvpVLlqevvellepkedgvyiDaTlGaggHsealLeqlselrligidrDksalarakerLkefekrltlvkgnfkel 81 H+pV+++evv+ +p++ ++++D+T+ a+L++ s+ l ++drD++a a +++ L +l + g+f++ tr|A0A1D5Q2S8|A0A1D5Q2S8_MACMU 57 LHIPVMVDEVVHWSSPQKRQIFLDMTF-------AILQKESDIVLWALDRDPTAYAFTEH-L----SELYPMLGQFSQA 123 6*************************9.......8*******************999877.4....4466789999999 PP HHHHHHTTTT..-EEEEEEE-S--HHHHH-GGG---SSS--BS---S-TTSS-...-HHHHHHH--HHHHHHHHHHHH- CS Methyltransf_5 82 kealaeanvt..kvdgillDLGvSslqldkaerGFsfaedapLDmRmdqsael...tAaevlnrlseeeLakilkeyge 155 ++ l +a+v+ ++dg+ +DLG S + Fs ++ +pLDmRmd s+ tAa+v n l++++La il++y e tr|A0A1D5Q2S8|A0A1D5Q2S8_MACMU 124 EALLMKAGVQpgTFDGVRMDLGCSYVD-------FSLRKGGPLDMRMDGSRYPdtpTAADVFNALDQQALACILRTYEE 195 9999999976559***********997.......79999*********99743335*********************** PP -TTHHHHHHHHHHHHHHS---BHHHHHHHHHHHS-------..-TTHHHHHHHHHHHHTHHHHHHHHHHHHHHHEEEEE CS Methyltransf_5 156 ekfakriakaivekrkkkpiettkelaevikkvvpakekrkkihpatrvFQAlRIavNdEleeleelleaalellakgg 234 ek+ak+ia iv++ + + + l+ + ++ + ++++h at++FQAl I+vN+El+el + l++a + l++gg tr|A0A1D5Q2S8|A0A1D5Q2S8_MACMU 196 EKHAKKIASPIVQQL-ASIVAGAFPLSAIYAR---KDLLQRSTHIATKTFQALHIFVNNELNELYTGLKTAQKFLRPGG 270 ************987.4455677777766655...34445578************************************ PP EEEEEESSHHHHHHHHHHHHH CS Methyltransf_5 235 rlavisFhSledrivkkvfke 255 l+ +sFhSl+dri+k++ + tr|A0A1D5Q2S8|A0A1D5Q2S8_MACMU 271 HLVTVSFHSLQDRIIKRTSSD 291 ****************97654 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (309 nodes) Target sequences: 21211 (10888574 residues searched) Passed MSV filter: 702 (0.033096); expected 424.2 (0.02) Passed bias filter: 382 (0.0180095); expected 424.2 (0.02) Passed Vit filter: 32 (0.00150865); 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.10 # Mc/sec: 33645.69 // [ok]