# 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/UPF0020.hmm.txt # target sequence database: proteomes/Anolis_carolinensis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: UPF0020 [M=197] Accession: PF01170.17 Description: Putative RNA methylase family UPF0020 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.5e-42 145.6 0.0 4.9e-42 143.9 0.0 1.8 1 tr|H9GLQ7|H9GLQ7_ANOCA THUMP domain containing 3 OS=Anolis c 7.8e-25 87.8 0.0 2.5e-24 86.1 0.0 1.8 1 tr|G1K8K6|G1K8K6_ANOCA THUMP domain containing 2 OS=Anolis c 6.9e-13 48.8 0.0 1.7e-12 47.5 0.0 1.6 1 tr|G1KAB8|G1KAB8_ANOCA tRNA methyltransferase 11 homolog OS= 1.3e-06 28.3 0.0 5.4e-06 26.3 0.0 2.0 1 tr|H9GPY3|H9GPY3_ANOCA Trimethylguanosine synthase 1 OS=Anol 1.1e-05 25.3 0.1 2.7e-05 24.0 0.1 1.8 1 tr|H9GNU1|H9GNU1_ANOCA Methyltransferase like 5 OS=Anolis ca 0.00036 20.3 0.0 0.0013 18.5 0.0 1.8 1 tr|R4GAD1|R4GAD1_ANOCA HemK methyltransferase family member 0.0025 17.6 0.0 0.0059 16.4 0.0 1.7 1 tr|H9G6W9|H9G6W9_ANOCA 2-methoxy-6-polyprenyl-1,4-benzoquino ------ inclusion threshold ------ 0.013 15.2 0.4 0.066 12.9 0.4 2.0 1 tr|H9GDC5|H9GDC5_ANOCA tRNA methyltransferase 2 homolog A OS Domain annotation for each sequence (and alignments): >> 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 ! 143.9 0.0 2.1e-45 4.9e-42 2 196 .. 305 500 .. 303 501 .. 0.92 Alignments for each domain: == domain 1 score: 143.9 bits; conditional E-value: 2.1e-45 SSS-SS---SS-HHHHHHHH............HHTT--TTS-EEETT-TTSHHHHHHHHHHTT--TTTT----SSHHHHHHHHHHHH CS UPF0020 2 gyRefngpapLketLAaalv............rlagwkpgepllDPmCGsGtilIEaallganiapgllrefvlelkaeaeeearae 76 ++++f gp++L++tLA++++ rl++ +p + ++DPmCG+G+i+IE+a + n+++++ +++++ tr|H9GLQ7|H9GLQ7_ANOCA 305 NITHF-GPTTLRSTLAYGMLrylcyvtkkqqrRLCDPQPADIIIDPMCGTGAIPIEGAAEWPNCYHIAG-------------DNNPQ 377 55555.5************944444444444448*******************************7666.............88999 PP --EEEEES-HHHHHHHHHHHHHTT-GGGEEEEE--GGG--STTS---EEEE----------HHHHHHHHHHHHHHHHHH-TT...-E CS UPF0020 77 lklygsDldrrvvqgareNaekagvgdliefsqadaakLrlkegevdvivtnpPYGerlgskkalekLYseflrelkrvlrgg.grl 162 + + ++++++++++ r++ ++ g+ i+ q+d ++L+l++g+vdvivt++P+G+r+gskk++++LY+++l e+ r++r + gr+ tr|H9GLQ7|H9GLQ7_ANOCA 378 AVKRTANNISSLLKQIRSKESSISAGKPIDSIQWDSCSLPLRTGSVDVIVTDMPFGKRIGSKKKNWDLYPSCLMEMGRICRPKtGRA 464 9999*****************************************************************************999*** PP EEEEESSHHHHHCC---..ESEEEEEEETTEEEEEE CS UPF0020 163 vlltsenkalekaleer..alrkvkeknvklglllk 196 vllt+++k+++kal+++ +rk + +v++g+l++ tr|H9GLQ7|H9GLQ7_ANOCA 465 VLLTQDRKCFAKALSKMghLWRKSHTIWVNVGGLHA 500 *****************989************9985 PP >> tr|G1K8K6|G1K8K6_ANOCA THUMP domain containing 2 OS=Anolis carolinensis OX=28377 GN=THUMPD2 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 86.1 0.0 1.1e-27 2.5e-24 10 195 .. 251 414 .. 243 416 .. 0.87 Alignments for each domain: == domain 1 score: 86.1 bits; conditional E-value: 1.1e-27 -SS-HHHHHHHHHHTT--TTS-EEETT-TTSHHHHHHHHHHTT--TTTT----SSHHHHHHHHHHHH--EEEEES-HHHHHHHHHHH CS UPF0020 10 apLketLAaalvrlagwkpgepllDPmCGsGtilIEaallganiapgllrefvlelkaeaeeearaelklygsDldrrvvqgareNa 96 L++t+A a+v la++k+g lDP CG Gtil+Eaa+ + ++ g D+ + +qga +N tr|G1K8K6|G1K8K6_ANOCA 251 IGLRSTIAWAMVSLAEIKAGAFVLDPTCGLGTILLEAAMEWPDVN------------------------YLGTDISDAQLQGAYTNI 313 579*****************************************9........................79**************** PP HHTT-GGGEEEEE--GGG--STTS---EEEE----------HHHHHHHHHHHHHHHHHH-TT..-EEEEEESS..HHHHHCC---.E CS UPF0020 97 ekagvgdliefsqadaakLrlkegevdvivtnpPYGerlgskkalekLYseflrelkrvlrgggrlvlltsen..kalekaleer.a 180 ++a++ d+i++ +a +++L+l+++++dv++++ P+G++ + +k++e L + + +el+rvl+ gg lvll+s + ++l++ +++ + tr|G1K8K6|G1K8K6_ANOCA 314 KAANLMDKIQLLKASVTALPLPSESIDVVLADIPFGKKFKITKDME-LLPDVFQELERVLCAGGTLVLLLSRElyMCLAHGINSGpE 399 *****************************************99986.6688889***************998622444554544466 PP SEEEEEEETTEEEEE CS UPF0020 181 lrkvkeknvklglll 195 +++ ++ v ++++ tr|G1K8K6|G1K8K6_ANOCA 400 SENHISVEVANDCTQ 414 666666777766665 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 ! 47.5 0.0 7.2e-16 1.7e-12 3 141 .. 202 333 .. 201 358 .. 0.83 Alignments for each domain: == domain 1 score: 47.5 bits; conditional E-value: 7.2e-16 SS-SS---SS-HHHHHHHHHHTT--TTS-EEETT-TTSHHHHHHHHHHTT--TTTT----SSHHHHHHHHHHHH--EEEEES-HHHH CS UPF0020 3 yRefngpapLketLAaalvrlagwkpgepllDPmCGsGtilIEaallganiapgllrefvlelkaeaeeearaelklygsDldrrvv 89 +R+f g++++ + LA + + ++kp+++ +DP+ G+G++lI a ga++ + + + + ++ r+++k +g D + tr|G1KAB8|G1KAB8_ANOCA 202 KRHFIGNTSMDACLAFIMANHGRVKPNDVVFDPFIGTGGLLIASAHFGAYVF---GTDIDYSTVHGLGKASRKNQKWRGPDENI--- 282 79**************************************************...566778888888899999999*****876... PP HHHHHHHHHTT-GGGE..EEEE--GGG--STTS---EEEE----------HHHH CS UPF0020 90 qgareNaekagvgdli..efsqadaakLrlkegevdvivtnpPYGerlgskkal 141 r N++++g+++ +f++ + + +k +++d+i+t+pPYG+r + ++ tr|G1KAB8|G1KAB8_ANOCA 283 ---RANLRQYGLEKHYldAFVSDASRPIWRKGAQFDAIITDPPYGIREATRRTG 333 ...57*******986532455555566667778************986666554 PP >> 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 ! 26.3 0.0 2.3e-09 5.4e-06 30 132 .. 784 859 .. 767 868 .. 0.87 Alignments for each domain: == domain 1 score: 26.3 bits; conditional E-value: 2.3e-09 S-EEETT-TTSHHHHHHHHHHTT--TTTT----SSHHHHHHHHHHHH--EEEEES-HHHHHHHHHHHHHTT-GGGEEEEE--GGG-- CS UPF0020 30 epllDPmCGsGtilIEaallganiapgllrefvlelkaeaeeearaelklygsDldrrvvqgareNaekagvgdliefsqadaakLr 116 + ++D +CG G+ I al g++++ ++D+d+ ++ a +Nae +gv+d+ief + d+ kL tr|H9GPY3|H9GPY3_ANOCA 784 DIIVDAFCGVGGNSIQFALAGKRVI--------------------------AIDIDPVKIRLAHNNAEVYGVADQIEFICGDFMKLA 844 5677777777777777777777777..........................9*********************************** PP STTS---EEEE----- CS UPF0020 117 lkegevdvivtnpPYG 132 ++ + d++ pP+G tr|H9GPY3|H9GPY3_ANOCA 845 -SSLKGDIVFLSPPWG 859 .88899*********9 PP >> tr|H9GNU1|H9GNU1_ANOCA Methyltransferase like 5 OS=Anolis carolinensis OX=28377 GN=METTL5 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 24.0 0.1 1.2e-08 2.7e-05 79 155 .. 77 148 .. 27 186 .. 0.77 Alignments for each domain: == domain 1 score: 24.0 bits; conditional E-value: 1.2e-08 EEEEES-HHHHHHHHHHHHHTT-GGGEEEEE--GGG--.STTS---EEEE----------HHHHHHHHHHHHHHHHHH CS UPF0020 79 lygsDldrrvvqgareNaekagvgdliefsqadaakLr.lkegevdvivtnpPYGerlgskkalekLYseflrelkrv 155 +g+D+d +++ +Nae+ ++++ i++ q ++++L+ ++ + +d+++ npP+G + ++ ++ +fl+++ ++ tr|H9GNU1|H9GNU1_ANOCA 77 CVGFDIDLDALDVFSKNAEEFELTN-IDMIQCNVCSLPdKMPKIFDTVIMNPPFGTKHNKGMDM-----AFLKTALQM 148 589*******************986.777888******867779***********875554443.....344444443 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 ! 18.5 0.0 5.6e-07 0.0013 76 131 .. 203 262 .. 178 265 .. 0.89 Alignments for each domain: == domain 1 score: 18.5 bits; conditional E-value: 5.6e-07 H--EEEEES-HHHHHHHHHHHHHTT-GGGEEEEE--GGG--....STTS---EEEE---- CS UPF0020 76 elklygsDldrrvvqgareNaekagvgdliefsqadaakLr....lkegevdvivtnpPY 131 + k++++D+ +++v+ ++eNae+ ++++++++ ++ + l+ g vd+i++npPY tr|R4GAD1|R4GAD1_ANOCA 203 HSKVIAIDKLEAAVNLTKENAERLNLQERVSVLHHEVSSSSwkylLPWGLVDTIISNPPY 262 56799*************************9988888776666777899*********** PP >> tr|H9G6W9|H9G6W9_ANOCA 2-methoxy-6-polyprenyl-1,4-benzoquinol methylase, mitochondrial OS=Anolis carolinensis OX=283 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 16.4 0.0 2.5e-06 0.0059 26 124 .. 86 195 .. 79 228 .. 0.68 Alignments for each domain: == domain 1 score: 16.4 bits; conditional E-value: 2.5e-06 --TTS-EEETT-TTSHHHHHHHHHHTT--TTTT----SS...........HHHHHHHHHHHH--EEEEES-HHHHHHHHHHHHHTT- CS UPF0020 26 wkpgepllDPmCGsGtilIEaallganiapgllrefvle...........lkaeaeeearaelklygsDldrrvvqgareNaekagv 101 pg++llD G+G i+ +++ ++r++ l +e e++ +++ D++ ++++ ++++a++ g tr|H9G6W9|H9G6W9_ANOCA 86 PYPGTQLLDVAGGTGDIAFRFLNSVRSHREHQVRQKLRShqnlswqeitkLYQEEEQQSLGGSHVVVCDINTEMLKVGKQKAQHLGY 172 55666666666777766654333333333333333333256666665432222222222244449999******************* PP GGGEEEEE--GGG--STTS---E CS UPF0020 102 gdliefsqadaakLrlkegevdv 124 ++ ++++ +a +L+ +++++d+ tr|H9G6W9|H9G6W9_ANOCA 173 SEGLSWVVGNAEELPFNDNNFDI 195 ***************66667775 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 ? 12.9 0.4 2.9e-05 0.066 14 81 .. 430 498 .. 421 524 .. 0.73 Alignments for each domain: == domain 1 score: 12.9 bits; conditional E-value: 2.9e-05 HHHHHHHHHHTT--TTS-EEETT-TTSHHHHHHHHHHTT--TTTT----SS.HHHHHHHHHHHH--EEE CS UPF0020 14 etLAaalvrlagwkpgepllDPmCGsGtilIEaallganiapgllrefvle.lkaeaeeearaelklyg 81 e L +a+ a++ p+++ lD CG+Gti I a +++++ + l ++++e +k +a+ + +++++y+ tr|H9GDC5|H9GDC5_ANOCA 430 EVLYSAIQDWAQVGPETTVLDVCCGTGTIGITLAKKVKKVIGIELCQEAVEdAKVNAQLNDLNNVEFYC 498 6677899999*********************99999999995544333333344444444445555555 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (197 nodes) Target sequences: 18528 (9796552 residues searched) Passed MSV filter: 479 (0.0258528); expected 370.6 (0.02) Passed bias filter: 367 (0.0198079); expected 370.6 (0.02) Passed Vit filter: 30 (0.00161917); expected 18.5 (0.001) Passed Fwd filter: 8 (0.000431779); expected 0.2 (1e-05) Initial search space (Z): 18528 [actual number of targets] Domain search space (domZ): 8 [number of targets reported over threshold] # CPU time: 0.18u 0.01s 00:00:00.19 Elapsed: 00:00:00.07 # Mc/sec: 27570.30 // [ok]