# 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_438.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_438 [L=638] Description: # 483308 # 485221 # -1 # ID=1_438;partial=00;start_type=ATG;rbs_motif=GGAG/GAGG;rbs_spacer=5-10bp;gc_cont=0.208 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 6.6e-08 32.2 3.4 1.6e-07 31.0 1.4 2.2 2 RPA Replication initiator protein A Domain annotation for each model (and alignments): >> RPA Replication initiator protein A # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 31.0 1.4 9.6e-12 1.6e-07 7 198 .. 82 286 .. 79 300 .. 0.78 2 ? 0.3 0.0 0.023 3.8e+02 55 99 .. 458 500 .. 444 542 .. 0.71 Alignments for each domain: == domain 1 score: 31.0 bits; conditional E-value: 9.6e-12 RPA 7 mevplfsLakkkrtavieyeakdvkieveas.......pevGlatifDkdiliaaisqlveaknkg.rpvsrllrltayell.k 81 + +p+f+ +k++++ +y+ +d ++++++ +v ++if + i+ ++is + e++nk ++ + ++++++++++ 23726_FusoPortal_Gene_438 82 IYFPFFTPQKQNKRINFKYTFEDLGVTMYSTlipkdkkDKVFQPSIFEEKIYTFLISMYQEKTNKKiDDSEVAIEFEISDFIvN 165 5689999999999999999999966666543221121145689***********************8888899******97526 PP RPA 82 atrrdtggkqykrLkealeRLqgttiettiraggk.....ersegfslidewevvtekegrleaieitlsdwlyravl.akevl 159 + + ++ +y ++++al+ L+ t ++ +i +++k ++++f+l++ ++++ k g+ +++l + ++++++ +++++ 23726_FusoPortal_Gene_438 166 FLGNKMNRAYYAKVEQALKNLKNTIYQFEISNHTKfgknkFEDSSFQLLN-YQKL--KVGKKTFYRVILNKNIVNKIKsKRYIK 246 99************************777777666222224555666776.4544..35899999*************778898 PP RPA 160 tldpdYf..rlrkgleRrlYrlaRkhagrqeagwqisvalL 198 +++ + ++ ++ r+Y+ + k ++++ +i+v++L 23726_FusoPortal_Gene_438 247 YNTKNLLeiMIKDPIASRIYKYISKIRYKNNK-GEINVRTL 286 888888855789********988777666666.56666666 PP == domain 2 score: 0.3 bits; conditional E-value: 0.023 RPA 55 isqlveaknkgrpvsrllrltayellkatrrdtggkqykrLkeal 99 is ++ kn g + s+ ++ t + + ++ g k+++++k+a 23726_FusoPortal_Gene_438 458 ISGIL--KNIGGKGSKKVKNTNMTIFENVSKGKGLKNKSQIKQAR 500 34443..45675668999999999999999999999999999874 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (638 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1852 (0.110819); expected 334.2 (0.02) Passed bias filter: 359 (0.0214816); expected 334.2 (0.02) Passed Vit filter: 22 (0.00131642); expected 16.7 (0.001) Passed Fwd filter: 1 (5.98372e-05); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.31u 0.22s 00:00:00.53 Elapsed: 00:00:00.32 # Mc/sec: 5795.90 // [ok]