Bio-FastParsers
view release on metacpan or search on metacpan
test/hmmer3.stdout view on Meta::CPAN
# hmmsearch :: search profile(s) against a sequence database
# HMMER 3.1b1 (May 2013); http://hmmer.org/
# Copyright (C) 2013 Howard Hughes Medical Institute.
# Freely distributed under the GNU General Public License (GPLv3).
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# query HMM file: test.hmm
# target sequence database: Meredith169AA-Abrocoma_bennettii.stock
# output directed to file: result_outclassic.res
# max ASCII text line length: unlimited
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: Meredith169AA [M=6705]
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
------- ------ ----- ------- ------ ----- ---- -- -------- -----------
0 12909.9 47.3 0 4943.5 15.5 3.0 3 Abrocoma_bennettii
Domain annotation for each sequence (and alignments):
>> Abrocoma_bennettii
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 3876.6 2.8 0 0 1 1789 [. 1 1788 [. 1 1795 [. 1.00
2 ! 4102.0 13.3 0 0 2038 4546 .. 2037 4459 .. 2026 4465 .. 0.99
3 ! 4943.5 15.5 0 0 4778 6705 .] 4691 6619 .] 4683 6619 .] 1.00
Alignments for each domain:
== domain 1 score: 3876.6 bits; conditional E-value: 0
Meredith169AA 1 cqfrsvmrmdymvyfsfftwiliplvvmcaiyldifyvirnklsqnfsssketgafygrefktakslflvlflfalswlplsiinciiyfpqvvlylgillshflilftifgnalvilavltsrslrapqnlflvslaaadilvatliipfslanellgywyfrrtwcevylaldvlfctssivhlcaisldrywavsraleynskrtprrikciiltvwliaa...
c+f +m+m+ymv+fsfftwi+iplvvmcaiyl+if+virnklsqn+s+sketgafygref+takslflvlflfal+wlplsiinc++yfp+ +++lgillsh i gnal ilavltsrslrapqnlflvslaaadilvatliipfslanellgywyfrrtwcevylaldvlfctssivhlc isldrywavsraleynskrtprrikciiltvwliaa...
Abrocoma_bennettii 1 CSFLPTMSMEYMVFFSFFTWIFIPLVVMCAIYLNIFHVIRNKLSQNLSASKETGAFYGREFRTAKSLFLVLFLFALCWLPLSIINCALYFPDDIMFLGILLSHXXXXXXILGNALGILAVLTSRSLRAPQNLFLVSLAAADILVATLIIPFSLANELLGYWYFRRTWCEVYLALDVLFCTSSIVHLCXISLDRYWAVSRALEYNSKRTPRRIKCIILTVWLIAA...
99******************************************************************************************************************************************************************************************************************************...
== domain 2 score: 4102.0 bits; conditional E-value: 0
Meredith169AA 2038 qkvnewfsrsdeiltsddshdggsekidllasdphdaliceservhskpvesniedkifgktyrrkaslpnlshiteitqerpltnklkrkrrttsglhpedfikkvdlavvqktpekinegtdqngqvmnitnnghenetkgdyvqkeknanpteslekesafrtkaepisssisnmelelnihsskapkknrlrrksstrhihalelvvsrnpsppnhtelq...
vnewfsrsde+ ts+d +d + ek+dl asdp++a++c++ r ++ + + ++edk+fgktyrrk sl nlsh + +tqe+plt+++krkrrt s+lhpedfikk+ la +qktpek ++g++q++qv++i nn+ enet+g vq+ +n++p++slekes+ + +ep+sss ++mele+n++ ap+ nrlr++ss+r+++ + spp h elq...
Abrocoma_bennettii 2037 XXVNEWFSRSDEMETSEDGCDWEXEKTDLKASDPQRAVTCATVRGCPSAGKGSVEDKVFGKTYRRKGSLANLSHSAXVTQECPLTSTVKRKRRTLSCLHPEDFIKKAGLAAAQKTPEKDTQGMKQSSQVTSIPNNN-ENETEGVAVQRXQNPSPVQSLEKESVPGMGTEPQSSSGNSMELEWNVCPPGAPQGNRLRKRSSARQVP--------HLSPPGHPELQ...
679**********************************************************************************************************************************985.9*******************************************************************........779********...
== domain 3 score: 4943.5 bits; conditional E-value: 0
Meredith169AA 4778 pinssgtksafvtvrvldtpsppvnlkvteitkdsvsitwepplldggskiknyivekreatrksYaAvvTnchKnsWKidqLqeGcsYyFRVtAeneYGiGlpartadpiKvaeVpqPpgKitvdDVTRnSVSLsWtKPEHDGGSKIiqYiVEMQAKhseKWSeCARVKsLEAVITNLtQGEEYLFRVvAVNEKGRSDPRSLAVPivAKDLVIEPDVRPaFns...
nssg+ksafvtvrvldtpsppvnlkvteitkdsvsitwepplldggskiknyivekreatrksYaAvvTnchKnsWKidqLqeGcsYyFRVtAeneYGiGlpartadpiKvaeVpqPpgKitvdDVTRnSVSLsWtKPEHDGGSKIiqYiVEMQAKhseKWSeCARVKsLEAVITNLtQGEEYLFRVvAVNEKGRSDPRSLAVPi+AKDLVIEPDVRPaF+s...
Abrocoma_bennettii 4691 XXNSSGMKSAFVTVRVLDTPSPPVNLKVTEITKDSVSITWEPPLLDGGSKIKNYIVEKREATRKSYAAVVTNCHKNSWKIDQLQEGCSYYFRVTAENEYGIGLPARTADPIKVAEVPQPPGKITVDDVTRNSVSLSWTKPEHDGGSKIIQYIVEMQAKHSEKWSECARVKSLEAVITNLTQGEEYLFRVVAVNEKGRSDPRSLAVPIIAKDLVIEPDVRPAFSS...
57******************************************************************************************************************************************************************************************************************************...
Internal pipeline statistics summary:
-------------------------------------
Query model(s): 1 (6705 nodes)
Target sequences: 1 (6619 residues searched)
Passed MSV filter: 1 (1); expected 0.0 (0.02)
Passed bias filter: 1 (1); expected 0.0 (0.02)
Passed Vit filter: 1 (1); expected 0.0 (0.001)
Passed Fwd filter: 1 (1); expected 0.0 (1e-05)
Initial search space (Z): 1 [actual number of targets]
Domain search space (domZ): 1 [number of targets reported over threshold]
# CPU time: 1.22u 0.11s 00:00:01.33 Elapsed: 00:00:01.32
# Mc/sec: 33.62
//
[ok]
( run in 1.047 second using v1.01-cache-2.11-cpan-df04353d9ac )