# 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_10.hmm.txt # target sequence database: proteomes/Anolis_carolinensis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Methyltransf_10 [M=299] Accession: PF05971.11 Description: Protein of unknown function (DUF890) 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.8e-130 433.3 0.0 6.9e-130 432.8 0.0 1.2 1 tr|G1KT98|G1KT98_ANOCA U6 small nuclear RNA (adenine-(43)-N( ------ inclusion threshold ------ 0.085 12.2 0.5 0.1 12.0 0.5 1.2 1 tr|H9G6A4|H9G6A4_ANOCA Uncharacterized protein OS=Anolis car Domain annotation for each sequence (and alignments): >> tr|G1KT98|G1KT98_ANOCA U6 small nuclear RNA (adenine-(43)-N(6))-methyltransferase OS=Anolis carolinensis OX=28377 GN # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 432.8 0.0 7.4e-134 6.9e-130 1 299 [] 1 292 [. 1 292 [. 0.98 Alignments for each domain: == domain 1 score: 432.8 bits; conditional E-value: 7.4e-134 XXXXXXXXXXXXXXXXX.XXXXXXXXXXXXXXXXXXXXXXXXXS-TTSHHHHHHHHHHHHHHHH----.---TTS----HHHHHHHH CS Methyltransf_10 1 aalksglhprnrhkdky.dfaalikvvPelkqfvilnpkGrvsidfadPlavkalnkallrefygvsewdiPdgfliPpvPgradyi 86 +al++++h+rnr+kdk+ dfa l++++Pe+kq+v++n++Grvs++f+dP+av+al+++ll+e++g+++ diP ++liP+vP+r++yi tr|G1KT98|G1KT98_ANOCA 1 MALNKSMHARNRYKDKPpDFASLASKYPEFKQHVQINLSGRVSVNFKDPEAVRALTCTLLKEDFGLQI-DIPLERLIPTVPLRLNYI 86 6889************99**************************************************.****************** PP HHHHHHHT--TCG.CS---EEEEES-TTTTHHHHHHHHHH--EEEEEES-HHHHHHHHHHHHHT.T-TTTEEEEE--ST-SSTTTST CS Methyltransf_10 87 hyvadllgkadsd.isaqvraldiGvGanCiyPllGvteyGWrfvgsevddvslasakaivesnsklsdaiklrkqkekkaifdgli 172 h+v+dl++++++ + + + ++diG+Ga+CiyPllG+t++GW+f+++evddv++++ak++ve+n +lsd+ik++k+++k+ ++d+l tr|G1KT98|G1KT98_ANOCA 87 HWVEDLIAQQETAaKDSVIWGIDIGTGASCIYPLLGATLNGWYFIATEVDDVCFSYAKKNVEQN-NLSDLIKVVKVPQKTLLLDALE 172 **********9984555899********************************************.********************** PP T--S.-EEEEEE-----SS------------...---.....---------TTTTHHHHTHHHHHHHHHHHHHHHGGGEEEEEEEES CS Methyltransf_10 173 eene.rydvtlCnPPfhaslaeakagserklknlgkanktralpslnfgGqkaelyceGGevafikklieeskqlakqvlWfttlvs 258 e++e +yd+++CnPPf+a+++eak++++r++ +++ +ps++++G+++e+++eGGe++f+k++i++s ql+k+++W++++++ tr|G1KT98|G1KT98_ANOCA 173 EKSEvIYDFCMCNPPFFANQLEAKGVNSRNP---RRP-----PPSSVNTGGITEIMAEGGELEFVKRIIHDSLQLKKRLRWYSCMLG 251 *******************************...776.....899****************************************** PP SGGGHHHHHHHHHHTT-SEEEEEEEEETTEEEEEEEEES-- CS Methyltransf_10 259 kkenlaplkeelkklGaskvtvvelaqGqklsrfiaWsfld 299 kk++laplk+el++ G++kvt +e++qG+++++++aWsf++ tr|G1KT98|G1KT98_ANOCA 252 KKCSLAPLKKELRIHGVPKVTHTEFCQGRTMRWALAWSFYE 292 ***************************************85 PP >> tr|H9G6A4|H9G6A4_ANOCA Uncharacterized protein 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 ? 12.0 0.5 1.1e-05 0.1 221 290 .. 28 101 .. 4 107 .. 0.72 Alignments for each domain: == domain 1 score: 12.0 bits; conditional E-value: 1.1e-05 --TTTTHHHHTHHHHHHHHHHHHHHHGGGEEEEE....EEESSGGGHHHHHHHHHHTT-SEEEEEEEEETTEEE CS Methyltransf_10 221 GqkaelyceGGevafikklieeskqlakqvlWft....tlvskkenlaplkeelkklGaskvtvvelaqGqkls 290 qk + +++ e + ++ li+ +kq++ + + + +l+ +++ l ++ ++l+++ a+ ++ + ++qG + tr|H9G6A4|H9G6A4_ANOCA 28 QQKVQALAKRAEKEAVRSLIQHTKQIKAENRQLRdellSLIKRAHGLDAIVRQLREQRAQLLREMRYCQGLADK 101 4555566667777777777777777765444332222299******************************6544 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (299 nodes) Target sequences: 18528 (9796552 residues searched) Passed MSV filter: 427 (0.0230462); expected 370.6 (0.02) Passed bias filter: 353 (0.0190522); expected 370.6 (0.02) Passed Vit filter: 22 (0.00118739); expected 18.5 (0.001) Passed Fwd filter: 2 (0.000107945); expected 0.2 (1e-05) Initial search space (Z): 18528 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.19u 0.01s 00:00:00.20 Elapsed: 00:00:00.07 # Mc/sec: 41845.27 // [ok]