# 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/Anolis_carolinensis.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.7e-42 145.0 0.0 4.1e-42 143.7 0.0 1.7 1 tr|H9GPY3|H9GPY3_ANOCA Trimethylguanosine synthase 1 OS=Anol 1e-06 28.6 0.1 2.5e-06 27.3 0.1 1.6 1 tr|H9GDC5|H9GDC5_ANOCA tRNA methyltransferase 2 homolog A OS 0.00028 20.7 0.0 0.00065 19.5 0.0 1.6 1 tr|H9GBK0|H9GBK0_ANOCA tRNA (guanine(37)-N1)-methyltransfera ------ inclusion threshold ------ 0.022 14.5 0.0 0.045 13.5 0.0 1.4 1 tr|H9G883|H9G883_ANOCA tRNA (guanine(26)-N(2))-dimethyltrans 0.03 14.1 0.0 0.052 13.3 0.0 1.3 1 tr|R4GAD1|R4GAD1_ANOCA HemK methyltransferase family member 0.03 14.1 0.0 0.051 13.3 0.0 1.5 1 tr|G1KY28|G1KY28_ANOCA NOP2/Sun RNA methyltransferase family 0.074 12.8 0.0 0.11 12.2 0.0 1.2 1 tr|G1KK62|G1KK62_ANOCA Ubiquinone biosynthesis O-methyltrans 0.098 12.4 0.0 0.19 11.5 0.0 1.4 1 tr|H9GLQ7|H9GLQ7_ANOCA THUMP domain containing 3 OS=Anolis c 0.18 11.6 0.1 0.73 9.6 0.0 2.1 2 tr|G1KAB8|G1KAB8_ANOCA tRNA methyltransferase 11 homolog OS= 1.1 9.0 5.1 6 6.6 0.1 3.1 3 tr|R4GCU4|R4GCU4_ANOCA Tissue factor pathway inhibitor OS=An Domain annotation for each sequence (and alignments): >> tr|H9GPY3|H9GPY3_ANOCA Trimethylguanosine synthase 1 OS=Anolis carolinensis OX=28377 GN=TGS1 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 143.7 0.0 2.2e-45 4.1e-42 3 164 .. 785 936 .. 783 937 .. 0.93 Alignments for each domain: == domain 1 score: 143.7 bits; conditional E-value: 2.2e-45 Methyltransf_15 3 vvlDlfgGaGGnviqfanvfekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklkleeikvdlvflSpPwGGpsYkkk 89 +++D+f+G GGn+iqfa+ ++viaidi++ +++la++NaevYgv+d++++++gd+++l ++lk d+vflSpPwGGp+Y+ tr|H9GPY3|H9GPY3_ANOCA 785 IIVDAFCGVGGNSIQFALAGKRVIAIDIDPVKIRLAHNNAEVYGVADQIEFICGDFMKLASSLK-----GDIVFLSPPWGGPEYTTA 866 89*******************************************************9998854.....488*************** PP Methyltransf_15 90 ekfdlekkleplklekllkesvklskniilfLPrnsdlkqlaaltaevlenaskvkvvklkenglvkallaylge 164 e+fd+++ + p +++k s k+++ni+ fLPrn+d++q+++l +++ kv++ + n+++k+++ay+g+ tr|H9GPY3|H9GPY3_ANOCA 867 EVFDVQTMICPD-GFEVFKLSQKITNNIVYFLPRNADIDQVTSL----AGPGGKVEIEQNFLNNRLKTITAYFGN 936 *****6666665.55689**********************9887....8************************96 PP >> tr|H9GDC5|H9GDC5_ANOCA tRNA methyltransferase 2 homolog A OS=Anolis carolinensis OX=28377 GN=TRMT2A PE=3 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 27.3 0.1 1.4e-09 2.5e-06 2 83 .. 446 527 .. 445 533 .. 0.84 Alignments for each domain: == domain 1 score: 27.3 bits; conditional E-value: 1.4e-09 Methyltransf_15 2 evvlDlfgGaGGnviqfanvfekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklkleeikvdlv.flSpPwGG 83 ++vlD+++G+G i++a+ +kvi+i+ ++ ++ a+ Na++ + ++v++ +g+ ++l+ +l + + +++ ++ pP G tr|H9GDC5|H9GDC5_ANOCA 446 TTVLDVCCGTGTIGITLAKKVKKVIGIELCQEAVEDAKVNAQLNDL-NNVEFYCGKAEDLVPQLVNSLALQNCItVVDPPRAG 527 579****************************************975.78*********9777666555555555267788777 PP >> tr|H9GBK0|H9GBK0_ANOCA tRNA (guanine(37)-N1)-methyltransferase OS=Anolis carolinensis OX=28377 GN=TRMT5 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 19.5 0.0 3.5e-07 0.00065 3 68 .. 302 367 .. 300 395 .. 0.86 Alignments for each domain: == domain 1 score: 19.5 bits; conditional E-value: 3.5e-07 Methyltransf_15 3 vvlDlfgGaGGnviqfanvfekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklkle 68 v+ D+f+G G +i a+ + v+a d n++ k q+N+++ v+ kv++ d e ++ +e tr|H9GBK0|H9GBK0_ANOCA 302 VLFDVFAGVGPFAIPAAKKNCIVFANDLNPESYKWLQRNCKLNKVDTKVQVFNMDGREFMKGPMKE 367 689************************************************999988865554444 PP >> tr|H9G883|H9G883_ANOCA tRNA (guanine(26)-N(2))-dimethyltransferase OS=Anolis carolinensis OX=28377 GN=TRMT1 PE=3 SV= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.5 0.0 2.4e-05 0.045 4 78 .. 142 219 .. 140 228 .. 0.89 Alignments for each domain: == domain 1 score: 13.5 bits; conditional E-value: 2.4e-05 Methyltransf_15 4 vlDlfgGaGGnviqfanv...fekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklkleeikvdlvflS 78 vl + + +G +i fa++ ++v+a d + + ++l+++N ++ +v + vt d ++ ++k +++ +d++ l tr|H9G883|H9G883_ANOCA 142 VLEALAASGLRSIRFAKEvpgLRSVVANDFSSKAVELITQNVQLNDVGHLVTPSFSDARMVMYQCKADKALFDVIDLD 219 677778889999*****9888889**********************************************9**99664 PP >> tr|R4GAD1|R4GAD1_ANOCA HemK methyltransferase family member 1 OS=Anolis carolinensis OX=28377 GN=HEMK1 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.3 0.0 2.8e-05 0.052 4 80 .. 183 261 .. 181 269 .. 0.81 Alignments for each domain: == domain 1 score: 13.3 bits; conditional E-value: 2.8e-05 Methyltransf_15 4 vlDlfgGaGGnviqfanv..fekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklkleeikvdlvflSpP 80 +l + +G+G + + + ++kviaid+ + ++l+++Nae + +++v + + +v + k l vd ++ pP tr|R4GAD1|R4GAD1_ANOCA 183 ILEVGCGSGAIALSLLKKlpHSKVIAIDKLEAAVNLTKENAERLNLQERVSVLHHEVSSSSWKYLLPWGLVDTIISNPP 261 6778888887777776656789*******************************99998876666666555666666666 PP >> tr|G1KY28|G1KY28_ANOCA NOP2/Sun RNA methyltransferase family member 6 OS=Anolis carolinensis OX=28377 GN=NSUN6 PE=3 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.3 0.0 2.8e-05 0.051 3 60 .. 211 270 .. 209 344 .. 0.78 Alignments for each domain: == domain 1 score: 13.3 bits; conditional E-value: 2.8e-05 Methyltransf_15 3 vvlDlfgGaGGnviqfanvfe...kviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfe 60 +lD+++ GG++ ++a + +via+d+ +++k +++Na + +++ + d ++ tr|G1KY28|G1KY28_ANOCA 211 RILDMCAAPGGKTTHLATLMHdegEVIALDKIADKVKKIKQNASSL-KLNCIKTFCYDGTK 270 79***************99553337******************988.34455555555555 PP >> tr|G1KK62|G1KK62_ANOCA Ubiquinone biosynthesis O-methyltransferase, mitochondrial OS=Anolis carolinensis OX=28377 GN # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.2 0.0 6.1e-05 0.11 3 78 .. 79 152 .. 77 155 .. 0.88 Alignments for each domain: == domain 1 score: 12.2 bits; conditional E-value: 6.1e-05 Methyltransf_15 3 vvlDlfgGaGGnviqfanvfekviaidinlsrlklaqhNaevY.gvsdkvtlvlgdvfellaklkleeikvdlvflS 78 +lD+ +G+G + + + +v +id +++++ a+ a + + +++++++ + +++++++ +e+ +d+v++S tr|G1KK62|G1KK62_ANOCA 79 KILDVGCGGGLLTEPLGRLGASVTGIDPLEDNIRTAEVHAAFDpDLAKRIQYKSCALEDIVEEA-SEQ--FDIVVAS 152 58*********************************99777774389**************9994.344..4899988 PP >> tr|H9GLQ7|H9GLQ7_ANOCA THUMP domain containing 3 OS=Anolis carolinensis OX=28377 GN=THUMPD3 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.5 0.0 0.0001 0.19 2 94 .. 344 445 .. 343 465 .. 0.83 Alignments for each domain: == domain 1 score: 11.5 bits; conditional E-value: 0.0001 Methyltransf_15 2 evvlDlfgGaGGnviqfanvfekviai..dinlsrlklaqhNaevYgvsdkvtlvlgdvfe...........llaklkleeikvdlv 75 +++D+++G+G i+ a ++++ + i d n++ +k + +N +++ + + l+++++ d++ tr|H9GLQ7|H9GLQ7_ANOCA 344 DIIIDPMCGTGAIPIEGAAEWPNCYHIagDNNPQAVKRTANNISSL--LKQIRSKESSISAgkpidsiqwdsCSLPLRTGSV--DVI 426 589*****************988764433999**********9998..9999999999888888888877554444555555..9** PP Methyltransf_15 76 flSpPwGGpsYkkkekfdl 94 ++ P G +kk+++dl tr|H9GLQ7|H9GLQ7_ANOCA 427 VTDMPFGKRIGSKKKNWDL 445 *****99999999999998 PP >> tr|G1KAB8|G1KAB8_ANOCA tRNA methyltransferase 11 homolog OS=Anolis carolinensis OX=28377 GN=TRMT11 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 9.6 0.0 0.00039 0.73 3 36 .. 230 263 .. 228 332 .. 0.85 2 ? -2.1 0.0 1.5 2.8e+03 121 146 .. 395 420 .. 389 430 .. 0.81 Alignments for each domain: == domain 1 score: 9.6 bits; conditional E-value: 0.00039 Methyltransf_15 3 vvlDlfgGaGGnviqfanvfekviaidinlsrlk 36 vv D+f G+GG i+ a+ v++ di+ s+++ tr|G1KAB8|G1KAB8_ANOCA 230 VVFDPFIGTGGLLIASAHFGAYVFGTDIDYSTVH 263 789**************************99985 PP == domain 2 score: -2.1 bits; conditional E-value: 1.5 Methyltransf_15 121 LPrnsdlkqlaaltaevlenaskvkv 146 +Pr +lk ++++ + + ++s+ + tr|G1KAB8|G1KAB8_ANOCA 395 IPRHPCLKLISNCEQTLSSHTSRYLI 420 69999999999999999999888654 PP >> tr|R4GCU4|R4GCU4_ANOCA Tissue factor pathway inhibitor OS=Anolis carolinensis OX=28377 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 6.6 0.1 0.0032 6 7 29 .. 65 87 .. 62 90 .. 0.86 2 ? 2.7 0.0 0.05 93 8 21 .. 142 155 .. 140 173 .. 0.85 3 ? -0.3 0.5 0.43 8e+02 8 21 .. 215 228 .. 214 236 .. 0.82 Alignments for each domain: == domain 1 score: 6.6 bits; conditional E-value: 0.0032 Methyltransf_15 7 lfgGaGGnviqfanvfekviaid 29 +++G+GGn ++fa + + ++ tr|R4GCU4|R4GCU4_ANOCA 65 IYSGCGGNKNNFATLVDCLYECE 87 689**********9988777766 PP == domain 2 score: 2.7 bits; conditional E-value: 0.05 Methyltransf_15 8 fgGaGGnviqfanv 21 ++G+GGn +++a + tr|R4GCU4|R4GCU4_ANOCA 142 YSGCGGNYNNYATY 155 79*******99975 PP == domain 3 score: -0.3 bits; conditional E-value: 0.43 Methyltransf_15 8 fgGaGGnviqfanv 21 fgG+ Gn+++f + tr|R4GCU4|R4GCU4_ANOCA 215 FGGCNGNSNNFGTI 228 9*********9765 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (165 nodes) Target sequences: 18528 (9796552 residues searched) Passed MSV filter: 576 (0.0310881); expected 370.6 (0.02) Passed bias filter: 471 (0.025421); expected 370.6 (0.02) Passed Vit filter: 62 (0.00334629); expected 18.5 (0.001) Passed Fwd filter: 11 (0.000593696); expected 0.2 (1e-05) Initial search space (Z): 18528 [actual number of targets] Domain search space (domZ): 10 [number of targets reported over threshold] # CPU time: 0.16u 0.01s 00:00:00.17 Elapsed: 00:00:00.08 # Mc/sec: 20205.39 // [ok]