# 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/PCIF1_WW.hmm.txt # target sequence database: proteomes/Homo_sapiens.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: PCIF1_WW [M=172] Accession: PF12237.8 Description: Phosphorylated CTD interacting factor 1 WW domain 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.8e-69 234.4 0.1 4.2e-69 233.8 0.1 1.3 1 sp|Q9H4Z3|PCIF1_HUMAN Phosphorylated CTD-interacting factor ------ inclusion threshold ------ 0.055 15.3 0.0 0.12 14.2 0.0 1.5 1 sp|Q8WV41|SNX33_HUMAN Sorting nexin-33 OS=Homo sapiens GN=SN Domain annotation for each sequence (and alignments): >> sp|Q9H4Z3|PCIF1_HUMAN Phosphorylated CTD-interacting factor 1 OS=Homo sapiens GN=PCIF1 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 233.8 0.1 9.4e-74 4.2e-69 1 172 [] 446 621 .. 446 621 .. 0.99 Alignments for each domain: == domain 1 score: 233.8 bits; conditional E-value: 9.4e-74 PCIF1_WW 1 kLkeLyesseeddekleaflsrvfalllrYdallg.....eaglqaalpeevfevLkkefgvsfecfAsPlncyfeqycsafpdtDay 83 kL Ly++s+ dd+ +e fl rv++ll+rY++++g ++glq++lp +vfe+L++ fgvsfecfAsPlncyf+qycsafpdtD+y sp|Q9H4Z3|PCIF1_HUMAN 446 KLWLLYRYSCIDDSAFERFLPRVWCLLRRYQMMFGvglyeGTGLQGSLPVHVFEALHRLFGVSFECFASPLNCYFRQYCSAFPDTDGY 533 7899***********************************99*********************************************** PP PCIF1_WW 84 FGsvGsffdfkpasgsfeanPPfveevleraaehieklLeaaeeesepLsfvvvvpewre..tealkkleesrykrkevvieakehey 169 FGs+G+++df+p sgsfeanPPf+ee+++++++h+e+lLe++ epLsf+v++pewre t al+++e+sr+kr++++++a ehey sp|Q9H4Z3|PCIF1_HUMAN 534 FGSRGPCLDFAPLSGSFEANPPFCEELMDAMVSHFERLLESS---PEPLSFIVFIPEWREppTPALTRMEQSRFKRHQLILPAFEHEY 618 ******************************************...***************999************************* PP PCIF1_WW 170 reg 172 r+g sp|Q9H4Z3|PCIF1_HUMAN 619 RSG 621 986 PP >> sp|Q8WV41|SNX33_HUMAN Sorting nexin-33 OS=Homo sapiens GN=SNX33 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.2 0.0 2.7e-06 0.12 92 123 .. 420 451 .. 408 460 .. 0.85 Alignments for each domain: == domain 1 score: 14.2 bits; conditional E-value: 2.7e-06 PCIF1_WW 92 dfkpasgsfeanPPfveevleraaehieklLe 123 +f+++s sf++ PPf++e l++a++h ++ e sp|Q8WV41|SNX33_HUMAN 420 AFQAISHSFQMDPPFCSEALNSAISHTGRTYE 451 58899*******************99877655 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (172 nodes) Target sequences: 88473 (35070517 residues searched) Passed MSV filter: 2044 (0.0231031); expected 1769.5 (0.02) Passed bias filter: 1597 (0.0180507); expected 1769.5 (0.02) Passed Vit filter: 96 (0.00108508); expected 88.5 (0.001) Passed Fwd filter: 2 (2.26058e-05); expected 0.9 (1e-05) Initial search space (Z): 88473 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.45u 0.02s 00:00:00.47 Elapsed: 00:00:00.20 # Mc/sec: 30160.64 // [ok]