# 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/Methyltr_RsmB-F.hmm.txt # target sequence database: proteomes/Drosophila_melanogaster.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Methyltr_RsmB-F [M=200] Accession: PF01189.16 Description: 16S rRNA methyltransferase RsmB/F 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.4e-70 236.7 0.3 4.1e-70 235.1 0.0 2.0 2 tr|A1Z909|A1Z909_DROME LD11307p OS=Drosophila melanogaster O 6.9e-39 133.2 0.0 1.2e-38 132.4 0.0 1.4 1 tr|Q8IPS4|Q8IPS4_DROME FI01826p OS=Drosophila melanogaster O 7.2e-36 123.3 0.0 1.1e-34 119.5 0.0 2.2 1 sp|Q9W4M9|NSUN2_DROME tRNA (cytosine(34)-C(5))-methyltransf 1.2e-30 106.3 0.1 2.4e-30 105.3 0.1 1.5 1 tr|Q9VDZ4|Q9VDZ4_DROME LD16340p OS=Drosophila melanogaster O 3.4e-27 95.0 0.0 6.1e-27 94.2 0.0 1.5 1 tr|Q9VPX3|Q9VPX3_DROME LD40271p OS=Drosophila melanogaster O ------ inclusion threshold ------ 0.02 14.2 0.0 0.049 12.9 0.0 1.6 1 sp|Q9W1V3|FBRL_DROME rRNA 2'-O-methyltransferase fibrillar Domain annotation for each sequence (and alignments): >> tr|A1Z909|A1Z909_DROME LD11307p OS=Drosophila melanogaster OX=7227 GN=CG8545-RA PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 235.1 0.0 1.8e-73 4.1e-70 1 200 [] 428 636 .. 428 636 .. 0.94 2 ? -2.6 0.3 1.2 2.8e+03 143 185 .. 649 690 .. 647 699 .. 0.66 Alignments for each domain: == domain 1 score: 235.1 bits; conditional E-value: 1.8e-73 HHHH---TT-EEEECS-TT-HHHHHHHHHTTTTSEEEEEESSCCCHHHHHHHHHHTT--SEEEEES-CCCCHHHHCC--EEEEEEEE CS Methyltr_RsmB-F 1 allLapqeeetilDlcaaPGgKTthiaelmknkgtvvAvdankeRlkrvaenlkrLgvtntivleaDarkidkvlggelfDriLlDa 87 +++Lapqe+e+ilD+c+aPGgK +hia +mkn+g+++A+d n+ R+k++++n +rLg+ n+ v++ D++k+ ++++g fDr+LlDa tr|A1Z909|A1Z909_DROME 428 VMALAPQENERILDMCSAPGGKGSHIASIMKNTGVLFANDSNRDRIKAIVANFHRLGIVNAVVSCEDGTKFRNIMTG--FDRVLLDA 512 589*************************************************************************9..******** PP --CCGGGTTT-TTHHHH.-CCHHCCHHHHHHHHHHHHHHHE....EEEEEEEEEES---CTCTHHHHHHHHHH-TTEEEESSBSTTS CS Methyltr_RsmB-F 88 PCsgtgvirrdpdvklrrkekdiaqlaklqkellsaaidll....kkgGvlvYsTCsvlveEneevieeaLkkredvelvetelsep 170 PC+gtgv+++dp+vk++++e d++++ +lq++ll +aid + +gG++vYsTCsvl+eEne+vi++aLkkr +v+lv+t l + tr|A1Z909|A1Z909_DROME 513 PCTGTGVVSKDPSVKTTKSEVDVQRCYNLQRKLLLTAIDCTdaksSTGGYIVYSTCSVLPEENEWVIDYALKKR-NVKLVPTGLDFG 598 ***************************************99998789**************************9.9******66544 PP TTSSECEEC........T-EEE-TTTTS.-SEEEEEEE CS Methyltr_RsmB-F 171 kglklgeea........gtlrllPhtenkdgffiaklk 200 ++ ++ + t+r++Pht+n+dgf++aklk tr|A1Z909|A1Z909_DROME 599 VEGFTKFRQhrfhpslnLTKRYYPHTHNMDGFYVAKLK 636 333333333689999999******************96 PP == domain 2 score: -2.6 bits; conditional E-value: 1.2 CTCTHHHHHHHHHH-TTEEEESSBSTTSTTSSECEECT-EEE- CS Methyltr_RsmB-F 143 veEneevieeaLkkredvelvetelsepkglklgeeagtlrll 185 eE+e++++ea+++ + +e +e ++ ++e+g ++ll tr|A1Z909|A1Z909_DROME 649 QEEDEKQLDEAIEADATTEAAAE-SEEVEEEVASKEKGPRKLL 690 68999****99999877776554.4454555555555556655 PP >> tr|Q8IPS4|Q8IPS4_DROME FI01826p OS=Drosophila melanogaster OX=7227 GN=Dmel\CG11109 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 132.4 0.0 5.3e-42 1.2e-38 2 200 .] 221 437 .. 220 437 .. 0.77 Alignments for each domain: == domain 1 score: 132.4 bits; conditional E-value: 5.3e-42 HHH---TT-EEEECS-TT-HHHHHHHHHTTTTSEEEEEESSCCCHHHHHHHHHHTT--SEEEEES-CCCCHHHHCC........... CS Methyltr_RsmB-F 2 llLapqeeetilDlcaaPGgKTthiaelmknkgtvvAvdankeRlkrvaenlkrLgvtntivleaDarkidkvlgg........... 77 +L+pq++e+ilD+caaPG KT+hiaelm ++g vvA+d++++R++++ l + v++ + + tr|Q8IPS4|Q8IPS4_DROME 221 RVLDPQPGERILDMCAAPGNKTSHIAELMGDRGSVVALDNSASRVRSMLPKLGHYKSITAHVFNSTEAVAADAPSApvgeftgppfp 307 58*******************************************999776655544444444433222222222223347779*** PP .--EEEEEEEE--CCGGGTTT-TTHH.HH.-CCHHCCHHHHHHHHHHHHHHHEEEEEEEEEEES---CTCTHHHHHHHHHH-TTEEE CS Methyltr_RsmB-F 78 .elfDriLlDaPCsgtgvirrdpdvk.lrrkekdiaqlaklqkellsaaidllkkgGvlvYsTCsvlveEneevieeaLkkredvel 162 e fDriLlDaPCsg+g +p++ ++ k ++++ + q++l++ a++l+++gG+lvYsTC+v+ +E e ++ aL+k + +l tr|Q8IPS4|Q8IPS4_DROME 308 cESFDRILLDAPCSGLGN---RPQLScSIKQVKVLSSYPHIQRRLFEQAVQLVRPGGILVYSTCTVTEDECECLVAWALRKFVELRL 391 ***************995...344430457888999**********************************************99999 PP ESSB..STTSTTSSECEEC....T-EEE-TTTTS.-..SEEEEEEE CS Methyltr_RsmB-F 163 vete..lsepkglklgeea....gtlrllPhtenkd..gffiaklk 200 ++ + +p g ea +r+ P n d gffiak++ tr|Q8IPS4|Q8IPS4_DROME 392 TDATprWGGPGLPLPGFEAcksrLLQRFGPSGANADtvGFFIAKFQ 437 8753113333333344444555477788899888876689999985 PP >> sp|Q9W4M9|NSUN2_DROME tRNA (cytosine(34)-C(5))-methyltransferase OS=Drosophila melanogaster OX=7227 GN=Nsun2 PE=2 SV # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 119.5 0.0 4.6e-38 1.1e-34 2 199 .. 170 429 .. 169 430 .. 0.86 Alignments for each domain: == domain 1 score: 119.5 bits; conditional E-value: 4.6e-38 HHH---TT-EEEECS-TT-HHHHHHHHHTTTT........SEEEEEESSCCCHHHHHHHHHHTT--SEEEEES-CCCCHHHHCC.... CS Methyltr_RsmB-F 2 llLapqeeetilDlcaaPGgKTthiaelmknk........gtvvAvdankeRlkrvaenlkrLgvtntivleaDarkidkvlgg.... 77 ++L+ ++ +++lD+caaPG+KT+ + e ++ + g v+A+d++++R +++ krL+ + v++ D + +++ ++ sp|Q9W4M9|NSUN2_DROME 170 IVLDVRPTDKVLDMCAAPGSKTAQLIEALHAApeehkippGFVLANDVDNNRCYMLVHQAKRLNSPCLLVTNHDSSVFPNLVTTkpdg 257 7899*******************99887543212223333789***********************************9999888999 PP ....--EEEEEEEE--CCGGGTTT-TTHHHH.-CCHHCCHHHHHHHHHHHHHHHEEEEEEEEEEES---CTCTHHHHHHHHHH-T.TE CS Methyltr_RsmB-F 78 ....elfDriLlDaPCsgtgvirrdpdvklrrkekdiaqlaklqkellsaaidllkkgGvlvYsTCsvlveEneevieeaLkkre.dv 160 +fD+iL D+PCsg g++r++pd+ l+ + ++ +l+ q +++ +++l gG lvYsTCs+++ Ene+v+++++k+ + sp|Q9W4M9|NSUN2_DROME 258 skaiLKFDKILCDVPCSGDGTLRKNPDIWLKWNLAQAYNLHGIQYRIVRRGAEMLEVGGRLVYSTCSLNPIENEAVLQRIIKDADgAL 345 999889**********************7788889999*******************************************9986234 PP EEESSBSTTSTTSSECEEC..............................................T-EEE-TTTTS.-SEEEEEE CS Methyltr_RsmB-F 161 elvetelsepkglklgeea..............................................gtlrllPhtenkdgffiakl 199 elv+ glk++ + + lr+lPh ++ gff+a l sp|Q9W4M9|NSUN2_DROME 346 ELVDAGHL-VPGLKYKPGMtdwklatkevdqiftrfeevpeslhtiirpgmfplpademakigleKCLRVLPHLQDSGGFFVAVL 429 66666332.2344444444559999999999999999999999999999999999999999999999999999999999999987 PP >> tr|Q9VDZ4|Q9VDZ4_DROME LD16340p OS=Drosophila melanogaster OX=7227 GN=Nsun5 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 105.3 0.1 1e-33 2.4e-30 3 199 .. 222 424 .. 220 425 .. 0.90 Alignments for each domain: == domain 1 score: 105.3 bits; conditional E-value: 1e-33 HH---TT-EEEECS-TT-HHHHHHHHHTTTTSEEEEEESSCCCHHHHHHHHHHTT--SEEEEES-CCCC.HHHHCC--EEEEEEEE- CS Methyltr_RsmB-F 3 lLapqeeetilDlcaaPGgKTthiaelmknkgtvvAvdankeRlkrvaenlkrLgvtntivleaDarki.dkvlggelfDriLlDaP 88 lLap + t+lD+caaPG KT hi ++m+nkg+++ v++++ R + + e +k g+ ++ + +Da + ++ ++ iL+D tr|Q9VDZ4|Q9VDZ4_DROME 222 LLAPPSGATVLDMCAAPGMKTVHICNVMQNKGCIYSVEQDHVRYNTLCEITKDAGCDIVKPILGDALNLtPERFP--DVEYILVDPS 306 799******************************************************999999998776144444..599******* PP -CCGGGTTT-TTHHHH.-CCHHCCHHHHHHHHHHHHHHHEEEEEEEEEEES---CTCTHHHHHHHHHH-TTEEEESSBST.TSTTSS CS Methyltr_RsmB-F 89 CsgtgvirrdpdvklrrkekdiaqlaklqkellsaaidllkkgGvlvYsTCsvlveEneevieeaLkkredvelvetels.epkglk 174 Csg g+ +r +++k +++l+ lq ++ls a+ + + Y+TCs+ eEne+v+++ L+ ++ +l + + k + tr|Q9VDZ4|Q9VDZ4_DROME 307 CSGSGMQNRMTVCDEPKEDKRLQKLQGLQIKILSHAMGAFPNVKRIAYCTCSLWKEENEQVVQRCLQLNPSFKLLSCKKAlRNKWHN 393 ****************************************9999*****************************99886441445566 PP ECEEC......T-EEE-TTTTS.-SEEEEEE CS Methyltr_RsmB-F 175 lgeea......gtlrllPhtenkdgffiakl 199 +g + + l P+++ +dg+f+a + tr|Q9VDZ4|Q9VDZ4_DROME 394 VGDKDypnigkNVLYCQPDSDLTDGIFLALF 424 6666689987777888899*********987 PP >> tr|Q9VPX3|Q9VPX3_DROME LD40271p OS=Drosophila melanogaster OX=7227 GN=l(2)10685 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 94.2 0.0 2.6e-30 6.1e-27 3 158 .. 285 446 .. 283 498 .. 0.83 Alignments for each domain: == domain 1 score: 94.2 bits; conditional E-value: 2.6e-30 HH---TT-EEEECS-TT-HHHHHHHHHTTTTSEEEEEESSCCCHHHHHHHHH........HTT--SEEEEES-CCCCHHHHCC--EE CS Methyltr_RsmB-F 3 lLapqeeetilDlcaaPGgKTthiaelmknkgtvvAvdankeRlkrvaenlk........rLgvtntivleaDarkidkvlggelfD 81 +L+ +++e++lD+ca+PGgK ++++ + ++v +d++++Rl+++ + ++ r + +i+++ Dar+ d++ fD tr|Q9VPX3|Q9VPX3_DROME 285 FLQVKPGERVLDACASPGGKSLLMLQT-LHLDHLVCNDIQESRLNKLRKVMQeylfdykeRWAGKRLILSQSDARNLDQYEQ---FD 367 79999********************99.6889999**********9987765222222226677799**********99866...** PP EEEEEE--CCGGGTTT-TT..HHHH.-CCHHCCHHHHHHHHHHHHHHHEEEEEEEEEEES---CTCTHHHHHHHHHH-T CS Methyltr_RsmB-F 82 riLlDaPCsgtgvirrdpd..vklrrkekdiaqlaklqkellsaaidllkkgGvlvYsTCsvlveEneevieeaLkkre 158 +iL+D+PC++ + ++d + + k+ ++ +lq +l + ll++gG lvYsTCs+ + n v++ aL+k tr|Q9VPX3|Q9VPX3_DROME 368 KILVDVPCTTDRHVLNEQDnnIFKPTRIKERLRIPELQAGILANCLRLLRPGGSLVYSTCSLSPIQNDGVVHMALQKVF 446 *********9765554444015444455677889******************************************965 PP >> sp|Q9W1V3|FBRL_DROME rRNA 2'-O-methyltransferase fibrillarin OS=Drosophila melanogaster OX=7227 GN=Fib PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.9 0.0 2.1e-05 0.049 10 87 .. 186 262 .. 179 278 .. 0.77 Alignments for each domain: == domain 1 score: 12.9 bits; conditional E-value: 2.1e-05 -EEEECS-TT-HHHHHHHHHTTTTSEEEEEESSCCCHHHHHHHHHHTT--SEEEEES-CCCCHHHHCC.--EEEEEEEE CS Methyltr_RsmB-F 10 etilDlcaaPGgKTthiaelmknkgtvvAvdankeRlkrvaenlkrLgvtntivleaDarkidkvlgg.elfDriLlDa 87 ++l l aa G +h+++++ ++g v+Av+ +++ + +++ k+ tn+i + Dar+++k+ ++D i D+ sp|Q9W1V3|FBRL_DROME 186 SKVLYLGAASGTTVSHVSDVVGPEGLVYAVEFSHRSGRDLINVAKK--RTNIIPIIEDARHPHKYRMLvGMVDTIFADV 262 6789999************************987665555544433..3456778889999998754435888888876 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (200 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 548 (0.0397505); expected 275.7 (0.02) Passed bias filter: 344 (0.0249529); expected 275.7 (0.02) Passed Vit filter: 43 (0.00311911); expected 13.8 (0.001) Passed Fwd filter: 6 (0.000435224); expected 0.1 (1e-05) Initial search space (Z): 13786 [actual number of targets] Domain search space (domZ): 6 [number of targets reported over threshold] # CPU time: 0.13u 0.01s 00:00:00.14 Elapsed: 00:00:00.06 # Mc/sec: 24649.88 // [ok]