# 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/Dictyostelium_discoideum.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.8e-29 101.5 0.3 4.6e-28 97.9 0.0 2.4 2 tr|Q54NV3|Q54NV3_DICDI Uncharacterized protein OS=Dictyostel 5e-22 78.3 0.0 2.5e-21 76.0 0.0 2.0 1 tr|Q54J74|Q54J74_DICDI Leucine carboxyl methyltransferase 1 Domain annotation for each sequence (and alignments): >> tr|Q54NV3|Q54NV3_DICDI Uncharacterized protein OS=Dictyostelium discoideum OX=44689 GN=DDB0218658 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -0.3 1.0 0.097 6.2e+02 16 47 .. 55 86 .. 37 125 .. 0.57 2 ! 97.9 0.0 7.3e-32 4.6e-28 19 187 .. 195 415 .. 187 416 .. 0.82 Alignments for each domain: == domain 1 score: -0.3 bits; conditional E-value: 0.097 HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHC CS LCM 16 adglidDplaeslvsaaglellerkadkelek 47 +++l++ + ++ ++a+++++++++++++ ++ tr|Q54NV3|Q54NV3_DICDI 55 QQQLLQQQQQKQAAEAEQKRIMQQQQQQQQQQ 86 23333333333333333333322222222222 PP == domain 2 score: 97.9 bits; conditional E-value: 7.3e-32 SSHTCCCCCCCHHHHHHHH...........................CHCCCHCHHCHGC........HCCCCT--HHHHHHHHHHHH CS LCM 19 lidDplaeslvsaaglell...........................erkadkelekeke........eekkekkapelnrgiavRtr 70 i+Dp+a +++ ++ ++l +++++++++++++ + + +k + + i Rt+ tr|Q54NV3|Q54NV3_DICDI 195 CIYDPFAIFFLDTEESKQLlldsfqlfhsvspqqqflnnssnniniNNNNNNNVTYKEFimkshpiqSWLEFQKGGNSINKIGFRTK 281 67777888777777777777777777888888877788777777775555555555553355653322333344445556689**** PP HHHHHHHHHHHCTT...........ESEEEEET-TT--HHHHHHHTTTEEEEEEE-HHHHHHHHHHHHHSHHHHH....EECSEEE- CS LCM 71 ffddllkaalsaag...........irqvViLGAGlDtrayRlpwlaglkvfEvDqpevlefKretlkelgaepk....aelryvpv 142 ++d+++ + ls +g qvV G+GlD+r++Rlp ++++vfE+D p v+++K+++l+++ +k ++ +y++ tr|Q54NV3|Q54NV3_DICDI 282 YIDQVIMDSLSIKGtvahfnasaffYTQVVFCGSGLDSRSLRLPIDKHVSVFEIDLPHVINYKKKVLEKATTIIKptsqCHINYIGG 368 **********87899999999999999*****************************************998777766668999**** PP -TTS..HHHHHHHHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHHHH CS LCM 143 Dlrd..ddwpeaLaaagfdpekptawvaeGvlvYLtpeasdalleai 187 D+ + d+w L ++gf+p+ t+wv+eGvl+YLt+e++ ll i tr|Q54NV3|Q54NV3_DICDI 369 DIVKegDRWIGSLLDCGFNPAMHTLWVLEGVLMYLTEEEIGTLLGSI 415 **99899*********************************9999766 PP >> tr|Q54J74|Q54J74_DICDI Leucine carboxyl methyltransferase 1 OS=Dictyostelium discoideum OX=44689 GN=DDB0219323 PE=3 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 76.0 0.0 3.9e-25 2.5e-21 5 186 .. 51 246 .. 47 248 .. 0.88 Alignments for each domain: == domain 1 score: 76.0 bits; conditional E-value: 3.9e-25 HHHHHHHHHHH..HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHH.HCTT..ES CS LCM 5 aaaralesrra..adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaal.saag..ir 86 ++ +a++++ + + g++ Dp+++++v++ ++++p +nrg++ R+ +++l+ ++ + + + tr|Q54J74|Q54J74_DICDI 51 TNDDAASCKLSavNVGYYSDPFVKYFVKHP----------------------IRRQPLINRGFFSRVECIEQLVSQFFtQYKDinKK 115 666777777644677777777777777766......................4579*********************9554447789 PP EEEEET-TT--HHHHHHHTT....TEEEEEEE-HHHHHHHHHHHHHSHHHHH...........................EECSEEE- CS LCM 87 qvViLGAGlDtrayRlpwla....glkvfEvDqpevlefKretlkelgaepk...........................aelryvpv 142 q++ LG+G+Dt +Rl + + +fEvD+ +v++ K++++++++++++ ++ r+ ++ tr|Q54J74|Q54J74_DICDI 116 QIISLGCGFDTYYFRLMNNKdikkDFIYFEVDYDQVISNKIKIIQNHKELQSmidqewdskydtneklasmvnhqrvssKSYRLGSI 202 ****************99998989999***********************99********************999998778888*** PP -TTSHHHHHHHHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHHH CS LCM 143 DlrdddwpeaLaaagfdpekptawvaeGvlvYLtpeasdallea 186 Dl++ + ++ + +d ++pt++++e+vlvY++ + + ++++ tr|Q54J74|Q54J74_DICDI 203 DLTNLETFKIFDELEIDYNVPTLFLSECVLVYIPTKCGNDVIQW 246 *********************************99887777776 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (188 nodes) Target sequences: 12739 (6856946 residues searched) Passed MSV filter: 734 (0.0576183); expected 254.8 (0.02) Passed bias filter: 341 (0.0267682); expected 254.8 (0.02) Passed Vit filter: 30 (0.00235497); expected 12.7 (0.001) Passed Fwd filter: 2 (0.000156998); expected 0.1 (1e-05) Initial search space (Z): 12739 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.15u 0.01s 00:00:00.16 Elapsed: 00:00:00.06 # Mc/sec: 21485.10 // [ok]