Mercurial > repos > iuc > hmmer3
diff test-data/MADE1.out @ 0:62479bdcc059 draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/hmmer3 commit 4164b44c651bcbdac6637eccce61b2a802c9b569
author | iuc |
---|---|
date | Tue, 12 May 2015 15:04:26 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/MADE1.out Tue May 12 15:04:26 2015 -0400 @@ -0,0 +1,86 @@ +# hmmscan :: search sequence(s) against a profile 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 sequence file: test-data/dna_target.fa +# target HMM database: test-data/MADE1.hmm +# per-seq hits tabular output: t.tbl +# per-dom hits tabular output: t.domtbl +# pfam-style tabular hit output: t.pfamtbl +# max ASCII text line length: unlimited +# Vit filter P threshold: <= 0.001 +# Fwd filter P threshold: <= 1e-05 +# random number seed set to: 4 +# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + +Query: humanchr1/239220001-239550000 [L=330000] +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 + ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- + 7.1e-18 51.5 28.5 2e-12 34.0 0.7 9.6 5 MADE1 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon + + +Domain annotation for each model (and alignments): +>> MADE1 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon + # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc + --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- + 1 ? -4.3 0.1 1 1 30 54 .. 80044 80068 .. 80030 80073 .. 0.80 + 2 ? -6.7 3.3 1 1 13 71 .. 154012 154072 .. 154011 154076 .. 0.75 + 3 ! 27.6 0.7 2e-10 2e-10 1 44 [. 174456 174514 .. 174456 174577 .. 0.62 + 4 ! 34.0 0.7 2e-12 2e-12 2 80 .] 302388 302466 .. 302387 302466 .. 0.86 + 5 ? 2.8 0.7 0.011 0.011 27 75 .. 304060 304107 .. 304021 304109 .. 0.61 + + Alignments for each domain: + == domain 1 score: -4.3 bits; conditional E-value: 1 + xxxxxxxxxxxxxxxxxxxxxxxxx RF + MADE1 30 ttgccattacttttaatggcaaaaa 54 + t g catt ttt aatggcaaa a + humanchr1/239220001-239550000 80044 TAGTCATTCATTTCAATGGCAAATA 80068 + 45789****************9966 PP + + == domain 2 score: -6.7 bits; conditional E-value: 1 + xxxxxxxxxxxxxxxxxxxxxxxxx......xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF + MADE1 13 aaaagtaattgcggtttttgccatt......acttttaatggcaaaaaccgcaattacttttgca 71 + aaaagta tt + ttttgc att a tttaa gcaaa a + tta tttgca + humanchr1/239220001-239550000 154012 AAAAGTAGTTTTCAATTTTGCAATTtgaccaATATTTAAATGCAAATATT----TTATATTTGCA 154072 + 78999999999999999999999984444444457777777899998876....77777888876 PP + + == domain 3 score: 27.6 bits; conditional E-value: 2e-10 + xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx...............x RF + MADE1 1 ttaggttggtgcaaaagtaattgcggtttttgccattactttt...............a 44 + ttaggtt gtgcaaaagtaattg+ggtttttg cattactttt a + humanchr1/239220001-239550000 174456 TTAGGTTAGTGCAAAAGTAATTGTGGTTTTTGTCATTACTTTTctgcatgctagaagtA 174514 + 79***************************************964443333333333330 PP + + == domain 4 score: 34.0 bits; conditional E-value: 2e-12 + xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF + MADE1 2 taggttggtgcaaaagtaattgcggtttttgccattacttttaatggcaaaaaccgcaattacttttgcaccaacctaa 80 + t ggt ggtgcaaaa aattg+ggtttttgccatt cttttaat gc a + a t ctttt caccaa ctaa + humanchr1/239220001-239550000 302388 TTGGTCGGTGCAAAATCAATTGTGGTTTTTGCCATTGCTTTTAATTGCTTTTAAAAGTAATGCTTTTACACCAATCTAA 302466 + 56899******************************************963333233345578**************997 PP + + == domain 5 score: 2.8 bits; conditional E-value: 0.011 + xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF + MADE1 27 tttttgccattacttttaatggcaaaaaccgcaattacttttgcaccaa 75 + tttt g c ta tt a tgg aaaaa ++ca tta ttttgca aa + humanchr1/239220001-239550000 304060 TTTTAGACTATA-GTTAAGTGGGAAAAAATACACTTATTTTTGCATTAA 304107 + 222222222222.3455779************************98765 PP + + + +Internal pipeline statistics summary: +------------------------------------- +Query sequence(s): 1 (330000 residues searched) +Target model(s): 1 (80 nodes) +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: 0.52u 0.03s 00:00:00.55 Elapsed: 00:00:00.54 +# Mc/sec: 48.89 +// +[ok]