# 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/Ciona_savignyi.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5.8e-38 129.6 0.0 9.3e-38 128.9 0.0 1.3 1 tr|H2ZN24|H2ZN24_CIOSA Uncharacterized protein OS=Ciona savi 1.8e-07 30.4 0.2 2.7e-07 29.8 0.2 1.4 1 tr|H2ZQD2|H2ZQD2_CIOSA Uncharacterized protein OS=Ciona savi 0.00034 19.7 0.2 0.0022 17.1 0.1 2.5 2 tr|H2Y5Z4|H2Y5Z4_CIOSA tRNA (guanine(37)-N1)-methyltransfera 0.0036 16.4 0.1 0.0074 15.4 0.0 1.5 2 tr|H2Z3B7|H2Z3B7_CIOSA Uncharacterized protein OS=Ciona savi ------ inclusion threshold ------ 0.02 14.0 0.0 0.022 13.9 0.0 1.1 1 tr|H2YEI9|H2YEI9_CIOSA Uncharacterized protein OS=Ciona savi 0.024 13.7 0.0 0.035 13.2 0.0 1.2 1 tr|H2YG63|H2YG63_CIOSA Uncharacterized protein OS=Ciona savi Domain annotation for each sequence (and alignments): >> tr|H2ZN24|H2ZN24_CIOSA Uncharacterized protein OS=Ciona savignyi OX=51511 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 128.9 0.0 4.8e-41 9.3e-38 2 164 .. 100 251 .. 99 252 .. 0.94 Alignments for each domain: == domain 1 score: 128.9 bits; conditional E-value: 4.8e-41 Methyltransf_15 2 evvlDlfgGaGGnviqfanvfekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklkleeikvdlvflSpPwGGpsYkk 88 +++D+f+G+GGn+iqfa + e+viaidi++ +l+ a hNa +Ygv+d++++++g +f+l +l k d++flSpPwGGp+Y + tr|H2ZN24|H2ZN24_CIOSA 100 DIIVDAFCGSGGNSIQFAFTCERVIAIDIDPVKLENARHNAAIYGVEDRIEFICGSFFDLAPTL-----KADVIFLSPPWGGPEYIH 181 689********************************************************98775.....45899************* PP Methyltransf_15 89 kekfdlekkleplklekllkesvklskniilfLPrnsdlkqlaaltaevlenaskvkvvklkenglvkallaylge 164 +++++ + ++ ek ++ + ++s+n+ +fLP+ns++++l +l +++ ++v++ + n+ +++l+ay+g tr|H2ZN24|H2ZN24_CIOSA 182 CDTYSI-AAMGEF-GEKAFQLARNISDNVAFFLPKNSNVDELVRL----AGQGKHVELEQNMLNNKISTLTAYYGD 251 ******.888888.58999***********************988....899999*****************9985 PP >> tr|H2ZQD2|H2ZQD2_CIOSA Uncharacterized protein OS=Ciona savignyi OX=51511 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 29.8 0.2 1.4e-10 2.7e-07 3 87 .. 51 137 .. 49 181 .. 0.77 Alignments for each domain: == domain 1 score: 29.8 bits; conditional E-value: 1.4e-10 Methyltransf_15 3 vvlDlfgGaGGnviqfanv.fekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklk.leeikvdlvflSpPwGGpsYk 87 v Dl +G+G +i + ++di+++ l + N+e Y+ ++ v++v+ dv +l + + l++ ++d v++ pP G s k tr|H2ZQD2|H2ZQD2_CIOSA 51 AVCDLGCGCGILSIGSSILgASYCLGVDIDKEALGIFLSNCEDYELESTVDCVHSDVVKLASCRNvLRKKQFDTVIMNPPFGTKSNK 137 5779999999888877666334589*********************************96543321233358*********976644 PP >> tr|H2Y5Z4|H2Y5Z4_CIOSA tRNA (guanine(37)-N1)-methyltransferase OS=Ciona savignyi OX=51511 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.8 0.0 0.72 1.4e+03 60 81 .. 8 29 .. 2 48 .. 0.74 2 ! 17.1 0.1 1.1e-06 0.0022 3 61 .. 256 321 .. 254 416 .. 0.86 Alignments for each domain: == domain 1 score: -1.8 bits; conditional E-value: 0.72 Methyltransf_15 60 ellaklkleeikvdlvflSpPw 81 +l+++l ++++ dl +lSpP+ tr|H2Y5Z4|H2Y5Z4_CIOSA 8 NLFRRLVVNKMDSDLHVLSPPT 29 6778888899999999****96 PP == domain 2 score: 17.1 bits; conditional E-value: 1.1e-06 Methyltransf_15 3 vvlDlfgGaGGnviqfanvfekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfe.........l 61 +v D+f+G G +i +a+ + v+a d n++ k hN ++ + +v++ +d + + tr|H2Y5Z4|H2Y5Z4_CIOSA 256 TVFDVFAGVGPFAIPIAKKGCLVYANDLNPDSYKWLVHNVKLN--KTEVECFNCDGRDfirdnlrnyI 321 589***************************************9..99999999998884444433332 PP >> tr|H2Z3B7|H2Z3B7_CIOSA Uncharacterized protein OS=Ciona savignyi OX=51511 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.4 0.0 3.8e-06 0.0074 39 98 .. 19 79 .. 12 100 .. 0.85 2 ? -3.7 0.0 2.8 5.3e+03 71 80 .. 172 181 .. 154 185 .. 0.66 Alignments for each domain: == domain 1 score: 15.4 bits; conditional E-value: 3.8e-06 Methyltransf_15 39 qhNaevYgvsdkvtlvlgdvfe...llaklkleeikvdlvflSpPwGGpsYkkkekfdlekkl 98 + N e+Y + k + + d++e l++ +k +e ++++++ pPw s ++ +k+d +l tr|H2Z3B7|H2Z3B7_CIOSA 19 KLNSEFYLIPAKSSFLMSDLKEigsLVEACKHKEK-FNCIVMDPPWENKSVRRGKKYDW-LSL 79 5699***************9998878999998887.9********************97.333 PP == domain 2 score: -3.7 bits; conditional E-value: 2.8 Methyltransf_15 71 kvdlvflSpP 80 +lv++S+P tr|H2Z3B7|H2Z3B7_CIOSA 172 IKELVICSVP 181 3567777777 PP >> tr|H2YEI9|H2YEI9_CIOSA Uncharacterized protein OS=Ciona savignyi OX=51511 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.9 0.0 1.1e-05 0.022 78 133 .. 41 99 .. 23 134 .. 0.83 Alignments for each domain: == domain 1 score: 13.9 bits; conditional E-value: 1.1e-05 Methyltransf_15 78 SpPwGGpsYkkk...ekfdlekkleplklekllkesvklskniilfLPrnsdlkqlaal 133 P GG s k ek +e k + ++ k++++++k k+i f P + +l+ +++ tr|H2YEI9|H2YEI9_CIOSA 41 ANPFGGASHAKGivlEKIGVEAKQPNSAIRKCVRVQLKNGKKITAFVPNDGCLNYIEEN 99 6799999999872225556888999*****************************99775 PP >> tr|H2YG63|H2YG63_CIOSA Uncharacterized protein OS=Ciona savignyi OX=51511 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.2 0.0 1.8e-05 0.035 4 39 .. 33 68 .. 31 103 .. 0.89 Alignments for each domain: == domain 1 score: 13.2 bits; conditional E-value: 1.8e-05 Methyltransf_15 4 vlDlfgGaGGnviqfanvfekviaidinlsrlklaq 39 lD+ + GG + ++a f+ vi+id++ +k ++ tr|H2YG63|H2YG63_CIOSA 33 ALDIGCAVGGASFTLASKFDDVIGIDYSHGFVKTCN 68 79*************************998888775 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (165 nodes) Target sequences: 11592 (4231397 residues searched) Passed MSV filter: 391 (0.0337302); expected 231.8 (0.02) Passed bias filter: 277 (0.0238958); expected 231.8 (0.02) Passed Vit filter: 32 (0.00276052); expected 11.6 (0.001) Passed Fwd filter: 6 (0.000517598); expected 0.1 (1e-05) Initial search space (Z): 11592 [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.03 # Mc/sec: 23272.68 // [ok]