# 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/Petromyzon_marinus.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.2e-74 249.0 0.0 4.2e-74 248.6 0.0 1.1 1 tr|S4R939|S4R939_PETMA Methyltransferase like 9 OS=Petromyzo 0.00073 18.3 0.3 0.037 12.7 0.0 2.5 3 tr|S4R4I3|S4R4I3_PETMA Phosphoethanolamine methyltransferase ------ inclusion threshold ------ 0.14 10.8 0.0 0.22 10.2 0.0 1.2 1 tr|S4RW71|S4RW71_PETMA Protein arginine methyltransferase 2 Domain annotation for each sequence (and alignments): >> tr|S4R939|S4R939_PETMA Methyltransferase like 9 OS=Petromyzon marinus OX=7757 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 248.6 0.0 9.2e-78 4.2e-74 61 265 .] 118 319 .. 113 319 .. 0.96 Alignments for each domain: == domain 1 score: 248.6 bits; conditional E-value: 9.2e-78 DREV 61 inGllgrGsmfvlseeqfrkllaiekespasdqedkklldlGaGdGeitkrvaptfeevyatelsetmrdrlkkkkynvlsalelqe 147 +n +lgrGsmfvls eqfr+ll i +s ++ lldlGaGdGe+t ++ +f++v+ te+s tm++rl+k+++ vl++ e++ tr|S4R939|S4R939_PETMA 118 LNPVLGRGSMFVLSAEQFRRLLGIP----DS-WRADSLLDLGAGDGEVTAIMGAHFKQVFVTEVSPTMQWRLRKRDFTVLGVTEWNI 199 6899********************8....44.88899************************************************85 PP DREV 148 t.dvkldlilclnlldrcfdplklledirsalep.ngrvivalvlPvshyvetntgksekpsellevkGrsyeeevarli.evlkka 231 d +dl++clnlldrc pl+ll+++r++l+p gr ++a+vlP+ +yve g++e+pse + v+G s+ee+ ++l+ +vl++ tr|S4R939|S4R939_PETMA 200 VyDRPYDLVSCLNLLDRCERPLTLLREVRQVLRPgSGRLLLAVVLPYRPYVERG-GRWEQPSECVPVHGDSWEEQLNSLViDVLEPN 285 43889*****************************99****************86.5**********************998****** PP DREV 232 gfeveaftrlPylceGdlnnsfyylddaifvlkk 265 gf+ve+ftr+PylceGdl+n++y+ldda++vlk+ tr|S4R939|S4R939_PETMA 286 GFQVESFTRVPYLCEGDLYNDHYVLDDAVLVLKP 319 ********************************97 PP >> tr|S4R4I3|S4R4I3_PETMA Phosphoethanolamine methyltransferase OS=Petromyzon marinus OX=7757 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 12.7 0.0 8.3e-06 0.037 83 144 .. 39 100 .. 15 140 .. 0.80 2 ! 3.1 0.0 0.0073 33 94 130 .. 283 320 .. 273 387 .. 0.82 3 ? -2.9 0.0 0.47 2.1e+03 71 114 .. 409 453 .. 398 467 .. 0.63 Alignments for each domain: == domain 1 score: 12.7 bits; conditional E-value: 8.3e-06 DREV 83 aiekespasdqedkklldlGaGdGeitkrvaptfeevyatelsetmrdrlkkkkynvlsale 144 + e s ++ +d k+l+lGaG G t ++a++ +v a + +t d+ + ++ + + + tr|S4R4I3|S4R4I3_PETMA 39 KPEILSLLPNIDDAKVLELGAGIGRFTGHLAKQARHVTAVDFMQTFLDKNRELNSDHVNIKF 100 455556678889********************************999998877766555544 PP == domain 2 score: 3.1 bits; conditional E-value: 0.0073 DREV 94 edkklldlGaGdGeitkrvaptfe.evyatelsetmrd 130 +k+ld+G G G +a+ f+ +v +ls++m + tr|S4R4I3|S4R4I3_PETMA 283 AGQKVLDVGCGIGGGNFYMAKAFKvDVLGMDLSSNMVE 320 5689**********999**9999736999****99976 PP == domain 3 score: -2.9 bits; conditional E-value: 0.47 DREV 71 fvlseeqfrkllaiekesp.asdqedkklldlGaGdGeitkrvap 114 ++s ++ k+l +++ ++ ++ + +d+ + e tk +++ tr|S4R4I3|S4R4I3_PETMA 409 ILYSVSKYGKFLEEAGFQNvRAEDRTDQFVDILNKELERTKSIQD 453 466777777777766666413445556666666665555555554 PP >> tr|S4RW71|S4RW71_PETMA Protein arginine methyltransferase 2 OS=Petromyzon marinus OX=7757 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.2 0.0 4.9e-05 0.22 94 129 .. 134 171 .. 123 186 .. 0.81 Alignments for each domain: == domain 1 score: 10.2 bits; conditional E-value: 4.9e-05 DREV 94 edkklldlGaGdGeitkrvaptfee..vyatelsetmr 129 + k +ld+G G G i+ a++ + vya e+se r tr|S4RW71|S4RW71_PETMA 134 RGKAVLDVGCGSGIISMLCAKHAQPsrVYAVEASELAR 171 5688**************999986445*******8655 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (265 nodes) Target sequences: 13539 (5048892 residues searched) Passed MSV filter: 198 (0.0146244); expected 270.8 (0.02) Passed bias filter: 153 (0.0113007); expected 270.8 (0.02) Passed Vit filter: 25 (0.00184652); expected 13.5 (0.001) Passed Fwd filter: 3 (0.000221582); expected 0.1 (1e-05) Initial search space (Z): 13539 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.09u 0.00s 00:00:00.09 Elapsed: 00:00:00.03 # Mc/sec: 44598.55 // [ok]