# 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_1715.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1715 [L=376] Description: # 1840693 # 1841820 # -1 # ID=1_1715;partial=00;start_type=ATG;rbs_motif=GGA/GAG/AGG;rbs_spacer=5-10bp;gc_cont=0.235 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.1e-15 58.3 1.2 1.1e-15 58.3 1.2 1.5 2 DUF87 Domain of unknown function DUF87 0.0014 18.8 0.0 0.0036 17.5 0.0 1.7 1 AAA_22 AAA domain 0.0015 17.4 0.0 0.0027 16.6 0.0 1.3 1 TrwB_AAD_bind Type IV secretion-system coupling protein DN 0.0054 16.6 0.6 0.015 15.2 0.3 1.9 2 Rad9_Rad53_bind Fungal Rad9-like Rad53-binding ------ inclusion threshold ------ 0.02 14.1 0.0 0.033 13.4 0.0 1.4 1 FtsK_SpoIIIE FtsK/SpoIIIE family 0.048 13.0 0.3 0.095 12.0 0.3 1.5 1 SMC_N RecF/RecN/SMC N terminal domain 0.11 12.2 0.1 0.24 11.1 0.1 1.5 1 AAA_29 P-loop containing region of AAA domain 0.16 12.4 1.6 0.32 11.4 0.4 1.9 1 ABC_tran ABC transporter Domain annotation for each model (and alignments): >> DUF87 Domain of unknown function DUF87 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -3.1 0.0 3.1 6.5e+03 146 178 .. 65 97 .. 56 126 .. 0.51 2 ! 58.3 1.2 5.1e-19 1.1e-15 11 173 .. 146 305 .. 137 360 .. 0.74 Alignments for each domain: == domain 1 score: -3.1 bits; conditional E-value: 3.1 DUF87 146 lleelkkekeelrklknqaiarlldkledllrk 178 ++++++ke++++ k+++++++ + ++ ++++ 23726_FusoPortal_Gene_1715 65 IKDNYQKESKNNYKENERYEKFIYARCIGFFEN 97 344444333344444444444444444444444 PP == domain 2 score: 58.3 bits; conditional E-value: 5.1e-19 DUF87 11 evpvyldlnklvtrHlAilgsTGsGKSntvavLleelle.....kkgatililDpHGEYgkalekllsekvkvitpspsa..k 86 ++++ l nk+++ H++i+g TGsGKSnt +L + ll+ +++ + +D +GEY + l + +++ + k 23726_FusoPortal_Gene_1715 146 GLKIDLPWNKIFNTHIGIFGNTGSGKSNTLTSLYKILLNheklnLENSKFVFIDFNGEYTGDK--VLCRNKIIYELNTRQdnK 226 46667777*******************************998887788*************99..444444466665555569 PP DUF87 87 dkelkinpwsldpeelldlleeddkpaaavqrsvleealdqakakeekeklssdepveklleelkkekeelrklknqaiarll 169 ++ i+ +s+ ee+l++l+ + +++q+s+l+++++ +++ +s + ++e+ l++++++ +e ++k++ + l+ 23726_FusoPortal_Gene_1715 227 GDKFPIKSSSFWDEEVLSILFH---ATEKTQKSFLKNLVK-----GRNKYISNKNSLENYLKKIYIRCLETLNPKKEFVEILI 301 99999****************9...44556*******999.....55666666666666666666333344444444444444 PP DUF87 170 dkle 173 ++l+ 23726_FusoPortal_Gene_1715 302 RTLK 305 4444 PP >> AAA_22 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 17.5 0.0 1.7e-06 0.0036 7 65 .. 160 220 .. 158 264 .. 0.81 Alignments for each domain: == domain 1 score: 17.5 bits; conditional E-value: 1.7e-06 AAA_22 7 ilvvtGesGtGKTtllrrlleqldeeakaeeerpevvfvtlpsstsskdl..lrellralg 65 ++ ++G +G+GK+ +l l + l +++k + e++++vf++ +++ + +++ ++++ +l+ 23726_FusoPortal_Gene_1715 160 HIGIFGNTGSGKSNTLTSLYKILLNHEKLNLENSKFVFIDFNGEYTGDKVlcRNKIIYELN 220 6779************************************999985444423555555555 PP >> TrwB_AAD_bind Type IV secretion-system coupling protein DNA-binding domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 16.6 0.0 1.3e-06 0.0027 15 57 .. 158 206 .. 147 216 .. 0.83 Alignments for each domain: == domain 1 score: 16.6 bits; conditional E-value: 1.3e-06 TrwB_AAD_bind 15 tqhilivGttGsGKtqalekllrairarg......draivlDltGafvs 57 + hi i+G tGsGK+++l++l + + + + + + +D++G+++ 23726_FusoPortal_Gene_1715 158 NTHIGIFGNTGSGKSNTLTSLYKILLNHEklnlenSKFVFIDFNGEYTG 206 68******************998877776677777999999****9986 PP >> Rad9_Rad53_bind Fungal Rad9-like Rad53-binding # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -3.1 0.0 3.2 6.6e+03 63 82 .. 206 221 .. 180 233 .. 0.68 2 ! 15.2 0.3 7.2e-06 0.015 4 61 .. 306 364 .. 303 368 .. 0.89 Alignments for each domain: == domain 1 score: -3.1 bits; conditional E-value: 3.2 Rad9_Rad53_bind 63 Gdavkvdmnkvkyivrglek 82 Gd+v ++ + i+ l++ 23726_FusoPortal_Gene_1715 206 GDKVLCRNK----IIYELNT 221 777777765....4444444 PP == domain 2 score: 15.2 bits; conditional E-value: 7.2e-06 Rad9_Rad53_bind 4 lltkkdiifedavWvqyklkyYpgkilskgtsqdkslvkfeegavevkked.vkyLdlr 61 l+ +di+ ++ W +++ k+Y gki ++ + d+sl k+ e +++v++ d L++r 23726_FusoPortal_Gene_1715 306 LIDYDDIFLKKLQWHSQNQKFYIGKIYYNQGKVDESLYKHLEDRIKVENTDqFTELEIR 364 567899*****************************************998735678877 PP >> FtsK_SpoIIIE FtsK/SpoIIIE family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.4 0.0 1.6e-05 0.033 32 59 .. 148 178 .. 127 206 .. 0.75 Alignments for each domain: == domain 1 score: 13.4 bits; conditional E-value: 1.6e-05 EEEGGGS-....SEEEE--TTSSHHHHHHHH CS FtsK_SpoIIIE 32 vldlkkma...shllvaGatgsGKStalntL 59 ++dl+ + h+ + G tgsGKS++l +L 23726_FusoPortal_Gene_1715 148 KIDLPWNKifnTHIGIFGNTGSGKSNTLTSL 178 455555555567****************998 PP >> SMC_N RecF/RecN/SMC N terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.0 0.3 4.6e-05 0.095 25 83 .. 159 217 .. 144 284 .. 0.63 Alignments for each domain: == domain 1 score: 12.0 bits; conditional E-value: 4.6e-05 SEEEEEESTTSSHHHHHHHHHHHHTT.TSSS--HHHHCCCT--.SSSSSTTCCEEEECE CS xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF SMC_N 25 sftaivGpNGSGKsnlldAilFvLGersakklrseklsdlihkskekasvkkaeveitf 83 + i G+ GSGKsn l ++ +L ++++ +l ++k+ + ++ + ++ +i++ 23726_FusoPortal_Gene_1715 159 THIGIFGNTGSGKSNTLTSLYKILLNHEKLNLENSKFVFIDFNGEYTGDKVLCRNKIIY 217 344699**************666665555455555555444445555566667777766 PP >> AAA_29 P-loop containing region of AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.1 0.1 0.00011 0.24 27 53 .. 163 189 .. 155 193 .. 0.85 Alignments for each domain: == domain 1 score: 11.1 bits; conditional E-value: 0.00011 AAA_29 27 ltGpsGsGKSTllDAiqalLvpakkvr 53 + G++GsGKS +l ++ +L+ ++k++ 23726_FusoPortal_Gene_1715 163 IFGNTGSGKSNTLTSLYKILLNHEKLN 189 89**************99999988665 PP >> ABC_tran ABC transporter # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.4 0.4 0.00015 0.32 12 36 .. 159 183 .. 157 201 .. 0.82 Alignments for each domain: == domain 1 score: 11.4 bits; conditional E-value: 0.00015 SEEEEEESTTSSHHHHHHHHTTSS- CS xxxxxxxxxxxxxxxxxxxxxxxxx RF ABC_tran 12 ekvaivGenGaGKStLlkllagllk 36 + ++i+G++G+GKS l l++ l 23726_FusoPortal_Gene_1715 159 THIGIFGNTGSGKSNTLTSLYKILL 183 5799***********9999998765 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (376 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1022 (0.0611537); expected 334.2 (0.02) Passed bias filter: 421 (0.0251915); expected 334.2 (0.02) Passed Vit filter: 50 (0.00299186); expected 16.7 (0.001) Passed Fwd filter: 8 (0.000478698); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 8 [number of targets reported over threshold] # CPU time: 0.42u 0.32s 00:00:00.74 Elapsed: 00:00:00.45 # Mc/sec: 2428.99 // [ok]