# 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_506.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_506 [L=1888] Description: # 554807 # 560470 # 1 # ID=1_506;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.265 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.7e-30 104.6 8.1 2.7e-30 104.6 8.1 7.5 4 Haemagg_act haemagglutination activity domain 1.2e-09 38.5 35.2 2.3e-07 31.0 8.8 14.3 13 Fil_haemagg_2 Haemagluttinin repeat ------ inclusion threshold ------ 1.3 9.2 3.3 7.1 6.9 0.2 3.3 3 Cytochrome_P460 Cytochrome P460 Domain annotation for each model (and alignments): >> Haemagg_act haemagglutination activity domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 104.6 8.1 4.9e-34 2.7e-30 3 115 .. 25 156 .. 23 160 .. 0.97 2 ? -60.0 68.0 3 1.7e+04 21 108 .. 614 707 .. 298 709 .. 0.83 3 ? -3.9 0.9 2.1 1.2e+04 7 35 .. 1169 1195 .. 1164 1197 .. 0.65 4 ? 1.4 0.3 0.047 2.6e+02 45 81 .. 1693 1730 .. 1663 1736 .. 0.84 Alignments for each domain: == domain 1 score: 104.6 bits; conditional E-value: 4.9e-34 TSEECCCEEEC-EEEETCS-ECCEEEEEETTSEEEEEECE--B-TTEEEEEE-S...................SSTTSEEEEEE CS Haemagg_act 3 iapdgtlvtnvkakkagngvpivniatpqkggnsfnsfseFnvpeggvvlfnaq...................ngaasvilnrV 67 i+pd+++ n +++ka+ngvp+vni +p ++g+s+n ++++nv+ g++l+n++ +++as+i+n+V 23726_FusoPortal_Gene_506 25 IVPDNAASRNLQVDKAANGVPLVNIEAPDNNGTSHNVYKDYNVDGRGAILNNSKdltnsqlggliygnpnlqnSKEASTIINEV 108 79**************************************************736***************************** PP -SSS-EEEECEEEEEECSCEEEEE-TT-EEEECCEEECECEEEEESSE CS Haemagg_act 68 tggspSkIqGaievlgakanlflaNpnGiifggnatiniggsfvttta 115 g ++S+I+G++e++g+kan++laNpnGi+++g ++in+g++++tt+ 23726_FusoPortal_Gene_506 109 SGVNRSRIEGYQEIAGKKANYILANPNGIYINGAGFINTGNVTLTTGS 156 **********************************************96 PP == domain 2 score: -60.0 bits; conditional E-value: 3 S-ECCEEEEEETTS.EEEEEECE--B-TTEEEEEE-S..........SSTTSEEEEEE.-SSS-EEEECEEEEEECSCEEE.EE CS Haemagg_act 21 gvpivniatpqkgg.nsfnsfseFnvpeggvvlfnaq..........ngaasvilnrV.tggspSkIqGaievlgakanlf.la 91 +++i++ + ++g+ ++ + + ++ +++++ + +++ ++i + + +s +I+G+i+ l+a+ n + 23726_FusoPortal_Gene_506 614 NTEITKLT--NDGNiVVKENLKAKDITNSNTI----KvgenlntdklQNSKTLIAKNInIEKSLNNINGKITSLNANINTSdIK 691 33444444..2222022233444444444444....034456677778888888888855677889******999888876599 PP -TT-EEEECCEEECECE CS Haemagg_act 92 NpnGiifggnatinigg 108 N nGii ++ ++ini+ 23726_FusoPortal_Gene_506 692 NNNGII-QAIKNINIKT 707 999996.7788888875 PP == domain 3 score: -3.9 bits; conditional E-value: 2.1 CCCEEEC-EEEETCS-ECCEEEEEETTSE CS Haemagg_act 7 gtlvtnvkakkagngvpivniatpqkggn 35 +++ +++ k+ +ngv++ i +++k++n 23726_FusoPortal_Gene_506 1169 NQEINST--KVDKNGVNLDSINIDKKNTN 1195 3333333..46788999999988888876 PP == domain 4 score: 1.4 bits; conditional E-value: 0.047 B-TTEEEEEE-.SSSTTSEEEEEE-SSS-EEEECEEEE CS Haemagg_act 45 vpeggvvlfna.qngaasvilnrVtggspSkIqGaiev 81 +++ v ++ +g+++vi + + g+ S I+G++++ 23726_FusoPortal_Gene_506 1693 AESNLYVKNQGiIGGDIKVIGSNLVLGNNSIINGKLTT 1730 5556667777766999*******************996 PP >> Fil_haemagg_2 Haemagluttinin repeat # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -8.5 24.6 3 1.7e+04 22 164 .. 258 393 .. 251 399 .. 0.58 2 ? -14.0 20.5 3 1.7e+04 13 101 .. 383 463 .. 377 494 .. 0.62 3 ? -17.5 28.4 3 1.7e+04 8 128 .. 443 562 .. 437 585 .. 0.67 4 ? -23.3 31.8 3 1.7e+04 13 167 .. 502 664 .. 495 667 .. 0.65 5 ? -16.8 29.4 3 1.7e+04 55 147 .. 638 726 .. 565 733 .. 0.56 6 ? -8.9 12.7 3 1.7e+04 34 130 .. 734 830 .. 716 871 .. 0.53 7 ? -3.5 0.1 1.6 9e+03 21 35 .. 899 914 .. 894 944 .. 0.71 8 ! 19.5 2.6 1.4e-07 0.00078 12 59 .. 1026 1072 .. 1018 1083 .. 0.91 9 ? 3.7 5.7 0.01 57 17 56 .. 1135 1180 .. 1113 1199 .. 0.48 10 ! 6.1 7.8 0.0018 10 16 66 .. 1417 1464 .. 1406 1468 .. 0.82 11 ! 31.0 8.8 4.2e-11 2.3e-07 5 76 .. 1529 1598 .. 1526 1601 .. 0.89 12 ! 19.7 8.0 1.2e-07 0.00069 3 82 .. 1601 1678 .. 1599 1686 .. 0.76 13 ! 13.6 13.3 9.4e-06 0.052 13 167 .. 1682 1826 .. 1676 1833 .. 0.80 Alignments for each domain: == domain 1 score: -8.5 bits; conditional E-value: 3 Fil_haemagg_2 22 inaggnvtikatg..gditlqGsqvsageevalnaandinleaaenteettstnksssasagvsvglggggkaagvgasasggk 103 i+++g+v +k t+ g+i+++ ++ + + l a+n i ++ + ++ + + n + ++s v + + +++ +++s++ + 23726_FusoPortal_Gene_506 258 ISSKGKVYLKDTQakGNIKITSTETEIAN--KLLAENLIDIKGKTTNSGQIQANNNITISGNVDSSNLIS-TNKDITISGNLTN 338 45677777777756678888888888777..555788888887777766666666666666665543333.3566777777777 PP Fil_haemagg_2 104 gkengdsttqtnshvaaggtvtlnsgkDtnlkGanvsgkkvnadvggdLnieSlqDtstyn 164 + e +++ t+++++++g++t+ ++n +++ +++ ++ vg +n + l +tst + 23726_FusoPortal_Gene_506 339 SGEVSTKN-LTTNNLDNKGNITVIN--NVN---SELITNNGKLLVGDTINSQNLTNTSTVQ 393 77763333.3334454555555443..333...4445555555555555555555555555 PP == domain 2 score: -14.0 bits; conditional E-value: 3 Fil_haemagg_2 13 ssttavgstinaggnvtikat.ggditlqGsqvsageevalnaandinleaaenteettstnksssasagvsvglggggkaagv 95 s++ +++st++ g++++ik++ +++ +++ ++s ++ +++ ++i ++ + t+e tn+ + +s +s ++++k ++ 23726_FusoPortal_Gene_506 383 SQNLTNTSTVQ-GKTLDIKSKvNSSGKIISDNISTKD---ISNSGNISSKSIT-TQE--LTNSGEIISNNLSSNNTNNSK--NI 457 44455677776.8999999996578888888888777...6678889888888.433..567777666666665555544..44 PP Fil_haemagg_2 96 gasasg 101 ++++ 23726_FusoPortal_Gene_506 458 FVNGNL 463 444433 PP == domain 3 score: -17.5 bits; conditional E-value: 3 Fil_haemagg_2 8 sestsssttavgstinaggnvtikat.ggditlqGsqvsageevalnaandinleaaenteettstnksssasagvsvglgggg 90 s++ ss++t ++ +i ++gn++i + +++ ++G +++ ++ ++++++i++++ ++ ++++ ++++ag+ ++++++ 23726_FusoPortal_Gene_506 443 SNNLSSNNTNNSKNIFVNGNLKIVNNlNNSGIIEGLELNTNS---IDNTGNITIKNKLTSQNLNNKKNTANINAGF-LDVQNKI 522 455588999999*************96456667888888555...8889999999999777777777778888886.4544444 PP Fil_haemagg_2 91 kaa.gv.gasasggkgkengdsttqtnshvaaggtvtlns 128 +++ ++ ++++++++ +++g+ t+ ++ ++ ++ ++ + 23726_FusoPortal_Gene_506 523 SSVgNIkAITLKTNNLDNSGNILTNSLTTTENINKGSIIA 562 3331121234444444444446666667776666665555 PP == domain 4 score: -23.3 bits; conditional E-value: 3 Fil_haemagg_2 13 ssttavgstinaggnvtikat..g....gditlqGsqvsage...evalnaandinlea..aenteettstnksssasagvsvg 85 ++++++ +ina g ++++ + ++itl+ ++++ + +l ++++in + a+n ++++ n+ s +s ++v 23726_FusoPortal_Gene_506 502 LNNKKNTANINA-GFLDVQNKisSvgniKAITLKTNNLDNSGnilTNSLTTTENINKGSiiAKNISSQNLVNSGSVISNNITVA 584 444555556665.44454444431335677777777776664122245667777776663355555555567999999999998 PP Fil_haemagg_2 86 lggggkaagvgasasggkgk.engdsttqtnshvaag.gtvtlnsgkDtnlkGanvsgkkvnadvggdLnieSlqDtstynsks 167 + ++ + ++ a+ + + +k +n+++ ++n+++++ ++ ++ ++++++ k ++ ++vg++Ln + lq+ +t +k+ 23726_FusoPortal_Gene_506 585 KNITN-TNSIFANEKISADKiSNSNKLVAKNTEITKLtNDGNIVVKENLKAK---DITNSNTIKVGENLNTDKLQNSKTLIAKN 664 88876.5556665555555436678899999999855488999998888744...44788999************998866555 PP == domain 5 score: -16.8 bits; conditional E-value: 3 Fil_haemagg_2 55 andinleaaenteettstnksssasagvsvglggggkaagvgasasggkgkengdsttqtnshvaaggtvtlnsgkDtnlkGan 138 +n+ +++ +en ++++ +n++ + +++ + + +++ ++ + + n+++ +++n+ ++a +++++++ +D++l+G 23726_FusoPortal_Gene_506 638 TNSNTIKVGENLNTDKLQNSKTLIAKNINIEKSLN----NINGKITSLNANINTSDIKNNNGIIQAIKNINIKTSNDLSLDGKY 717 23333444444444444445555555555555554....355555555555555567899999999999999999999999965 PP Fil_haemagg_2 139 vsgkkvnad 147 + n++ 23726_FusoPortal_Gene_506 718 TANDSLNIN 726 544555555 PP == domain 6 score: -8.9 bits; conditional E-value: 3 Fil_haemagg_2 34 g.......gditlqGsqvsageevalnaandinleaaenteettstnksssasagvsvglggggkaagvgasasggkgkengds 110 g +++l+G+ ++ ++ +++ +++ ++a e ++ + + + + +s+++++ ++ ++ ++ + +k ++++ 23726_FusoPortal_Gene_506 734 GnlendgkINLNLTGNLINNKR---ISSSGNLDITAKEV--LNNGKDSAIGSETNLSITANSLKNEGNLLFGVRTD-NKLKTTG 811 0333332346667777777555...44455555555442..22233333333444555555544322222222222.2222344 PP Fil_haemagg_2 111 ttqtnshvaaggtvtlnsgk 130 +++ +++ g++++++ k 23726_FusoPortal_Gene_506 812 NITNKGVIGSLGKLSIEA-K 830 455566666666666666.3 PP == domain 7 score: -3.5 bits; conditional E-value: 1.6 Fil_haemagg_2 21 tinaggnvtikatg.g 35 i++++n++i+a++ + 23726_FusoPortal_Gene_506 899 DIESEKNIKIEAKDiK 914 6888889999888542 PP == domain 8 score: 19.5 bits; conditional E-value: 1.4e-07 Fil_haemagg_2 12 sssttavgstinaggnvtikatggditlqGsqvsageevalnaandin 59 +s+ t++ sti+aggn+t+kat+ di s++ a+ +v++ a+n in 23726_FusoPortal_Gene_506 1026 ESNYTSEMSTIKAGGNITLKATN-DIENLESNILANGNVDITAKNLIN 1072 678899****************9.*******************99887 PP == domain 9 score: 3.7 bits; conditional E-value: 0.01 Fil_haemagg_2 17 avgstinaggnvtikatg........gditlqGsqvsageevalnaan 56 ++ +i+aggn++i+a++ +++++ + + ++ + + ++n 23726_FusoPortal_Gene_506 1135 GKNAKISAGGNIKIEANKvgngleikENVSIKSKNQEINS-TKV-DKN 1180 2334455555555555532223333333333333333333.333.333 PP == domain 10 score: 6.1 bits; conditional E-value: 0.0018 Fil_haemagg_2 16 tavgstinaggnvtikatggditlqGsqvsageevalnaandinleaaent 66 a++ +i g++v+i+a++ +++ +G++++a+e++++ +a+di +++ t 23726_FusoPortal_Gene_506 1417 GANKASID-GNTVNITAKN-SVDNIGADIKAKENLTV-TAKDISNLSTKRT 1464 56777887.8********9.*************9999.8888888888766 PP == domain 11 score: 31.0 bits; conditional E-value: 4.2e-11 Fil_haemagg_2 5 ksksestsssttavgstinaggnvtikatggditlqGsqvsageevalnaandinleaaenteettstnkss 76 s + +t ++++ +s+i+a +n++i+a++ di+++Gs+v a++e +++a++d+n+ +a+++ +++s+ 23726_FusoPortal_Gene_506 1529 NSNNYTTIDNKSNIASNIEA-KNINIEAKK-DIDIKGSNVVAKNEANIKADGDVNIVSATDSRFYAHKETSK 1598 567889999**********9.********8.*****************************988766666655 PP == domain 12 score: 19.7 bits; conditional E-value: 1.2e-07 Fil_haemagg_2 3 asksksestssstta.vgstinaggnvtikatggditlqGsqvsageevalnaandinleaaenteettstnksssasagv 82 ++ks+se + t v+s+i g++v+i++++ ++++ Gs+v a+++ +++a++ i+ aa++ +++ ++k +++ +g+ 23726_FusoPortal_Gene_506 1601 FGKSSSEENIAYATRnVASNI-IGDKVNITSGK-NVNIFGSNVGAKDTGNISAKGTITEAAAKEINYS-YHKKTKTGFMGL 1678 678888888666655167777.69********8.*************99999999**99999977764.444444444444 PP == domain 13 score: 13.6 bits; conditional E-value: 9.4e-06 Fil_haemagg_2 13 ssttavgstinaggnvtikatg...gditlqGsqvsageevalnaandinleaaenteettstnksssasagvsvglggggk 91 ss ++ +++na+ n+ +k +g gdi+++Gs++ g++ ++ +++++ + + +++ + ++ + g+s ++ggg 23726_FusoPortal_Gene_506 1682 SSAEKIRQELNAESNLYVKNQGiigGDIKVIGSNLVLGNNSII--NGKLTTDSNQ----LHNSYSYEESKKGFSGSIGGG-- 1755 5666778899***********7788***********9995555..4444444444....444454444567888888888.. PP Fil_haemagg_2 92 aagvgasasggkgkengdsttqtnshvaaggtvtlnsgkDtnlkGanvsgkkvnadvggdLnieSlqDtstynsks 167 g++++++++++ +++s+ + +s++ g+ + ln g D + +n++ +k+ ++ gd++ ++Dt++ + s 23726_FusoPortal_Gene_506 1756 --GFSVGYGKSESGLKEKSVINAKSNLVLGDGTVLNKGAD--ITATNLTHGKITVNN-GDVKFGARKDTKDVETYS 1826 ..599*******************************9888..5689**********9.***********7665332 PP >> Cytochrome_P460 Cytochrome P460 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -2.5 0.1 1.1 5.9e+03 36 59 .. 550 573 .. 545 616 .. 0.66 2 ? 0.6 0.2 0.11 6.4e+02 40 62 .. 612 641 .. 575 659 .. 0.65 3 ? 6.9 0.2 0.0013 7.1 71 98 .. 989 1016 .. 966 1029 .. 0.80 Alignments for each domain: == domain 1 score: -2.5 bits; conditional E-value: 1.1 Cytochrome_P460 36 rkggefpdGtvlvkelykakkske 59 ++++ +G++++k + +++ ++ 23726_FusoPortal_Gene_506 550 TTTENINKGSIIAKNISSQNLVNS 573 556667888888888777665543 PP == domain 2 score: 0.6 bits; conditional E-value: 0.11 Cytochrome_P460 40 e........fpdGtvlvkelykakkskeasg 62 dG ++vke +kak+ ++ s+ 23726_FusoPortal_Gene_506 612 AknteitklTNDGNIVVKENLKAKDITN-SN 641 023334457899***********99884.33 PP == domain 3 score: 6.9 bits; conditional E-value: 0.0013 Cytochrome_P460 71 gvevmvKdskrfaatggWgfarfsaekk 98 ++v Kd+++f+ g W + +++++ 23726_FusoPortal_Gene_506 989 SFFVLKKDDEKFHTYGAWDWEGYKSKAG 1016 5788999****************97754 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (1888 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1058 (0.0633078); expected 334.2 (0.02) Passed bias filter: 180 (0.0107707); expected 334.2 (0.02) Passed Vit filter: 21 (0.00125658); expected 16.7 (0.001) Passed Fwd filter: 7 (0.000418861); 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.72u 0.35s 00:00:01.06 Elapsed: 00:00:00.45 # Mc/sec: 12196.61 // [ok]