# 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_4.hmm.txt # target sequence database: proteomes/Drosophila_melanogaster.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Methyltransf_4 [M=173] Accession: PF02390.16 Description: Putative 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.1e-61 206.9 0.0 1.3e-61 206.7 0.0 1.1 1 sp|O77263|TRMB_DROME tRNA (guanine-N(7)-)-methyltransferas 2.1e-07 30.2 0.0 4.7e-07 29.1 0.0 1.5 1 tr|Q9VVV7|Q9VVV7_DROME Uncharacterized protein OS=Drosophila 0.00026 20.2 0.0 0.00091 18.4 0.0 2.0 1 sp|Q8INR6|DOT1L_DROME Histone-lysine N-methyltransferase, H 0.00036 19.7 0.2 0.0011 18.1 0.0 1.8 2 sp|Q9VJ34|GSTCD_DROME Glutathione S-transferase C-terminal 0.00043 19.4 0.0 0.00094 18.3 0.0 1.6 1 sp|Q7K3B9|MET16_DROME U6 small nuclear RNA (adenine-(43)-N( 0.0018 17.4 0.1 0.0036 16.4 0.0 1.6 2 sp|Q9VNH1|BN3D2_DROME Probable RNA methyltransferase CG1239 0.0022 17.1 0.0 0.0045 16.1 0.0 1.5 1 tr|Q9VMD3|Q9VMD3_DROME Uncharacterized protein, isoform A OS 0.0044 16.1 0.0 0.0068 15.5 0.0 1.3 1 tr|Q9VQJ8|Q9VQJ8_DROME Protein-lysine N-methyltransferase CG 0.0062 15.6 0.0 0.016 14.4 0.0 1.7 1 sp|Q7K480|BN3D1_DROME Probable RNA methyltransferase bin3 O 0.0065 15.6 0.0 0.11 11.6 0.0 2.2 2 tr|Q9VQF8|Q9VQF8_DROME GH13185p2 OS=Drosophila melanogaster ------ inclusion threshold ------ 0.022 13.9 0.0 0.044 12.9 0.0 1.5 1 tr|Q9VFB3|Q9VFB3_DROME Arginine methyltransferase 3, isoform 0.027 13.5 0.0 0.049 12.7 0.0 1.4 1 sp|Q9VAQ5|DIM1_DROME Probable dimethyladenosine transferas 0.056 12.5 0.0 0.12 11.5 0.0 1.5 1 tr|Q9VBJ3|Q9VBJ3_DROME Fire dancer OS=Drosophila melanogaste Domain annotation for each sequence (and alignments): >> sp|O77263|TRMB_DROME tRNA (guanine-N(7)-)-methyltransferase OS=Drosophila melanogaster OX=7227 GN=CG4045 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 206.7 0.0 1.2e-64 1.3e-61 1 171 [. 71 249 .. 71 251 .. 0.98 Alignments for each domain: == domain 1 score: 206.7 bits; conditional E-value: 1.2e-64 S-EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHC........TTSSEEEEES-CCCHHHHCS-TTSEEEEEEES- CS Methyltransf_4 1 kapvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkieal........elkNlrilrgdalkvllelledgslekifilFP 81 +++ +++iGcG+G+flv+l +++P+k+ +G+Eirvkv+++v+++i+al ++N+++lr++a+k+l++++ +g+lek+f+l+P sp|O77263|TRMB_DROME 71 QQVSFADIGCGYGGFLVTLGEMFPEKLSIGMEIRVKVSDYVVDRIAALrrrcadtgAYQNIACLRTNAMKYLPNYFVKGQLEKMFFLYP 159 5899************************************************************************************* PP ----SGGGGGGSTTSHHHHHHHHHHEEEEEEEEEEESHHHHHHHHHHHHHCSTTEEESSCCCSCCTCHCHHCCCSSTTCCCCCCTTSS- CS Methyltransf_4 82 DPwfKkrhhKrRllqpellkelarvLkkgGvlelaTDveeyaeemlkvleeeplferlskeeeleeeplsvlriatkyeqkvqskgkei 170 DP+fK+ +hK+R+++++ll+e+a++L+kgG+l+++TDve++++++++++ee+pl+erl++ee+++++++ +l++++++++kv+++++++ sp|O77263|TRMB_DROME 160 DPHFKRAKHKWRIINQALLSEYAYILRKGGLLYTMTDVEDLHKWIVTHMEEHPLYERLTEEEANADPITPKLYQSSEEGAKVVRNKGDH 248 **************************************************************************************998 PP E CS Methyltransf_4 171 y 171 + sp|O77263|TRMB_DROME 249 F 249 7 PP >> tr|Q9VVV7|Q9VVV7_DROME Uncharacterized protein OS=Drosophila melanogaster OX=7227 GN=cg3808 PE=1 SV=3 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 29.1 0.0 4.4e-10 4.7e-07 3 70 .. 417 482 .. 415 486 .. 0.91 Alignments for each domain: == domain 1 score: 29.1 bits; conditional E-value: 4.4e-10 EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHCTTSSEEEEES-CCCHHHHCS-T CS Methyltransf_4 3 pvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkiealelkNlrilrgdalkvllelled 70 ++l+i+cG G++++a Ak+ k+++G+Ei + ++ + ++ea+++kN + + g+a +++++++++ tr|Q9VVV7|Q9VVV7_DROME 417 TTMLDICCGTGTITLAFAKHC--KKVMGVEIVPDAIKDAEFNAEANGIKNAKFFTGNADDFIKSMVRE 482 5799********999999875..889********************************9999999876 PP >> sp|Q8INR6|DOT1L_DROME Histone-lysine N-methyltransferase, H3 lysine-79 specific OS=Drosophila melanogaster OX=7227 G # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 18.4 0.0 8.6e-07 0.00091 3 36 .. 163 196 .. 161 204 .. 0.88 Alignments for each domain: == domain 1 score: 18.4 bits; conditional E-value: 8.6e-07 EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HH CS Methyltransf_4 3 pvvleiGcGkGeflvalAkknPdknflGiEirvk 36 +++++G+G G++++++A +P k+++GiE+ + sp|Q8INR6|DOT1L_DROME 163 DTFIDLGSGVGQVVLQMAGSFPLKTCIGIEKADT 196 689***************************9766 PP >> sp|Q9VJ34|GSTCD_DROME Glutathione S-transferase C-terminal domain-containing protein homolog OS=Drosophila melanogas # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -2.7 0.0 2.5 2.6e+03 106 148 .. 150 193 .. 144 204 .. 0.61 2 ! 18.1 0.0 1e-06 0.0011 5 63 .. 390 447 .. 387 455 .. 0.92 Alignments for each domain: == domain 1 score: -2.7 bits; conditional E-value: 2.5 HEEEEEEEEEEES.HHHHHHHHHHHHHCSTTEEESSCCCSCCTC CS Methyltransf_4 106 vLkkgGvlelaTD.veeyaeemlkvleeeplferlskeeeleee 148 +L+ g vle++ v + ++m + ++ ++ ++++ ++ ++ e+ sp|Q9VJ34|GSTCD_DROME 150 KLHSGQVLEFVPQeVVRFEQHMNEPVRMHNIYKQAREQANQTEN 193 55666666665544666777777777777666666666555554 PP == domain 2 score: 18.1 bits; conditional E-value: 1e-06 EEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHCTTSSEEEEES-CCCH CS Methyltransf_4 5 vleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkiealelkNlrilrgdalkv 63 +++ ++G G++ + lA k P+++++ +E + ++ ++ k+ ++l+l+N + + ++ sp|Q9VJ34|GSTCD_DROME 390 IVDFCSGTGHLAILLALKLPNCTIIVMENKAFSLLQAQKRSNELGLTNCVFYQCNI-DY 447 699*********************************************99999886.44 PP >> sp|Q7K3B9|MET16_DROME U6 small nuclear RNA (adenine-(43)-N(6))-methyltransferase OS=Drosophila melanogaster OX=7227 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 18.3 0.0 8.9e-07 0.00094 6 72 .. 109 177 .. 104 191 .. 0.83 Alignments for each domain: == domain 1 score: 18.3 bits; conditional E-value: 8.9e-07 EEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHCTTSSEE..EEES-CCCHHHHCS-TTS CS Methyltransf_4 6 leiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkiealelkNlr..ilrgdalkvllelledgs 72 ++iGcG++ + + l k+ + +l++E + + ++++ +++++++++ l + + d ++++++++e+++ sp|Q7K3B9|MET16_DROME 109 IDIGCGSSCIYSLLGAKKNGWHMLALESKPQNIEYAKENVKRNHMESLIevYAQPDNTNIFKSYFEQDQ 177 8******************************************87776422445566799999999865 PP >> sp|Q9VNH1|BN3D2_DROME Probable RNA methyltransferase CG1239 OS=Drosophila melanogaster OX=7227 GN=CG1239 PE=2 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -3.6 0.1 4.6 4.9e+03 131 154 .. 20 43 .. 11 55 .. 0.54 2 ! 16.4 0.0 3.4e-06 0.0036 3 48 .. 111 156 .. 109 182 .. 0.87 Alignments for each domain: == domain 1 score: -3.6 bits; conditional E-value: 4.6 HCSTTEEESSCCCSCCTCHCHHCC CS Methyltransf_4 131 eeeplferlskeeeleeeplsvlr 154 +++ +++++ +e+++e + l++ + sp|Q9VNH1|BN3D2_DROME 20 AKKRKCVITLDEKQVESKRLKKEE 43 445555555555555444444333 PP == domain 2 score: 16.4 bits; conditional E-value: 3.4e-06 EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHC CS Methyltransf_4 3 pvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkieal 48 + +l+iGc +G++ +++A k+ k+++G++i + +++ + k++ +l sp|Q9VNH1|BN3D2_DROME 111 KQLLDIGCNSGHLSIQIARKFEVKSLVGLDIDRGLINDAQKTVSHL 156 6789************************************999877 PP >> tr|Q9VMD3|Q9VMD3_DROME Uncharacterized protein, isoform A OS=Drosophila melanogaster OX=7227 GN=HemK1 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 16.1 0.0 4.3e-06 0.0045 3 73 .. 147 220 .. 145 232 .. 0.79 Alignments for each domain: == domain 1 score: 16.1 bits; conditional E-value: 4.3e-06 EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHCTTSS.EEEEES-C..CCHHHHCS-TTSE CS Methyltransf_4 3 pvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkiealelkN.lrilrgda..lkvllelledgsl 73 + +le+GcG+G + + + + P++ +iE ++ + + ++++ l+l N ++++++ k+l+e+l+d++ tr|Q9VMD3|Q9VMD3_DROME 147 VDLLEVGCGSGAMSLSMLHSLPQVVATAIERSKAATVLAAENAKMLGLLNrFEVHNHTMeeDKYLPEALKDKKY 220 6789*****************************99999999999996666333333322225788888888765 PP >> tr|Q9VQJ8|Q9VQJ8_DROME Protein-lysine N-methyltransferase CG9643 OS=Drosophila melanogaster OX=7227 GN=Dmel\CG9643 P # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.5 0.0 6.4e-06 0.0068 3 49 .. 61 107 .. 59 122 .. 0.88 Alignments for each domain: == domain 1 score: 15.5 bits; conditional E-value: 6.4e-06 EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHCT CS Methyltransf_4 3 pvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkieale 49 + vl++GcG+G flv lA++ + ++ G++ + k v+ + + +e+++ tr|Q9VQJ8|Q9VQJ8_DROME 61 SRVLDLGCGNGMFLVGLANEGFTGDLTGVDYSPKAVELAQNIAEDNK 107 5689****************************999998888777664 PP >> sp|Q7K480|BN3D1_DROME Probable RNA methyltransferase bin3 OS=Drosophila melanogaster OX=7227 GN=bin3 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 14.4 0.0 1.5e-05 0.016 2 42 .. 821 861 .. 820 865 .. 0.91 Alignments for each domain: == domain 1 score: 14.4 bits; conditional E-value: 1.5e-05 -EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHH CS Methyltransf_4 2 apvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvl 42 ++ +l+iGc G+++++ A + + k+++Gi+i +++va++ sp|Q7K480|BN3D1_DROME 821 NKDILDIGCNVGHMTITVARHLAPKTIVGIDIDRELVARAR 861 56799*******************************99875 PP >> tr|Q9VQF8|Q9VQF8_DROME GH13185p2 OS=Drosophila melanogaster OX=7227 GN=HemK2 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 11.6 0.0 0.0001 0.11 3 61 .. 44 105 .. 42 120 .. 0.84 2 ? 1.9 0.0 0.099 1e+02 94 117 .. 156 179 .. 152 220 .. 0.75 Alignments for each domain: == domain 1 score: 11.6 bits; conditional E-value: 0.0001 EEEEEET-TTSHHHHHHHHHSTT.SEEEEEES-HHHHHHHHHHHHHC..TTSSEEEEES-CC CS Methyltransf_4 3 pvvleiGcGkGeflvalAkknPd.knflGiEirvkvvakvlkkieal..elkNlrilrgdal 61 ++++e+G+G+G ++ alAkk + ++l+ +i+ k ++ +++ ++ +l +r+ dal tr|Q9VQF8|Q9VQF8_DROME 44 SLCVELGSGSGVIITALAKKLAGfSLCLATDINPKACNATRRTATRNgaRLDSIRCSLADAL 105 689*****************98725799****999999999888888778888888888884 PP == domain 2 score: 1.9 bits; conditional E-value: 0.099 TTSHHHHHHHHHHEEEEEEEEEEE CS Methyltransf_4 94 llqpellkelarvLkkgGvlelaT 117 +++ llk+l ++L ++Gvl+l tr|Q9VQF8|Q9VQF8_DROME 156 RVTDILLKQLDDILSPRGVLYLLL 179 4567899*************9865 PP >> tr|Q9VFB3|Q9VFB3_DROME Arginine methyltransferase 3, isoform A OS=Drosophila melanogaster OX=7227 GN=Art3 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.9 0.0 4.1e-05 0.044 3 58 .. 245 298 .. 243 303 .. 0.86 Alignments for each domain: == domain 1 score: 12.9 bits; conditional E-value: 4.1e-05 EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHCTTSSEEEEES CS Methyltransf_4 3 pvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkiealelkNlrilrg 58 ++vl++GcG G + + A k ++Gi+ +++v ++++ i++++++N+++++g tr|Q9VFB3|Q9VFB3_DROME 245 KTVLDVGCGTGIL-SIFASKAGAARVVGIDN-SDIVYTAMDIIRKNKVENVELIKG 298 689*******965.55788888899****96.689999*************99976 PP >> sp|Q9VAQ5|DIM1_DROME Probable dimethyladenosine transferase OS=Drosophila melanogaster OX=7227 GN=CG11837 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.7 0.0 4.6e-05 0.049 3 51 .. 51 97 .. 49 113 .. 0.85 Alignments for each domain: == domain 1 score: 12.7 bits; conditional E-value: 4.6e-05 EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHCTTS CS Methyltransf_4 3 pvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkiealelk 51 vvleiG G G+++v++ ++ k++++ Ei +++ a++ k+++a l+ sp|Q9VAQ5|DIM1_DROME 51 DVVLEIGPGTGNMTVRMLER--AKKVIACEIDTRLAAELQKRVQATPLQ 97 589**************876..5889*****************998554 PP >> tr|Q9VBJ3|Q9VBJ3_DROME Fire dancer OS=Drosophila melanogaster OX=7227 GN=fid PE=4 SV=3 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.5 0.0 0.00011 0.12 3 55 .. 163 213 .. 161 230 .. 0.81 Alignments for each domain: == domain 1 score: 11.5 bits; conditional E-value: 0.00011 EEEEEET-TTSHHHHHHHHHSTTSEEEEEES-HHHHHHHHHHHHHC.TTSSEEE CS Methyltransf_4 3 pvvleiGcGkGeflvalAkknPdknflGiEirvkvvakvlkkieal.elkNlri 55 +vv+++GcG+G+ l + nP + +G+E ++++ + +k ++ + Nl++ tr|Q9VBJ3|Q9VBJ3_DROME 163 SVVCDVGCGSGRYLT---QCNPAICTIGVERCYRLSKVAHEKGGEVaLCDNLEL 213 689*********976...679**********99999999988887734466665 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (173 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 531 (0.0385173); expected 275.7 (0.02) Passed bias filter: 308 (0.0223415); expected 275.7 (0.02) Passed Vit filter: 38 (0.00275642); expected 13.8 (0.001) Passed Fwd filter: 13 (0.000942986); expected 0.1 (1e-05) Initial search space (Z): 13786 [actual number of targets] Domain search space (domZ): 13 [number of targets reported over threshold] # CPU time: 0.12u 0.01s 00:00:00.13 Elapsed: 00:00:00.05 # Mc/sec: 25586.58 // [ok]