# 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/Dictyostelium_discoideum.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 8.4e-87 291.1 0.0 1.1e-86 290.6 0.0 1.2 1 tr|Q54PT4|Q54PT4_DICDI Uncharacterized protein OS=Dictyostel Domain annotation for each sequence (and alignments): >> tr|Q54PT4|Q54PT4_DICDI Uncharacterized protein OS=Dictyostelium discoideum OX=44689 GN=DDB0185958 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 290.6 0.0 9e-91 1.1e-86 3 306 .. 58 352 .. 56 355 .. 0.95 Alignments for each domain: == domain 1 score: 290.6 bits; conditional E-value: 9e-91 XX--TTTHHHHHHT---TT-EEEES--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT.....---EEEEES-GGGHHHHHHHTT CS Methyltransf_5 3 eHvpVLlqevvellepkedgvyiDaTlGaggHsealLeqlselrligidrDksalarakerLkefekrltlvkgnfkelkealaean 89 H+pV+++ev+e+l+p++++vyiDaT+G ggHs+ +Le ++ +i++drD+++ +++k +ef++rl +++++f+++++al+e+n tr|Q54PT4|Q54PT4_DICDI 58 IHKPVFIKEVLEILNPQPNKVYIDATFGMGGHSRGILETCKDCYVIAFDRDPKVFEMTKPLRSEFKDRLITINSTFSSMRSALKEKN 144 5*************************************************************************************9 PP ..TT-EEEEEEE-S--HHHHH-GGG---SSS..--BS---S-TTSS-.-HHHHHHH--HHHHHHHHHHHH--TTHHHHHHHHHHHHH CS Methyltransf_5 90 ..vtkvdgillDLGvSslqldkaerGFsfae..dapLDmRmdqsael.tAaevlnrlseeeLakilkeygeekfakriakaivekrk 171 k+ gil+D GvSs+q+d+a+rGFsf++ ++pLDmRm+q++++ tA++++n++se++L +i++ ygeek++k+ia++i +krk tr|Q54PT4|Q54PT4_DICDI 145 ldHLKISGILYDFGVSSYQIDTADRGFSFKKelEGPLDMRMNQNSSQiTAYDIINTFSEKDLRDIFFYYGEEKHTKNIANEILRKRK 231 655699***********************973369********99988*************************************** PP HS---BHHHHHHHHHHHS-------..-TTHHHHHHHHHHHHTHHHHHHHHHHHHHHHEEEEEEEEEEESSHHHHHHHHHHHHHTT- CS Methyltransf_5 172 kkpiettkelaevikkvvpakekrkkihpatrvFQAlRIavNdEleeleelleaalellakggrlavisFhSledrivkkvfkeask 258 +pi++t++l++++++ +p + k rvF +lR++vNdEl+e++ l++ +ll+++ l isFhS+edrivk++++++++ tr|Q54PT4|Q54PT4_DICDI 232 IEPIKNTSDLVKIVESSIPYPMALK---TISRVFRSLRMFVNDELNEIKIGLTQGEQLLEPNSPLLAISFHSIEDRIVKQYLNKYTT 315 *******************999866...889*******************************************************9 PP ------------------EEEEE.E.E---HHHHHH-GGGTT-EEEEE CS Methyltransf_5 259 keevPaglpvaseelkaklklitkkpilpseeeikeNprsrsakLrvi 306 ee + k++ k pi ps+ee++ Nprs+sa Lr + tr|Q54PT4|Q54PT4_DICDI 316 -----------IEEPTFKFSPTYKSPIFPSKEELEWNPRSKSAILRTA 352 ...........2444568888999**********************87 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (309 nodes) Target sequences: 12739 (6856946 residues searched) Passed MSV filter: 654 (0.0513384); expected 254.8 (0.02) Passed bias filter: 187 (0.0146793); expected 254.8 (0.02) Passed Vit filter: 16 (0.00125599); expected 12.7 (0.001) Passed Fwd filter: 1 (7.84991e-05); expected 0.1 (1e-05) Initial search space (Z): 12739 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.14u 0.01s 00:00:00.15 Elapsed: 00:00:00.06 # Mc/sec: 35313.27 // [ok]