# 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/Gcd10p.hmm.txt # target sequence database: proteomes/Trichoplax_adhaerens.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Gcd10p [M=300] Accession: PF04189.12 Description: Gcd10p 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 7e-40 136.4 5.3 4.1e-23 81.3 1.6 2.5 2 tr|B3S6F9|B3S6F9_TRIAD Uncharacterized protein OS=Trichoplax Domain annotation for each sequence (and alignments): >> tr|B3S6F9|B3S6F9_TRIAD Uncharacterized protein OS=Trichoplax adhaerens OX=10228 GN=TRIADDRAFT_59791 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 81.3 1.6 3.6e-27 4.1e-23 2 157 .. 7 145 .. 6 146 .. 0.91 2 ! 59.3 0.1 1.9e-20 2.1e-16 184 289 .. 143 238 .. 142 253 .. 0.82 Alignments for each domain: == domain 1 score: 81.3 bits; conditional E-value: 3.6e-27 Gcd10p 2 IkpnqhvllkLpsenlkivqvkpnttisLgKfgsfplnliigrpygltfeildkreeeeksrlrvvpaaeleaeslaeeeaeeeeee 88 I ++ v+++ +++k ++++++ + +++ f+l++++g y++tf++ +++ +l+ e++++++++ + tr|B3S6F9|B3S6F9_TRIAD 7 IDYGDKVIIE-KGDKFKYITIRKGRKACIDRLN-FYLDNAVGCYYNTTFDVRNDHLY------------KLHYEEISSSTSKP--AV 77 788999****.777*******************.******************99943............56777777776554..34 PP Gcd10p 89 eeardnreiiddgarQkLtkeeIeeLKkegasagkeiIaklleshtafdqKTaFSqeKYlkrKkkKYlk 157 e+ dn + dd ++QkL++++Ie+ Kk+g+ g ++++kl++++++fd KT+FS+ KY+k+K k +l+ tr|B3S6F9|B3S6F9_TRIAD 78 IEGIDNSTLADDPSSQKLSRSQIEDYKKSGVI-GMDLVEKLVKNSATFDAKTEFSKVKYIKKKLKNHLR 145 469*****************************.********************************9987 PP == domain 2 score: 59.3 bits; conditional E-value: 1.9e-20 Gcd10p 184 elreetlglllslanvhaggryLvvDdtgGLlvaalaeRmgifessakegtitliheneqpnls.llkyfnydaaepehplkkhlkt 269 +lr +tl++ll++a+v+a+++ ++ ++++GL++aa++eRmg+ eg i+++h+++ p + + + n+ p+ +++ +++ tr|B3S6F9|B3S6F9_TRIAD 143 HLRYDTLAQLLTMADVKAKSKLVIAESCQGLVTAAVMERMGG------EGVIIQVHAGDVPTVGySITNCNF----PQAYMD-MILY 218 79****************************************......************999846899999....776666.7888 PP Gcd10p 270 lswlqllepeedetyeeepe 289 ls+ ++ ++e+d+ + + + tr|B3S6F9|B3S6F9_TRIAD 219 LSLRKVNSQESDNMTGLSSK 238 88888887777765543333 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (300 nodes) Target sequences: 11518 (5222075 residues searched) Passed MSV filter: 875 (0.0759681); expected 230.4 (0.02) Passed bias filter: 239 (0.0207501); expected 230.4 (0.02) Passed Vit filter: 23 (0.00199687); expected 11.5 (0.001) Passed Fwd filter: 1 (8.68206e-05); expected 0.1 (1e-05) Initial search space (Z): 11518 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.14u 0.00s 00:00:00.14 Elapsed: 00:00:00.06 # Mc/sec: 26110.38 // [ok]