# 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/DREV.hmm.txt # target sequence database: proteomes/Drosophila_melanogaster.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: DREV [M=265] Accession: PF05219.11 Description: DREV 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5.8e-132 438.3 0.0 6.8e-132 438.1 0.0 1.0 1 tr|Q9W1H1|Q9W1H1_DROME GH17155p OS=Drosophila melanogaster O ------ inclusion threshold ------ 0.016 13.9 0.0 0.028 13.2 0.0 1.3 1 sp|Q9VAQ5|DIM1_DROME Probable dimethyladenosine transferas 0.04 12.6 0.0 0.15 10.8 0.0 1.8 2 sp|Q9VH48|CARM1_DROME Probable histone-arginine methyltrans 0.082 11.6 0.0 0.1 11.3 0.0 1.4 1 tr|Q9VMD3|Q9VMD3_DROME Uncharacterized protein, isoform A OS Domain annotation for each sequence (and alignments): >> tr|Q9W1H1|Q9W1H1_DROME GH17155p OS=Drosophila melanogaster OX=7227 GN=Dmel\CG5339 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 438.1 0.0 2e-135 6.8e-132 1 265 [] 32 302 .. 32 302 .. 0.99 Alignments for each domain: == domain 1 score: 438.1 bits; conditional E-value: 2e-135 DREV 1 WyvlnrelseslqalfvaslpdaetqsflenskklsanlfiqlyhalassilslfmsktdinGllgrGsmfvlseeqfrkllaieke 87 Wy+ln+el++++q++fv+++pda+t+s+le++k+lsan++++l+hala+sil++fm++tdinG+l+rGsmf+lseeqfrkll+++++ tr|Q9W1H1|Q9W1H1_DROME 32 WYCLNHELPPRFQSKFVPLEPDATTVSWLERTKALSANIWTHLWHALARSILQFFMTQTDINGFLKRGSMFILSEEQFRKLLVAGGF 118 *************************************************************************************** PP DREV 88 spasdqedkklldlGaGdGeitkrvaptfee........vyatelsetmrdrlkkkkynvlsale.lqetdvkldlilclnlldrcf 165 +pas+qe+++lld+GaGdGei+ rva+t++e v+ate+s+tmrdrlkk+++nv+++++ lq +l+lilcln+ldrcf tr|Q9W1H1|Q9W1H1_DROME 119 QPASGQEPVTLLDIGAGDGEISLRVANTVSElsgsaglrVFATEASWTMRDRLKKLNFNVITEIDgLQ----NLELILCLNVLDRCF 201 *****************************************************************999....9************** PP DREV 166 dplklledirsalepngrvivalvlPvshyvetntgksekpsellevkG..rsyeeevarlievlkkagfeveaftrlPylceGdln 250 d++klledi++al+pngr++valvlP++hyvetnt +++++++lle++G +s+eee+ar++e+l+++g++vea+t++PylceGdl+ tr|Q9W1H1|Q9W1H1_DROME 202 DSFKLLEDIHLALAPNGRAVVALVLPYMHYVETNT-SHLPLRPLLENNGrqASFEEEAARFMELLENCGYRVEAWTKAPYLCEGDLH 287 ***********************************.*************8889********************************** PP DREV 251 nsfyylddaifvlkk 265 +sfy+l+d+i+v++k tr|Q9W1H1|Q9W1H1_DROME 288 QSFYWLIDLIVVISK 302 ************986 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 ? 13.2 0.0 8.1e-06 0.028 88 135 .. 43 90 .. 13 94 .. 0.83 Alignments for each domain: == domain 1 score: 13.2 bits; conditional E-value: 8.1e-06 DREV 88 spasdqedkklldlGaGdGeitkrvaptfeevyatelsetmrdrlkkk 135 + a+ + + +l++G G G++t r+ + ++v a e+ +++ l+k+ sp|Q9VAQ5|DIM1_DROME 43 EKAALRATDVVLEIGPGTGNMTVRMLERAKKVIACEIDTRLAAELQKR 90 44555677889********************************99986 PP >> sp|Q9VH48|CARM1_DROME Probable histone-arginine methyltransferase CARMER OS=Drosophila melanogaster OX=7227 GN=Art4 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.8 0.0 4.2e-05 0.15 90 128 .. 175 214 .. 165 233 .. 0.82 2 ? -1.6 0.0 0.26 8.9e+02 181 222 .. 245 286 .. 226 308 .. 0.81 Alignments for each domain: == domain 1 score: 10.8 bits; conditional E-value: 4.2e-05 DREV 90 asdqedkklldlGaGdGeitk.rvaptfeevyatelsetm 128 a d++dk +ld+GaG G ++ v+ +vya e+s+ sp|Q9VH48|CARM1_DROME 175 AVDFQDKIVLDVGAGSGILSFfAVQAGAAKVYAIEASNMA 214 67899*************996257888999******9754 PP == domain 2 score: -1.6 bits; conditional E-value: 0.26 DREV 181 ngrvivalvlPvshyvetntgksekpsellevkGrsyeeeva 222 + vi++ + ++ y e +++ ++ l+ +G+ y + + sp|Q9VH48|CARM1_DROME 245 KVDVIISEPMGYMLYNERMLETYLHARKWLKPQGKMYPTHGD 286 67899999999******99999999******99999987765 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 ? 11.3 0.0 2.9e-05 0.1 76 139 .. 131 191 .. 62 252 .. 0.78 Alignments for each domain: == domain 1 score: 11.3 bits; conditional E-value: 2.9e-05 DREV 76 eqfrkllaiekespasdqedkklldlGaGdGeitkrvaptfeevyatelse.tmrdrlkkkkynv 139 e f +l+ i++ + + +++ ll++G G G ++ + + +v at + + l + ++ tr|Q9VMD3|Q9VMD3_DROME 131 EEFMRLV-IDD-HKN--AKHVDLLEVGCGSGAMSLSMLHSLPQVVATAIERsKAATVLAAENAKM 191 5677765.552.333..4799*************************9875313444455555555 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (265 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 303 (0.0219788); expected 275.7 (0.02) Passed bias filter: 235 (0.0170463); expected 275.7 (0.02) Passed Vit filter: 28 (0.00203105); 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.14u 0.01s 00:00:00.15 Elapsed: 00:00:00.06 # Mc/sec: 32661.09 // [ok]