# 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_15.hmm.txt # target sequence database: proteomes/Petromyzon_marinus.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Methyltransf_15 [M=165] Accession: PF09445.9 Description: RNA cap guanine-N2 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 0.00013 21.4 0.0 0.00031 20.1 0.0 1.7 2 tr|S4RQZ8|S4RQZ8_PETMA THUMP domain containing 2 OS=Petromyz 0.00017 21.0 0.0 0.00021 20.6 0.0 1.3 1 tr|S4RXG9|S4RXG9_PETMA Methyltransferase like 5 OS=Petromyzo 0.00095 18.5 0.0 0.0013 18.0 0.0 1.2 1 tr|S4RE18|S4RE18_PETMA NOL1/NOP2/Sun domain family, member 6 0.0028 17.0 0.0 0.0055 16.0 0.0 1.4 1 tr|S4RHT4|S4RHT4_PETMA tRNA (guanine(26)-N(2))-dimethyltrans ------ inclusion threshold ------ 0.095 12.0 0.0 0.14 11.5 0.0 1.2 1 tr|S4RMF5|S4RMF5_PETMA HemK methyltransferase family member 0.8 9.0 5.3 1.3 8.3 0.1 2.5 3 tr|S4RW35|S4RW35_PETMA Tissue factor pathway inhibitor OS=Pe Domain annotation for each sequence (and alignments): >> tr|S4RQZ8|S4RQZ8_PETMA THUMP domain containing 2 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 ? -2.8 0.0 1.5 3.3e+03 22 48 .. 167 193 .. 158 195 .. 0.71 2 ! 20.1 0.0 1.4e-07 0.00031 3 91 .. 270 355 .. 268 379 .. 0.81 Alignments for each domain: == domain 1 score: -2.8 bits; conditional E-value: 1.5 Methyltransf_15 22 fekviaidinlsrlklaqhNaevYgvs 48 + + +++ + +++ +N + Y+ + tr|S4RQZ8|S4RQZ8_PETMA 167 ICATTDVESFRVSFRCSGNNSQHYNSQ 193 555666777777889999999999876 PP == domain 2 score: 20.1 bits; conditional E-value: 1.4e-07 Methyltransf_15 3 vvlDlfgGaGGnviqfanvfekv..iaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklkleeikvdlvflSpPwGGpsYk 87 vv+D+++G G + a+ +++v + d++ ++l+ a++N ++ + ++++ l ++ v+ l ++i d++++ +P G +Y+ tr|S4RQZ8|S4RQZ8_PETMA 270 VVMDPMCGMGTLLLEAAMAWKTVrcLGTDVDAQQLAHAKENIHFAELDHRIGLLQASVTALPIP--CSAI--DVIICDIPFG-QKYS 351 799************999999872246799*********************9999988766544..4455..9*******96.6787 PP Methyltransf_15 88 kkek 91 +e+ tr|S4RQZ8|S4RQZ8_PETMA 352 SDEE 355 7665 PP >> tr|S4RXG9|S4RXG9_PETMA Methyltransferase like 5 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 ! 20.6 0.0 9.5e-08 0.00021 6 112 .. 58 164 .. 53 203 .. 0.85 Alignments for each domain: == domain 1 score: 20.6 bits; conditional E-value: 9.5e-08 Methyltransf_15 6 DlfgGaGGnviqfanvfe.kviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklkleeikvdlvflSpPwGGpsYkkkek 91 Dl +G+G +i a+ ++ d++++ l+ + hN ++v v++v+ dv l + + +++ + v++ pP G + k + tr|S4RXG9|S4RXG9_PETMA 58 DLGCGCGVLSIGAAMLGAgLCVGFDVDEEALETCLHNVDEFDVG-VVDVVQSDVCALHQP-RWRKA-FHTVVMNPPFGTKQNKGMDM 141 899********999996515789******************996.69******9976555.66666.799********999999988 PP Methyltransf_15 92 fdlek..kleplklekllkesvk 112 l++ +l +++++ l k++++ tr|S4RXG9|S4RXG9_PETMA 142 VFLQTglELANTSVYSLHKTATR 164 88888677888888888887554 PP >> tr|S4RE18|S4RE18_PETMA NOL1/NOP2/Sun domain family, member 6 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 ! 18.0 0.0 5.8e-07 0.0013 3 63 .. 242 304 .. 240 336 .. 0.89 Alignments for each domain: == domain 1 score: 18.0 bits; conditional E-value: 5.8e-07 Methyltransf_15 3 vvlDlfgGaGGnviqfanv...fekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfella 63 vlD+++ GG++ ++a ++ia+d+ +++ +q+Na + +++ ++d ++++a tr|S4RE18|S4RE18_PETMA 242 RVLDMCAAPGGKTTHIATLmgdTGEIIALDKTAAKVQKIQNNASRL-NLKSIQTFCCDSTQAIA 304 58***************99766568******************999.5688*******999776 PP >> tr|S4RHT4|S4RHT4_PETMA tRNA (guanine(26)-N(2))-dimethyltransferase 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 ! 16.0 0.0 2.5e-06 0.0055 4 78 .. 145 222 .. 143 239 .. 0.89 Alignments for each domain: == domain 1 score: 16.0 bits; conditional E-value: 2.5e-06 Methyltransf_15 4 vlDlfgGaGGnviqfanv...fekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklkleeikvdlvflS 78 vl + + +G +i fa++ +++i+ d +++ ++ +++N e+ +v++ vt + d l+ k k ++ ++d+v l tr|S4RHT4|S4RHT4_PETMA 145 VLEALAASGLRSIRFAKEipgIKSIITNDFSKNAVAAINKNIEMNQVQHLVTASYSDASMLMYKSKNHKESFDVVDLD 222 567778889999*****97877899999****************************************999**99664 PP >> tr|S4RMF5|S4RMF5_PETMA HemK methyltransferase family member 1 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 ? 11.5 0.0 6.1e-05 0.14 4 69 .. 62 129 .. 60 148 .. 0.73 Alignments for each domain: == domain 1 score: 11.5 bits; conditional E-value: 6.1e-05 Methyltransf_15 4 vlDlfgGaGGnviqfanv..fekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklklee 69 vlD+ +G+G + ++ + +a+d++ + ++la +Nae +d+ +v+ d++ + + + e tr|S4RMF5|S4RMF5_PETMA 62 VLDIGCGSGAISLAILSAspLASCVAVDKSADAVALARENAEDLSLQDRCSVVQLDITAVGRAAHATE 129 79999999976655543355899*******************************99988544444333 PP >> tr|S4RW35|S4RW35_PETMA Tissue factor pathway inhibitor 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 ? -3.0 0.0 1.8 4e+03 8 22 .. 31 45 .. 30 60 .. 0.72 2 ? 8.3 0.1 0.0006 1.3 8 58 .. 102 154 .. 100 182 .. 0.72 3 ? 2.9 0.4 0.027 60 8 25 .. 196 213 .. 194 220 .. 0.80 Alignments for each domain: == domain 1 score: -3.0 bits; conditional E-value: 1.8 Methyltransf_15 8 fgGaGGnviqfanvf 22 +gG+ Gn ++f ++ tr|S4RW35|S4RW35_PETMA 31 YGGCRGNQNNFLEEA 45 689999999997654 PP == domain 2 score: 8.3 bits; conditional E-value: 0.0006 Methyltransf_15 8 fgGaGGnviqfanvfekviaidinlsrlkla.qhNaevYgv.sdkvtlvlgdv 58 +gG+GGn ++f ++ e ++++ +l +l+ + +h + Y v ++ +l+ + tr|S4RW35|S4RW35_PETMA 102 YGGCGGNFNNFLEETECMVTCGRSLGTLRTCrYHRTLNYLVvNNSCFLISSCW 154 79***********************9998663577777766244444444433 PP == domain 3 score: 2.9 bits; conditional E-value: 0.027 Methyltransf_15 8 fgGaGGnviqfanvfekv 25 +gG+GGn ++f ++ e v tr|S4RW35|S4RW35_PETMA 196 YGGCGGNFNNFLEETECV 213 79********98876655 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (165 nodes) Target sequences: 13539 (5048892 residues searched) Passed MSV filter: 338 (0.0249649); expected 270.8 (0.02) Passed bias filter: 296 (0.0218628); expected 270.8 (0.02) Passed Vit filter: 33 (0.0024374); expected 13.5 (0.001) Passed Fwd filter: 6 (0.000443164); expected 0.1 (1e-05) Initial search space (Z): 13539 [actual number of targets] Domain search space (domZ): 6 [number of targets reported over threshold] # CPU time: 0.08u 0.00s 00:00:00.08 Elapsed: 00:00:00.04 # Mc/sec: 20826.68 // [ok]