# 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/LCM.hmm.txt # target sequence database: proteomes/Monosiga_brevicollis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: LCM [M=188] Accession: PF04072.14 Description: Leucine carboxyl 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3e-25 88.3 0.0 8.8e-25 86.8 0.0 1.7 1 tr|A9V7T0|A9V7T0_MONBE Predicted protein OS=Monosiga brevico 7e-25 87.1 0.0 1.5e-24 86.0 0.0 1.5 1 tr|A9UPK8|A9UPK8_MONBE Leucine carboxyl methyltransferase 1 Domain annotation for each sequence (and alignments): >> tr|A9V7T0|A9V7T0_MONBE Predicted protein OS=Monosiga brevicollis OX=81824 GN=10965 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 86.8 0.0 1.9e-28 8.8e-25 50 186 .. 59 271 .. 20 273 .. 0.87 Alignments for each domain: == domain 1 score: 86.8 bits; conditional E-value: 1.9e-28 CHCCCCT--HHHHHHHHHHHHHHHHHHHHHH.HCTT.ESEEEEET-TT--HHHHHHHTT....TEEEEEEE-HHHHHHHHHHHHHSH CS LCM 50 eeekkekkapelnrgiavRtrffddllkaal.saag.irqvViLGAGlDtrayRlpwla....glkvfEvDqpevlefKretlkelg 130 +++ +++p +nrg+++R++++d+++ a+l ++a + q+ LGAG+D++++Rl+ g+++fEvD+p ++++K+++++++ tr|A9V7T0|A9V7T0_MONBE 59 MVQRVTRRSPIINRGYYIRAKAVDHIMSAFLqQHADaAPQILSLGAGFDATFFRLKAAGllppGSVYFEVDFPTLIQRKTALIHQSS 145 245667889**********************65444799******************99999***********************99 PP HHHH.....................................................................EECSEEE--TTS.H CS LCM 131 aepk.....................................................................aelryvpvDlrd.d 147 ++ ++ r v+vDl + + tr|A9V7T0|A9V7T0_MONBE 146 RLSAhlpeapitssdalsdlntpaaqealtqlreqhldraakkaaannaaqasdadsasepdlatllghallyEDYRAVGVDLCKvK 232 99999******************************************************************9999************ PP HHHHHHHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHHH CS LCM 148 dwpeaLaaagfdpekptawvaeGvlvYLtpeasdallea 186 +++ +L+aag+ p++pt++++e+vl Y p+ + a++++ tr|A9V7T0|A9V7T0_MONBE 233 QLRLRLQAAGLLPDRPTLVLSECVLTYVGPKPAFAVVKW 271 ****************************99988888776 PP >> tr|A9UPK8|A9UPK8_MONBE Leucine carboxyl methyltransferase 1 OS=Monosiga brevicollis OX=81824 GN=22106 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 86.0 0.0 3.4e-28 1.5e-24 53 187 .. 86 219 .. 61 220 .. 0.87 Alignments for each domain: == domain 1 score: 86.0 bits; conditional E-value: 3.4e-28 CCCT--HHHHHHHHHHHHHHHHHHHHHHHCTT.ESEEEEET-TT--HHHHHHHTT..TEEEEEEE-HHHHHHHHHHHHHSHHHHHEE CS LCM 53 kkekkapelnrgiavRtrffddllkaalsaag.irqvViLGAGlDtrayRlpwla..glkvfEvDqpevlefKretlkelgaepkae 136 k++ +p ++rg++vR+ +++ ++++++++g +q+ LGAG+Dt+++Rl + +++EvD+p+v+e Kr++ +ka tr|A9UPK8|A9UPK8_MONBE 86 KTHPCSPLIHRGYFVRVQAVERTVRSFIAQHGaQAQILSLGAGFDTLFFRLMNDNvtLRRYVEVDFPAVVEKKRQICS-----SKAP 167 34556899********************5577688*****************99997789**************9985.....6677 PP CSEEE--TTS.HHHHHHHHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHHHH CS LCM 137 lryvpvDlrd.ddwpeaLaaagfdpekptawvaeGvlvYLtpeasdalleai 187 lr v Dl + +++++ L + +dp p ++++e+vlvY+ pea+++l++++ tr|A9UPK8|A9UPK8_MONBE 168 LRSVLSDLASpSEVASLLDRGRLDPGLPLLVISECVLVYMAPEAVANLIQWT 219 899*******99*************************************996 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (188 nodes) Target sequences: 9188 (5511963 residues searched) Passed MSV filter: 448 (0.0487593); expected 183.8 (0.02) Passed bias filter: 209 (0.0227471); expected 183.8 (0.02) Passed Vit filter: 17 (0.00185024); 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.10u 0.00s 00:00:00.10 Elapsed: 00:00:00.05 # Mc/sec: 20724.98 // [ok]