# 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/Drosophila_melanogaster.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.7e-82 277.0 0.0 2.5e-82 276.5 0.0 1.2 1 sp|Q9VGY5|MET15_DROME Probable methyltransferase-like prote 0.00012 21.4 0.0 0.00031 20.0 0.0 1.6 2 sp|Q9VIK9|MET13_DROME Methyltransferase-like protein 13 OS= 0.0095 15.1 0.0 0.025 13.7 0.0 1.6 2 tr|Q9VBH7|Q9VBH7_DROME tRNA (adenine(58)-N(1))-methyltransfe ------ inclusion threshold ------ 0.12 11.5 0.0 0.16 11.0 0.0 1.2 1 sp|Q9VZD2|BN3D3_DROME Probable RNA methyltransferase CG1134 Domain annotation for each sequence (and alignments): >> sp|Q9VGY5|MET15_DROME Probable methyltransferase-like protein 15 homolog OS=Drosophila melanogaster OX=7227 GN=CG146 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 276.5 0.0 7.2e-86 2.5e-82 2 308 .. 25 354 .. 24 355 .. 0.88 Alignments for each domain: == domain 1 score: 276.5 bits; conditional E-value: 7.2e-86 XXX--TTTHHHHHHT---TT-EEEES--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT.......---EEEEES-GGGHHHHHHH CS Methyltransf_5 2 leHvpVLlqevvellepkedgvyiDaTlGaggHsealLeqlselrligidrDksalaraker..Lkefekrltlvkgnfkelkealae 87 + HvpVL + ++e+l+p ++g+y D+T+GaggH+++lLe+ +e +++++drD+ a + a+++ +ef+ rl + g+f++l + +e sp|Q9VGY5|MET15_DROME 25 APHVPVLCDTAIEYLQPVPGGTYFDMTFGAGGHTRRLLEKCPEAKVYALDRDPLAHQLARDMseSEEFKGRLIPLLGKFSDLPKLFKE 112 57***************************************************999999875115689999***************** PP TT..TT-EEEEEEE-S--HHHHH-GGG---SSS--BS---S-..TTSS--HHHHHHH--HHHHHHHHHHHH--TTHHHHHHHHHHHHH CS Methyltransf_5 88 an..vtkvdgillDLGvSslqldkaerGFsfaedapLDmRmd..qsaeltAaevlnrlseeeLakilkeygeekfakriakaivekrk 171 ++ ++vdg+l+D G Ss+q+d+a rGFs +d+pLDmRmd +s +tAa+vl +++e +L+kil+ ygeek ak+ia+ +v++r+ sp|Q9VGY5|MET15_DROME 113 HGlaKNSVDGMLFDFGCSSMQFDEAVRGFSISRDGPLDMRMDggHSGGVTAADVLANVEEGDLVKILRMYGEEKAAKKIARGLVDARN 200 **87789***********************************3335567**********************************99665 PP .HS---BHHHHHHHHHHHS-------......-TTHHHHHHHHHHHHTHHHHHHHHHHHHHHHEEEEEEEEEEESSHHHHHHHHHHHH CS Methyltransf_5 172 .kkpiettkelaevikkvvpakekrk....kihpatrvFQAlRIavNdEleeleelleaalellakggrlavisFhSledrivkkvfk 254 +iettk+la++i++++ +++ h at++FQA+RI+vN+El+e++ + a e+l+ +grl+ i+FhSled ivk++++ sp|Q9VGY5|MET15_DROME 201 aLFKIETTKQLADIIENIMDGGTRKDklrrPAHSATKTFQAIRIFVNNELNEINYGMVLANEILRVDGRLVTITFHSLEDTIVKRHIN 288 1579*************998665554233368*******************************************************9 PP HTT------------------.............-EEEEE.E.E---HHHHHH-GGGTT-EEEEEEE CS Methyltransf_5 255 easkkeevPaglpvaseelka.............klklitkkpilpseeeikeNprsrsakLrviek 308 + + lp++ +++ a k + + +i+p ++e+++N rsrsakLr++ k sp|Q9VGY5|MET15_DROME 289 GNVL-AGAANALPLKYSSHYAidepdilesltkkSWKQLHRHVIVPDADEVARNTRSRSAKLRAAVK 354 9887.46677777777666555556644444444667777779********************9865 PP >> sp|Q9VIK9|MET13_DROME Methyltransferase-like protein 13 OS=Drosophila melanogaster OX=7227 GN=CG2614 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -3.9 0.0 1.7 6e+03 212 244 .. 341 373 .. 338 379 .. 0.79 2 ! 20.0 0.0 9.1e-08 0.00031 27 100 .. 475 548 .. 451 559 .. 0.88 Alignments for each domain: == domain 1 score: -3.9 bits; conditional E-value: 1.7 HHTHHHHHHHHHHHHHHHEEEEEEEEEEESSHH CS Methyltransf_5 212 vNdEleeleelleaalellakggrlavisFhSl 244 v + l+e+++ l+ ++++l + g i + Sl sp|Q9VIK9|MET13_DROME 341 VYSTLDEVKQELADSIKNLSPAGLTDQIPYLSL 373 667789999999999999999988877777776 PP == domain 2 score: 20.0 bits; conditional E-value: 9.1e-08 S--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT.....---.EEEEES-GGGHHHHHHHTTTT-EEEEEEE- CS Methyltransf_5 27 aTlGaggHsealLeqlselrligidrDksalarakerLkefekr.ltlvkgnfkelkealaeanvtkvdgillDL 100 + lG+gg l +l++ r++++++D+ +l++a++ ++ ++++ +++v + ++ e ++++++ ++d++l+D sp|Q9VIK9|MET13_DROME 475 VGLGGGGLCSFLHAALPQARITAVEIDPIMLEVAEQYFELKQDKrFHVVIDDGLDFVERCRNEDI-HFDAVLFDV 548 5699*****999999**********************99888666****999999******9999.9******95 PP >> tr|Q9VBH7|Q9VBH7_DROME tRNA (adenine(58)-N(1))-methyltransferase catalytic subunit TRMT61A OS=Drosophila melanogaste # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -2.8 0.0 0.77 2.6e+03 89 106 .. 9 26 .. 5 30 .. 0.83 2 ! 13.7 0.0 7.4e-06 0.025 12 100 .. 103 192 .. 99 196 .. 0.84 Alignments for each domain: == domain 1 score: -2.8 bits; conditional E-value: 0.77 TTT-EEEEEEE-S--HHH CS Methyltransf_5 89 nvtkvdgillDLGvSslq 106 +++k d ++l L vSs++ tr|Q9VBH7|Q9VBH7_DROME 9 HIEKGDVVILYLSVSSMH 26 4889999*********98 PP == domain 2 score: 13.7 bits; conditional E-value: 7.4e-06 HHHHT---TT-EEEES--TTSHHHHHHHTT--TT.EEEEEES-HHHHHHHTT.....---..EEEEES-GGGHHHHHHHTTTT-EEE CS Methyltransf_5 12 vvellepkedgvyiDaTlGaggHsealLeqlsel.rligidrDksalarakerLkefekr..ltlvkgnfkelkealaeanvtkvdg 95 +++ le+++++v + + G+g s+++L++l+ + +l ++d + ++a++ ++++ +t+++ +l + +k+d+ tr|Q9VBH7|Q9VBH7_DROME 103 IIHQLEVRPGAVVVESGTGSGSLSHYFLRALKPTgHLHTFDFHEARADQARDEFRRHGIAdfVTVYHRDVCNLGFGAELD--GKADA 187 7999***************************9999*****************9999874335788887776666554433..57888 PP EEEE- CS Methyltransf_5 96 illDL 100 ++lDL tr|Q9VBH7|Q9VBH7_DROME 188 VFLDL 192 88887 PP >> sp|Q9VZD2|BN3D3_DROME Probable RNA methyltransferase CG11342 OS=Drosophila melanogaster OX=7227 GN=CG11342 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.0 0.0 4.8e-05 0.16 40 116 .. 75 157 .. 48 166 .. 0.77 Alignments for each domain: == domain 1 score: 11.0 bits; conditional E-value: 4.8e-05 TT--TT.EEEEEES-HHHHHHHTT.....---....EEEEE.S-GGGHHHHHHHTTTT-EEEEEEE-S--HHHHH-GGG---S CS Methyltransf_5 40 eqlsel.rligidrDksalarakerLkefekr....ltlvk.gnfkelkealaeanvtkvdgillDLGvSslqldkaerGFsf 116 e+l++ +++g+d+D++ ++ra e ++ ++ + ++ + f+++k++++ +n +k+d+i++ ++l+++++G f sp|Q9VZD2|BN3D3_DROME 75 ERLHRSvKVLGVDIDPRLIQRASEENESPKDVsyacVDVLDdEAFESVKTYMEVNNLEKFDAICCYSITMWIHLNHHDQGLRF 157 445555699*************99776665545555555551568******************98665566788888888776 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (309 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 459 (0.0332946); expected 275.7 (0.02) Passed bias filter: 234 (0.0169737); expected 275.7 (0.02) Passed Vit filter: 22 (0.00159582); expected 13.8 (0.001) Passed Fwd filter: 4 (0.000290149); expected 0.1 (1e-05) Initial search space (Z): 13786 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.17u 0.01s 00:00:00.18 Elapsed: 00:00:00.06 # Mc/sec: 38084.06 // [ok]