# 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_5.hmm.txt # target sequence database: proteomes/Danio_rerio.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Methyltransf_5 [M=309] Accession: PF01795.18 Description: MraW methylase family 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.3e-86 291.5 0.0 1.1e-85 288.4 0.0 1.9 1 tr|A0A1D5NT05|A0A1D5NT05_DANRE Methyltransferase-like 15 OS=Dan ------ inclusion threshold ------ 0.15 12.0 0.0 0.22 11.5 0.0 1.1 1 tr|F1R9J8|F1R9J8_DANRE Zgc:194242 OS=Danio rerio OX=795 Domain annotation for each sequence (and alignments): >> tr|A0A1D5NT05|A0A1D5NT05_DANRE Methyltransferase-like 15 OS=Danio rerio OX=7955 GN=mettl15 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 288.4 0.0 8.3e-90 1.1e-85 3 306 .. 66 399 .. 64 402 .] 0.88 Alignments for each domain: == domain 1 score: 288.4 bits; conditional E-value: 8.3e-90 XX--TTTHHHHHHT---TT-EEEES--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT.....---EEEEES-GGGH CS Methyltransf_5 3 eHvpVLlqevvellepkedgvyiDaTlGaggHsealLeqlselrligidrDksalarakerLkefekrltlvkgnfkel 81 H+pV+++ev+++l++k++++ +D T+GaggHs+a+L+ ++ ++++ drD++a +a++ +e++ +++ v g+f+el tr|A0A1D5NT05|A0A1D5NT05_DANRE 66 FHTPVMVKEVLQFLNIKPGQIVLDLTFGAGGHSKAILQSVPGVTVVAADRDPTAFRMAQQLAEEYQGKVKPVLGRFSEL 144 6****************************************************************************** PP HHHHHHTT..TT-EEEEEEE-S--HHHHH-GGG---SSS--BS---S-TTSS-...-HHHHHHH--HHHHHHHHHHHH- CS Methyltransf_5 82 kealaean..vtkvdgillDLGvSslqldkaerGFsfaedapLDmRmdqsael...tAaevlnrlseeeLakilkeyge 155 ++ l + ++vd+ llD G Ss+q+d+aerGFs +d+pLDmRmd ++ Aa+v+n l++++La++l yge tr|A0A1D5NT05|A0A1D5NT05_DANRE 145 NNLLPALGlgPGTVDAALLDAGCSSMQMDSAERGFSLSKDGPLDMRMDGDRYPdmpCAADVVNALDQQALASVLAAYGE 223 *****9996567************************************987533369********************** PP -TTHHHHHHHHHHHHHHS---BHHHHHHHHHHHS-------.........-TTHHHHHHHHHHHHTHHHHHHHHHHHHH CS Methyltransf_5 156 ekfakriakaivekrkkkpiettkelaevikkvvpakekrk.......kihpatrvFQAlRIavNdEleeleelleaal 227 e+++++ia+ai+++r+ pi +t +la++++ ++pa++ + h at++FQAlRI+vNdEl+el+ l+ a tr|A0A1D5NT05|A0A1D5NT05_DANRE 224 ERHTRKIAAAIAQARSVYPIGRTLQLASIVAGAFPASALYArrdrlqrPSHVATKTFQALRIFVNDELNELHAGLRVAQ 302 ***********************************998776667777789***************************** PP HHEEEEEEEEEEESSHHHHHHHHHHHHHTT------------------...................-EEEEE..E.E- CS Methyltransf_5 228 ellakggrlavisFhSledrivkkvfkeaskkeevPaglpvaseelka...................klklit.kkpil 286 +ll+++grl+v++FhSledr++k++++ + P++ +++++ + ++ ++ kk+i tr|A0A1D5NT05|A0A1D5NT05_DANRE 303 TLLRPKGRLCVLTFHSLEDRLTKRFLRGEDL--AAPPRRSLRQHARARpagqeeefeeeeegdeggnSVHWVSlKKVIK 379 **************************97655..2332222222221112233334444555777777666666599*** PP --HHHHHH-GGGTT-EEEEE CS Methyltransf_5 287 pseeeikeNprsrsakLrvi 306 pseee++eNpr+rsakLr + tr|A0A1D5NT05|A0A1D5NT05_DANRE 380 PSEEEVQENPRARSAKLRTA 399 ******************87 PP >> tr|F1R9J8|F1R9J8_DANRE Zgc:194242 OS=Danio rerio OX=7955 GN=zgc:194242 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 1.7e-05 0.22 4 74 .. 29 103 .. 27 124 .. 0.81 Alignments for each domain: == domain 1 score: 11.5 bits; conditional E-value: 1.7e-05 X--TTTHHHHHHT---TT-EEEES--T..TSHHHHHHHTT--TTEEEEEES-HHHHHHHTT.....---..EEEE CS Methyltransf_5 4 HvpVLlqevvellepkedgvyiDaTlG..aggHsealLeqlselrligidrDksalarakerLkefekr..ltlv 74 H VL +++v++ +++++ + + G g H L +e +l g+d+ + + ++a +r+ke + r ++l+ tr|F1R9J8|F1R9J8_DANRE 29 HNQVLEENAVKFCNIEPNNTVLELCHGpgFGLHYASQLLTGPEGKLFGVDHSQYMHEMASKRMKEQISRgkVELY 103 999************99999999887522678998888899************************9988444443 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (309 nodes) Target sequences: 25687 (13798380 residues searched) Passed MSV filter: 985 (0.0383462); expected 513.7 (0.02) Passed bias filter: 471 (0.0183361); expected 513.7 (0.02) Passed Vit filter: 36 (0.00140149); expected 25.7 (0.001) Passed Fwd filter: 2 (7.78604e-05); expected 0.3 (1e-05) Initial search space (Z): 25687 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.29u 0.01s 00:00:00.30 Elapsed: 00:00:00.10 # Mc/sec: 42636.99 // [ok]