# 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/Monosiga_brevicollis.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- ------ inclusion threshold ------ 0.019 13.5 0.0 0.048 12.2 0.0 1.6 2 tr|A9VAX4|A9VAX4_MONBE Predicted protein OS=Monosiga brevico 0.077 11.5 0.0 0.13 10.8 0.0 1.3 1 tr|A9UYM4|A9UYM4_MONBE Predicted protein OS=Monosiga brevico Domain annotation for each sequence (and alignments): >> tr|A9VAX4|A9VAX4_MONBE Predicted protein OS=Monosiga brevicollis OX=81824 GN=34319 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.2 0.0 1e-05 0.048 18 97 .. 164 244 .. 144 248 .. 0.82 2 ? -3.0 0.0 0.44 2e+03 71 97 .. 496 523 .. 479 524 .. 0.72 Alignments for each domain: == domain 1 score: 12.2 bits; conditional E-value: 1e-05 --TT-EEEES--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT.....---...EEEEES-GGGHHHHHHHTTTT-EEEEE CS Methyltransf_5 18 pkedgvyiDaTlGaggHsealLeqlselrligidrDksalarakerLkefekr...ltlvkgnfkelkealaeanvtkvdgil 97 +k++ +D+ G+g ++a+ ++++ + +g+d D+ +++ra++ +++ + ++++ ++ ++++++ + + +k++ il tr|A9VAX4|A9VAX4_MONBE 164 QKPGCRVLDVGCGTGESTRAIALAYPNAKCTGVDPDTASIERARAAATNMGTKgnqMSFFDQTIEAFADEAEGH--GKFNVIL 244 47778889************************************988887765444888888888887776655..5566655 PP == domain 2 score: -3.0 bits; conditional E-value: 0.44 EEEEES-G.GGHHHHHHHTTTT-EEEEE CS Methyltransf_5 71 ltlvkgnf.kelkealaeanvtkvdgil 97 t +++f ++l++ +e +v++vdgi+ tr|A9VAX4|A9VAX4_MONBE 496 TTADQQSFnERLHALAKECGVEDVDGII 523 5555677878889999999999999996 PP >> tr|A9UYM4|A9UYM4_MONBE Predicted protein OS=Monosiga brevicollis OX=81824 GN=25176 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.8 0.0 2.8e-05 0.13 10 86 .. 170 246 .. 167 271 .. 0.86 Alignments for each domain: == domain 1 score: 10.8 bits; conditional E-value: 2.8e-05 HHHHHHT---TT-EEEES--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT.....---EEEEES-GGGHHHHHH CS Methyltransf_5 10 qevvellepkedgvyiDaTlGaggHsealLeqlselrligidrDksalarakerLkefekrltlvkgnfkelkeala 86 ++++++++ +++ +D G+g H +l ++ s+l+l g+d+ +++a++r + +k+++++ g+ ++l + + tr|A9UYM4|A9UYM4_MONBE 170 DAILNAVALPDQAAVLDLCCGQGRHLLELSQRRSDLQLFGLDQSAFLIQTAQARADASQKNVHFTVGTCDALPYEAN 246 567777888899999**************99999*********************9999999999888877766555 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (309 nodes) Target sequences: 9188 (5511963 residues searched) Passed MSV filter: 316 (0.0343927); expected 183.8 (0.02) Passed bias filter: 184 (0.0200261); expected 183.8 (0.02) Passed Vit filter: 19 (0.00206791); expected 9.2 (0.001) Passed Fwd filter: 2 (0.000217675); expected 0.1 (1e-05) Initial search space (Z): 9188 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.11u 0.00s 00:00:00.11 Elapsed: 00:00:00.04 # Mc/sec: 42579.91 // [ok]