# hmmscan :: search sequence(s) against a profile 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 sequence file: AA/23726_FusoPortal_Gene_627.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_627 [L=721] Description: # 705211 # 707373 # 1 # ID=1_627;partial=00;start_type=ATG;rbs_motif=3Base/5BMM;rbs_spacer=13-15bp;gc_cont=0.256 Scores for complete sequence (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Model Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5e-75 251.5 1.4 4.1e-71 238.7 0.0 3.4 3 Peptidase_U32 Peptidase family U32 2.5e-21 76.2 4.3 2.5e-21 76.2 4.3 2.8 3 DUF3656 Collagenase ------ inclusion threshold ------ 0.46 10.6 8.8 0.075 13.2 1.5 2.3 2 DHquinase_I Type I 3-dehydroquinase Domain annotation for each model (and alignments): >> Peptidase_U32 Peptidase family U32 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 238.7 0.0 7.3e-75 4.1e-71 4 231 .. 75 280 .. 72 282 .. 0.94 2 ! 6.3 0.0 0.00078 4.4 36 79 .. 567 611 .. 558 619 .. 0.89 3 ! 3.6 0.1 0.0052 29 126 188 .. 638 699 .. 626 712 .. 0.81 Alignments for each domain: == domain 1 score: 238.7 bits; conditional E-value: 7.3e-75 Peptidase_U32 4 lrkyleklaelgvdavivqDpgvlelareeapelelhlStqanitnseavrflkelGasrvvlaRElsleeikaikekvpdvel 87 l+ +l++l+++g da+ivqD+g ++++++++p++e+h Stq++++n+ ++++lkelG++r+vl+REls+eeik+i+e++ d+el 23726_FusoPortal_Gene_627 75 LYPTLKDLYNYGLDAIIVQDIGYADYLHKNFPSIEIHGSTQMTVANYYEINYLKELGFKRIVLPRELSFEEIKEIREHT-DIEL 157 67899**************************************************************************.**** PP Peptidase_U32 88 EvfvhGalcvaysGrcllssylsgrsanrgecaqpcrldyrlieekekgeeepieedekgtyllspkd.lclieelpelleaGv 170 Evfv G+lc+++sG+c++ss+++grs+nrg+caqpcr++y+ ++ y+lspkd l +++e+++l e+G+ 23726_FusoPortal_Gene_627 158 EVFVSGSLCISFSGNCYMSSFIGGRSGNRGMCAQPCRKEYKTSCGEK----S---------YFLSPKDqLYGLDEIKKLQEIGI 228 *****************************************444433....3.........9********************** PP Peptidase_U32 171 dslKIeGRmkspeyvatvvrvYreaidallagreevkeelkeeleklanrgltsgfllgrk 231 +s+K+eGRmk+ +yv+++v+++r +i+ + +++ + kl+nrg+++g+++++ 23726_FusoPortal_Gene_627 229 ESIKVEGRMKDVSYVYETVSYFRSLINGIDKEE---------NTYKLFNRGYSKGYFYNNG 280 ****************************86555.........77899**********9975 PP == domain 2 score: 6.3 bits; conditional E-value: 0.00078 Peptidase_U32 36 elelhlStqanitnseavrflkel.GasrvvlaRElsleeikaik 79 + + l+ +ni+ns ++++l++ + +++ E+++ +i+++ 23726_FusoPortal_Gene_627 567 NSTVMLNWNMNIVNSYTINVLEKIeKLESFIISPEINFSKIRELG 611 667788899***************9*****************985 PP == domain 3 score: 3.6 bits; conditional E-value: 0.0052 Peptidase_U32 126 dyrlieekekgeeepieedekgtyllspkdlclieelpelleaGvdslKIeGRmkspeyvatv 188 ++++i +ke+++ ++i+ +e gt ++ k l++i+ ++++ + +vd + +e + + e +++v 23726_FusoPortal_Gene_627 638 KNEVITNKENDKFNIIK-NEYGTEIFLDKPLNIINIMKDIKKLNVDIVVLEFTTETIEDIKKV 699 66788888877777665.8***************************98877666665555555 PP >> DUF3656 Collagenase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 76.2 4.3 4.5e-25 2.5e-21 1 106 [. 352 474 .. 352 475 .. 0.93 2 ? -0.4 0.6 0.3 1.7e+03 11 64 .. 492 555 .. 478 562 .. 0.57 3 ? -1.4 0.1 0.63 3.5e+03 22 67 .. 625 679 .. 602 694 .. 0.59 Alignments for each domain: == domain 1 score: 76.2 bits; conditional E-value: 4.5e-25 DUF3656 1 vyrtsdkklteelekslekseerkipvdlelaakegeplkltvtde.................leeAkkrplteenirkqLskl 67 ++r+++k+l++e++k++++ +++k+++++ +++k +e+l l+ e ++A+kr+++ee+i+++L++ 23726_FusoPortal_Gene_627 352 IFRNYNKRLNDEISKKIKS-TDKKLEINFDFTGKLDENLILKTYLEdengnrilnleeisetlTQKAQKRAISEEDIKEKLTEI 434 79************99998.899**************9988887667799999999999************************* PP DUF3656 68 GgtpFeleeieidld.glfiPvseLnelRRealekLeear 106 G++ F++++i+id+d ++fiP+seL++l+R+a+ek+ e++ 23726_FusoPortal_Gene_627 435 GDSEFTVKDIKIDIDeNIFIPLSELKNLKRNAVEKFREKI 474 **************8566*******************997 PP == domain 2 score: -0.4 bits; conditional E-value: 0.3 DUF3656 11 eelekslekseerkipvdlelaakegeplkltvtde..........leeAkkrplteenirkqL 64 e ++ ++ek+e + +++++ ++ +e++++ +++de Ak++ l++ n+ ++L 23726_FusoPortal_Gene_627 492 EYFKLEIEKDEPKDLEIKVIVSNEEQKNFLENIKDEynikkvyyrtYDIAKQSMLSQHNLDNKL 555 2444445555555566666665555555544444446666666666666666666666666665 PP == domain 3 score: -1.4 bits; conditional E-value: 0.63 DUF3656 22 erkipvdlelaakegeplkltvtde.........leeAkkrplteenirkqLskl 67 + + +d+ +a++++e ++ + +d+ e +pl+ ni k ++kl 23726_FusoPortal_Gene_627 625 LKGMTIDVDIADNKNEVITNKENDKfniikneygTEIFLDKPLNIINIMKDIKKL 679 3335555555555554332222222444444567777777777777777777777 PP >> DHquinase_I Type I 3-dehydroquinase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 1.2 0.9 0.061 3.4e+02 92 169 .. 481 560 .. 414 579 .. 0.70 2 ? 13.2 1.5 1.3e-05 0.075 87 169 .. 624 703 .. 582 713 .. 0.83 Alignments for each domain: == domain 1 score: 1.2 bits; conditional E-value: 0.061 EEEEEGG..HHG..HHHHHHHHHHHTT-EEEEEEE...ESS----HH.HHHHHHHHHHHTT-SEEEEEE--SSHHHHHHHHHHH CS DHquinase_I 92 yvDvElteapee..skkrlieekreggtkvIvSyh...dfegtpsae.elkallekmqalgaDivKiavmaksieDvlrllslt 169 +D El e ++e + + ++e ++ ++kvIvS+ +f + ++e ++k+++ + Di+K + + D+ +l+ 23726_FusoPortal_Gene_627 481 DLDSELKENNQEyfKLEIEKDEPKDLEIKVIVSNEeqkNFLENIKDEyNIKKVYYRT----YDIAKQSMLSQHNLDNKLASNLY 560 4677775434446433333455778999*****995555555555556666666655....48888888777777776666666 PP == domain 2 score: 13.2 bits; conditional E-value: 1.3e-05 HCT-SEEEEEGG..HHGHHHHHHHHHHHTT-EEEEEEEESS----HH.HHHHHHHHHHHTT-SEEEEEE--SSHHHHHHHHHHH CS DHquinase_I 87 elgadyvDvElteapeeskkrlieekreggtkvIvSyhdfegtpsae.elkallekmqalgaDivKiavmaksieDvlrllslt 169 +l+ +Dv+++ + +k+++i +k++ + +I +e ++ +++++++ +++l+ Div + + ++ieD ++l+ 23726_FusoPortal_Gene_627 624 KLKGMTIDVDIA---D-NKNEVITNKENDKFNIIKNEYGTEIFLDKPlNIINIMKDIKKLNVDIVVLEFTTETIEDIKKVLKQL 703 566688*****4...3.5679*****99999999999988888888889*******************************9876 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (721 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1347 (0.0806008); expected 334.2 (0.02) Passed bias filter: 372 (0.0222595); expected 334.2 (0.02) Passed Vit filter: 35 (0.0020943); expected 16.7 (0.001) Passed Fwd filter: 3 (0.000179512); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.32u 0.23s 00:00:00.55 Elapsed: 00:00:00.30 # Mc/sec: 6986.57 // [ok]