# 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/Danio_rerio.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.6e-70 235.2 0.2 6.7e-70 234.7 0.2 1.3 1 tr|A0A0R4IKJ1|A0A0R4IKJ1_DANRE PDX1 C-terminal-inhibiting facto ------ inclusion threshold ------ 0.012 15.8 0.1 0.32 11.1 0.0 2.0 2 tr|Q5EFQ8|Q5EFQ8_DANRE Interleukin-10 OS=Danio rerio OX Domain annotation for each sequence (and alignments): >> tr|A0A0R4IKJ1|A0A0R4IKJ1_DANRE PDX1 C-terminal-inhibiting factor 1 OS=Danio rerio OX=7955 GN=pcif1 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 234.7 0.2 5.3e-74 6.7e-70 1 172 [] 451 626 .. 451 626 .. 0.99 Alignments for each domain: == domain 1 score: 234.7 bits; conditional E-value: 5.3e-74 PCIF1_WW 1 kLkeLyesseeddekleaflsrvfalllrYdallg.....eaglqaalpeevfevLkkefgvsfecfAsPlncyfeqyc 74 kL+ Ly++s+ dd+++e+flsrv++l++rY+ ++g ++glq++lp vfe+L+k+fgv+fecfAsPlncyf+q+c tr|A0A0R4IKJ1|A0A0R4IKJ1_DANRE 451 KLELLYRYSCIDDPRFEKFLSRVWCLIKRYQVMFGsgvneGSGLQGSLPVPVFEALNKQFGVTFECFASPLNCYFKQFC 529 89*********************************99*999************************************** PP PCIF1_WW 75 safpdtDayFGsvGsffdfkpasgsfeanPPfveevleraaehieklLeaaeeesepLsfvvvvpewre..tealkkle 151 safpd+D++FGs+G+f++f+pasgsfeanPPf+ee+++++++h+e+lL ++ sepLsf+++vpewr+ t al+++e tr|A0A0R4IKJ1|A0A0R4IKJ1_DANRE 530 SAFPDIDGFFGSRGPFLSFSPASGSFEANPPFCEELMDAMVTHFEDLLGRS---SEPLSFIIFVPEWRDppTPALTRME 605 ***************************************************...****************99******* PP PCIF1_WW 152 esrykrkevvieakeheyreg 172 +sr++r++++++a eheyr+g tr|A0A0R4IKJ1|A0A0R4IKJ1_DANRE 606 ASRFRRHQMTVPAFEHEYRSG 626 ******************986 PP >> tr|Q5EFQ8|Q5EFQ8_DANRE Interleukin-10 OS=Danio rerio OX=7955 GN=il10 PE=2 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.1 0.0 2.5e-05 0.32 53 132 .. 23 104 .. 16 108 .. 0.82 2 ? 2.9 0.0 0.008 1e+02 86 120 .. 113 152 .. 106 169 .. 0.79 Alignments for each domain: == domain 1 score: 11.1 bits; conditional E-value: 2.5e-05 PCIF1_WW 53 kefgvsfec..fAsPlncyfeqycsafpdtDayFGsvGsffdfkpasgsfeanPPfveevleraaeh.ieklLeaaeeesepL 132 ++++++++c f + ++++ sa++++ +++ s ++ + ++ +++n P+ +v++++++ +e++L +a +++pL tr|Q5EFQ8|Q5EFQ8_DANRE 23 RRVECKTDCcsFVEGFPLRLRELRSAYKEIQKFYESNDDLEPLLNEDIKHNINSPYGCHVMNEILHFyLETILPTA-LQKNPL 104 6788899983345789999********************************************9875145556555.555555 PP == domain 2 score: 2.9 bits; conditional E-value: 0.008 PCIF1_WW 86 svGsff.....dfkpasgsfeanPPfveevleraaehiek 120 s+G++f d+ + + f ++ Pf + l++ +e++++ tr|Q5EFQ8|Q5EFQ8_DANRE 113 SIGNIFqelkrDMVKCKRYFSCQNPFEVNSLKNSYEKMKE 152 7888886443344666799************999999876 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (172 nodes) Target sequences: 25687 (13798380 residues searched) Passed MSV filter: 781 (0.0304045); expected 513.7 (0.02) Passed bias filter: 533 (0.0207498); expected 513.7 (0.02) Passed Vit filter: 28 (0.00109005); expected 25.7 (0.001) Passed Fwd filter: 2 (7.78604e-05); expected 0.3 (1e-05) Initial search space (Z): 25687 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.19u 0.01s 00:00:00.20 Elapsed: 00:00:00.08 # Mc/sec: 29666.52 // [ok]