# hmmsearch :: search profile(s) against a sequence 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 HMM file: /home/srlab/.dammit/databases/buscodb/metazoa/hmms/BUSCOmEOG7R2BM0.hmm # target sequence database: ./run_Trinity.fasta.metazoa.busco.results//translated_proteins/Transcript_36917.fas # output directed to file: temp # per-dom hits tabular output: ./run_Trinity.fasta.metazoa.busco.results/hmmer_output/BUSCOmEOG7R2BM0.out.3 # number of worker threads: 1 # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: BUSCOmEOG7R2BM0 [M=941] Scores for complete sequences (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Sequence Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 7.8e-157 512.9 64.5 9e-157 512.8 64.5 1.0 1 transcript_3 Domain annotation for each sequence (and alignments): >> transcript_3 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 512.8 64.5 1.5e-157 9e-157 194 774 .. 2 652 .. 1 680 [. 0.87 Alignments for each domain: == domain 1 score: 512.8 bits; conditional E-value: 1.5e-157 BUSCOmEOG7R2BM0 194 edssedddgeedset...deeeseedeesevagserqsGDGqesedd........eekklDdDEDRkNPqYIPkkGtFYeHDdRtaedeeedee 276 e+++e + ge +se ++ees+++e+++ +g erqsGDGqes+dd + +klDdDEDRkNP+YIP+kG+F+eHD+R++ + transcript_3 2 EEEQEPERGEPESEGggqSKEESKAEEKENLTG-ERQSGDGQESTDDpetkaggkPGQKLDDDEDRKNPAYIPRKGLFFEHDVRGH-------A 87 334555666666666778899************.****************************************************.......8 PP BUSCOmEOG7R2BM0 277 qeeesrpkskkkklwkdkedrWdHDrfneeeQaPKSreELiaiYGYDIRneegp....PrarRrrrygrgPnkydrnwe...dEeaygkk.... 359 qeee rpk++++klwkd e+rW+HD+f+eeeQaPKSreELiaiYGYDIRn+ g +++Rr+r+ +P++ d++w+ +E++ + + transcript_3 88 QEEE-RPKGRNRKLWKD-EGRWEHDKFREEEQAPKSREELIAIYGYDIRNGGGGpgdrSYRQRRPRQSMSPSR-DKRWReggGERSVR-Swhgg 177 8777.9***********.*********************************99989999**************.99999666666633.24666 PP BUSCOmEOG7R2BM0 360 ..........rnkerkkqprgedfPalq........snknssqtvrqpprssssrreqnkkenrpplrsrednnkskersqdntr......... 426 + + + + s ++s++ +r+pp++s r+ ++++++ p++s+ +n +s+ ++ ++++ transcript_3 178 ggglnrgvppPSSIHTSS---------TsssslplsSAQRSNNPSRPPPSRS---RPSHQSQMHLPPQSHYRNDESNAPQMHPRErqgpkahse 259 655444333211133332.........133333344899*************...999***********999999999999666669******* PP BUSCOmEOG7R2BM0 427 ......................................................snqnkapnrskssaekkrsaqreeenadssgdaeekvnee 466 + q ++p r+++++ ++sa+r+ + ++s+ ++++ ++ transcript_3 260 lacdrgvasrggrgsggrggpsvvvensctsqpgveekdlsssavktvalsqpgGYQSASPRRQQKEQ--RGSADRSVSGLAASSSTSDPSLQS 351 *****************************999999999999999999988888877777777777777..899********9999999998884 PP BUSCOmEOG7R2BM0 467 sprenpqkqrpqqerqsesksysrgrrqrqkaqdnqqqqqeveekqkgrgpqprprknrnnnspqe.........eeeteelsqdlkqvnennd 551 s++ n ++++p+ er++e+ksys +rr+r+++ d +++q++vee g+ ++p + r++ + + + +el+qd++++++ + transcript_3 352 SSAAN-RESSPPTERPVERKSYSLARRTRSRPTDLGSKQPSVEETAGGNASSPGSVGGRSWTGAGDgpsqtggggVGGLTELDQDVARLSLA-G 443 44444.678***************************************************************9999999*************.* PP BUSCOmEOG7R2BM0 552 gswsqmsqsaqqne.rgvpskqhvssivpPrfqqqrqeqqqssgnrPKrYSsqRqrpvpqppappqqvqvqegy...pqqlqgqpqqhqeqpap 641 ++w q+++s+ ++e rg ps h ++ pP+f + +e++++ +nr+KrYSsqRqrpvp+p++p ++ v++g+ p+++qg+++ h+++p+p transcript_3 444 QTWNQSPTSYIRSEmRGLPSPLHIPGG-PPQF--SGMEEMGVGSNRAKRYSSQRQRPVPEPAPPM-HLGVMDGHyyePMSYQGPIYAHGDGPTP 533 ****88*********************.****..88***********************999666.6799**999******************* PP BUSCOmEOG7R2BM0 642 vpPqg...qPeialPpp.ahphqqpapsqqpalyappPvsv..pPpqqllppaylpapPavinyvpapsgpayppngyqayppaalppPpPqly 729 pPqg qPe+++P+p +hphq++ap+++palya+pPvs+ + pqqllpp+++ +pP+v+++ +y + y P Pq transcript_3 534 IPPQGmlvQPEMHIPHPgLHPHQSGAPIANPALYAGPPVSLspGQPQQLLPPPFY-PPPGVMTF-------PYT----TMY-------PTPQ-- 606 *****************************************98668999866664.*******9.......222....223.......4567.. PP BUSCOmEOG7R2BM0 730 qaqpqvtvGGitYYnpqqQqqqqpqaspqrRpkaaIPIlpPpde...v 774 +q qvt+GG+tYY+++qQ q+qp++sp+rR++++++++pPp+e transcript_3 607 -GQAQVTYGGVTYYDTMQQ-QAQPKPSPPRRTSQPVTVKPPPPEvpfA 652 .9*****************.***********************73331 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (941 nodes) Target sequences: 6 (4184 residues searched) Passed MSV filter: 3 (0.5); expected 0.1 (0.02) Passed bias filter: 1 (0.166667); expected 0.1 (0.02) Passed Vit filter: 1 (0.166667); expected 0.0 (0.001) Passed Fwd filter: 1 (0.166667); expected 0.0 (1e-05) Initial search space (Z): 6 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.04u 0.00s 00:00:00.04 Elapsed: 00:00:00.08 # Mc/sec: 49.21 // [ok]