# 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_318.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_318 [L=202] Description: # 344032 # 344637 # -1 # ID=1_318;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.279 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1e-48 165.0 17.0 4e-48 163.1 17.0 1.8 1 GrpE GrpE ------ inclusion threshold ------ 0.12 11.8 15.8 0.071 12.5 14.4 1.4 1 Golgin_A5 Golgin subfamily A member 5 0.82 10.0 9.2 5.1 7.4 8.6 2.6 1 Syntaxin_2 Syntaxin-like protein 1.3 8.2 4.1 1.7 7.8 4.1 1.2 1 DUF4407 Domain of unknown function (DUF4407) 8.5 5.8 12.3 17 4.9 12.3 1.5 1 Bap31 B-cell receptor-associated protein 31-like 9.1 6.3 10.7 13 5.7 10.7 1.2 1 DUF262 Protein of unknown function DUF262 Domain annotation for each model (and alignments): >> GrpE GrpE # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 163.1 17.0 1.4e-51 4e-48 3 166 .] 38 200 .. 2 200 .. 0.93 Alignments for each domain: == domain 1 score: 163.1 bits; conditional E-value: 1.4e-51 ....CTTTHHHHCCCCCCTTCHHHCTTSCCCCCCTTTTTTTTHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHC..HHH CS GrpE 3 neeeeeeeedeeeeleeeleelekevkelkekllrllaefenlrkrtekekeeakkfaiekfakdLLpvlDnleralksakene 86 +++++e++ + + +ee+e+l++e++e+k+++lr++a+f+n++kr eke +e+kkfa ek+++++L lDn+era++s++e++ 23726_FusoPortal_Gene_318 38 EHKNDEHACCGKHNHKEEIEKLKAEIEEWKNSFLRKQADFQNFTKRKEKEVDELKKFASEKIITQFLGSLDNFERAIESSSESK 121 455566666677788888****************************************************************** PP -HHHHHHHHHHHHHHHHHHHHHTT-EE---TTSB--TTTEEEEEEECSSC...TBEEEEEE--EEETTEEEE-EEEEEEE CS GrpE 87 eekslleGvemtlkqllkiLekhgvekiepkgepFDpnlheavstvesdekeegtvvevvqkGYklkdrvlRpAmVkVsk 166 + +sll+Gvem+ ++l++i++++ ve+i+++ +F+p++h+av ++ s++k+e+++v+v+qkGY++k++v+RpAmV+V+k 23726_FusoPortal_Gene_318 122 DFDSLLQGVEMIVRNLKDIMSSEDVEEIPTE-GAFNPEYHHAVGVEASEDKKEDEIVKVLQKGYMMKGKVIRPAMVIVCK 200 ******************************5.57********************************************86 PP >> Golgin_A5 Golgin subfamily A member 5 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.5 14.4 2.6e-05 0.071 68 158 .. 6 98 .. 1 120 [. 0.84 Alignments for each domain: == domain 1 score: 12.5 bits; conditional E-value: 2.6e-05 Golgin_A5 68 qqieqlrtel.qdleaqqqaeaeslreqlqeleeqleee.rkakeeveaelerlkeelryleeelvrskatlqsrikdrekeie 149 ++ e l+ e+ ++++++++ e e+ +e+ +e e++ +e k++ ++e+e+lk e++ ++ r++a++q+ k++eke++ 23726_FusoPortal_Gene_318 6 IKEEVLKEEInKEVNEKKKCECEEGKEEAHEHEHKNDEHaCCGKHNHKEEIEKLKAEIEEWKNSFLRKQADFQNFTKRKEKEVD 89 5557799999779***********99999888777665515679999************************************* PP Golgin_A5 150 klrnqlttk 158 +l++ + k 23726_FusoPortal_Gene_318 90 ELKKFASEK 98 **9866665 PP >> Syntaxin_2 Syntaxin-like protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 7.4 8.6 0.0018 5.1 32 96 .. 56 120 .. 3 123 .. 0.85 Alignments for each domain: == domain 1 score: 7.4 bits; conditional E-value: 0.0018 Syntaxin_2 32 reklhelieeinelikktnellkklsqldnsssdsqqklqkekLkkdfkealkeFqkiqrqyaek 96 ekl+ iee ++ + + +++++++++ ++++ d+ +k++ ek+ ++f +l +F+++ ++++e 23726_FusoPortal_Gene_318 56 IEKLKAEIEEWKNSFLRKQADFQNFTKRKEKEVDELKKFASEKIITQFLGSLDNFERAIESSSES 120 46777777777777888888888888888888888***********************9998886 PP >> DUF4407 Domain of unknown function (DUF4407) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 7.8 4.1 0.00062 1.7 105 169 .. 12 90 .. 3 130 .. 0.64 Alignments for each domain: == domain 1 score: 7.8 bits; conditional E-value: 0.00062 DUF4407 105 ekeIdaelaeeqqeeakelkeqlaaa..............ltaelaeleaqiaalqaeidaaearldalraeaqkeldg 169 ++eI++e++e+++ e +e ke+ ++ ++e+++l+a+i+++++ + ++a++++ ++ +ke+d 23726_FusoPortal_Gene_318 12 KEEINKEVNEKKKCECEEGKEEAHEHehkndehaccgkhnHKEEIEKLKAEIEEWKNSFLRKQADFQNFTKRKEKEVDE 90 5566666666666666666666666666666666666666666666666666666666666665555555555555443 PP >> Bap31 B-cell receptor-associated protein 31-like # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 4.9 12.3 0.0059 17 126 194 .. 18 86 .. 4 93 .. 0.64 Alignments for each domain: == domain 1 score: 4.9 bits; conditional E-value: 0.0059 Bap31 126 lakleaklealkkqaesaskaakklleekeee.....kkekeeveklkkelkekekelkkakkdlealkkqaeg 194 ++ ++k+e + +++a+++ ++++e+ +++kee+eklk e++e +++ ++++d++++ k e+ 23726_FusoPortal_Gene_318 18 EVNEKKKCECEEG-----KEEAHEHEHKNDEHaccgkHNHKEEIEKLKAEIEEWKNSFLRKQADFQNFTKRKEK 86 2333333332222.....222332222222224555599999***********************999988776 PP >> DUF262 Protein of unknown function DUF262 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 5.7 10.7 0.0048 13 95 199 .. 24 149 .. 3 149 .. 0.48 Alignments for each domain: == domain 1 score: 5.7 bits; conditional E-value: 0.0048 DUF262 95 k.........................ekeeededrlkrreresdkellqkllrseklsleekeelaelnsltkkdlakeekrne 153 k ++++e+ ++lk +++e ++++l+k + +++++ +++e el+++ +++ 23726_FusoPortal_Gene_318 24 KceceegkeeahehehkndehaccgkHNHKEEIEKLKAEIEEWKNSFLRKQADFQNFTKRKEKEVDELKKFASEKI-------- 99 2333333333322222222233333444444444444444444444444444444444444333333333333333........ PP DUF262 154 eknlikayevieeilksa...........llsalktklslieieleseddedqayei 199 + + + +++ + +++l + +++i+ +l+++++ ++++ei 23726_FusoPortal_Gene_318 100 ---ITQFLGSLDN----FeraiessseskDFDSLLQGVEMIVRNLKDIMSSEDVEEI 149 ...3333222222....123345556666555666999************9999987 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (202 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 2594 (0.155218); expected 334.2 (0.02) Passed bias filter: 669 (0.0400311); expected 334.2 (0.02) Passed Vit filter: 77 (0.00460747); expected 16.7 (0.001) Passed Fwd filter: 6 (0.000359023); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 6 [number of targets reported over threshold] # CPU time: 0.28u 0.23s 00:00:00.51 Elapsed: 00:00:00.32 # Mc/sec: 1835.06 // [ok]