# HG changeset patch # User iuc # Date 1523087515 14400 # Node ID 1a83249ddfffbb027e3df86902ee7e85b4ccec3c # Parent 25fa92aa34c0d2ab000643c958470635310c9e61 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/hmmer3 commit fa7dec5f222510d58f566f4799a04e3731fa03f6 diff -r 25fa92aa34c0 -r 1a83249ddfff hmmscan.xml --- a/hmmscan.xml Wed Apr 04 14:02:12 2018 -0400 +++ b/hmmscan.xml Sat Apr 07 03:51:55 2018 -0400 @@ -1,17 +1,13 @@ - - search sequence(s) against a profile database + + search protein sequence(s) against a protein profile database macros.xml '$output' -#if $input_hmm_conditional.input_hmm_source == "history": - ## Remove pressed database. Eventually will be replaced with a dedicated tool/datatype - && rm '${input_hmm_conditional.hmmfile}.h3m' '${input_hmm_conditional.hmmfile}.h3i' '${input_hmm_conditional.hmmfile}.h3f' '${input_hmm_conditional.hmmfile}.h3p'; -#end if - ]]> + ]]> @@ -42,15 +34,15 @@ - - - 'tblout' in oformat + + + oformat and 'tblout' in oformat - - 'domtblout' in oformat + + oformat and 'domtblout' in oformat - - 'pfamtblout' in oformat + + oformat and 'pfamtblout' in oformat @@ -61,9 +53,9 @@ - - - + + + + ]]> diff -r 25fa92aa34c0 -r 1a83249ddfff macros.xml --- a/macros.xml Wed Apr 04 14:02:12 2018 -0400 +++ b/macros.xml Sat Apr 07 03:51:55 2018 -0400 @@ -2,11 +2,11 @@ - hmmer + hmmer - 0.1 + 3.1b2 @@ -21,64 +21,64 @@ -E $E --domE $domE -#if $T: --T $T +#if str($T): + -T $T #end if -#if $domT: ---domT $domT +#if str($domT): + --domT $domT #end if -#if $incE: ---incE $incE +#if str($incE): + --incE $incE #end if -#if $incT: ---incT $incT +#if str($incdomE): + --incdomE $incdomE #end if -#if $incdomE: ---incdomE $incdomE +#if str($incT): + --incT $incT #end if -#if $incdomT: ---incdomT $incdomT +#if str($incdomT): + --incdomT $incdomT #end if - - - - + + + + - - - - + + + + -E $E -#if $T: --T $T +#if str($T): + -T $T #end if -#if $incE: ---incE $incE +#if str($incE): + --incE $incE #end if -#if $incT: ---incT $incT +#if str($incT): + --incT $incT #end if - - + + - - + + $max @@ -86,15 +86,14 @@ --F2 $F2 --F3 $F3 $nobias - - - - - - + + + + + --EmL $EmL @@ -107,90 +106,90 @@ - - - - - - - + + + + + + + #if 'tblout' in str($oformat): - --tblout $tblout + --tblout '$tblout' #end if #if 'domtblout' in str($oformat): - --domtblout $domtblout + --domtblout '$domtblout' #end if $acc $noali $notextw - + - - - + + + #if 'tblout' in str($oformat): - --tblout $tblout + --tblout '$tblout' #end if #if 'domtblout' in str($oformat): - --domtblout $domtblout + --domtblout '$domtblout' #end if #if 'pfamtblout' in str($oformat): - --pfamtblout $pfamtblout + --pfamtblout '$pfamtblout' #end if $acc $noali $notextw - + - - - + + + - + #if 'tblout' in str($oformat): - --tblout $tblout + --tblout '$tblout' #end if #if 'dfamtblout' in str($oformat): - --dfamtblout $dfamtblout + --dfamtblout '$dfamtblout' #end if #if 'aliscoresout' in str($oformat): - --aliscoresout $aliscoresout + --aliscoresout '$aliscoresout' #end if $acc $noali $notextw - + - - - + + + #if $hssi.hssi_select == "singlemx": @@ -205,11 +204,10 @@ + - - - - + + @@ -222,7 +220,7 @@ --seed $seed - + @@ -230,19 +228,19 @@ $nonull2 -#if $Z: --Z $Z +#if str($Z): + -Z $Z #end if -#if $domZ: ---domZ $domZ +#if str($domZ): + --domZ $domZ #end if - - - + + + $input_format_select @@ -264,7 +262,7 @@ $arsws.arsws_select #if $arsws.arsws_select == "--wblosum": ---wid $arsws.wid + --wid $arsws.wid #end if @@ -282,7 +280,7 @@ - + @@ -315,13 +313,13 @@ - - - + + + - - + + @@ -333,17 +331,19 @@ $cut_tc - - - + + + --$mcs.model_construction_strategy_select #if $mcs.model_construction_strategy_select == "fast": ---symfrac $mcs.symfrac + --symfrac $mcs.symfrac #end if - +#if str($fragthresh) + --fragthresh $fragthresh +#end if @@ -353,12 +353,11 @@ - + - - + $aps_select @@ -376,32 +375,32 @@ -#if $w_beta: ---w_beta $w_beta +#if str($w_beta): + --w_beta $w_beta #end if -#if $w_length: ---w_length $w_length +#if str($w_length): + --w_length $w_length #end if - - - + + - + + ]]> - + @@ -413,12 +412,12 @@ - + - + > MADE1 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon + score bias Evalue hmmfrom hmm to alifrom ali to envfrom env to mod len acc + ------ ----- --------- ------- ------- --------- --------- --------- --------- --------- ---- + ! 38.6 7.4 1.2e-10 4 80 .] 302390 302466 .. 302387 302466 .. 80 0.87 + + Alignment: + score: 38.6 bits + xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx....xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF + MADE1 4 ggttggtgcaaaagtaattgcggtttttgccattacttttaatggc....aaaaaccgcaattacttttgcaccaacctaa 80 + ggt ggtgcaaaa aattg ggtttttgccatt cttttaat gc a aaa g a t ctttt caccaa ctaa + humanchr1/239220001-239550000 302390 GGTCGGTGCAAAATCAATTGTGGTTTTTGCCATTGCTTTTAATTGCttttA-AAA--GTA-ATGCTTTTACACCAATCTAA 302466 + 899******************************************955533.443..334.4689***********99986 PP + +>> MADE1 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon + score bias Evalue hmmfrom hmm to alifrom ali to envfrom env to mod len acc + ------ ----- --------- ------- ------- --------- --------- --------- --------- --------- ---- + ! 29.6 8.3 7.8e-08 1 43 [. 174456 174498 .. 174456 174518 .. 80 0.92 + + Alignment: + score: 29.6 bits + xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF + MADE1 1 ttaggttggtgcaaaagtaattgcggtttttgccattactttt 43 + ttaggtt gtgcaaaagtaattg ggtttttg cattactttt + humanchr1/239220001-239550000 174456 TTAGGTTAGTGCAAAAGTAATTGTGGTTTTTGTCATTACTTTT 174498 + 589************************************9975 PP + +>> MADE1 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon + score bias Evalue hmmfrom hmm to alifrom ali to envfrom env to mod len acc + ------ ----- --------- ------- ------- --------- --------- --------- --------- --------- ---- + ! 28.9 6.0 1.2e-07 1 77 [. 302466 302390 .. 302466 302387 .. 80 0.74 + + Alignment: + score: 28.9 bits + xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx................xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF + MADE1 1 ttaggttggtgcaaaagtaattgcggtttttgccattactttt................aatggcaaaaaccgcaattacttttgcaccaacc 77 + ttag ttggtg aaaag cattactttt aatggcaaaaacc caatt ttttgcacc acc + humanchr1/239220001-239550000 302466 TTAGATTGGTGTAAAAG----------------CATTACTTTTaaaagcaattaaaagcAATGGCAAAAACCACAATTGATTTTGCACCGACC 302390 + 68999999999999998................5666777776222222222222222268****************************9998 PP + +>> MADE1 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon + score bias Evalue hmmfrom hmm to alifrom ali to envfrom env to mod len acc + ------ ----- --------- ------- ------- --------- --------- --------- --------- --------- ---- + ! 23.3 7.0 7.2e-06 43 80 .] 174493 174456 .. 174513 174456 .. 80 0.91 + + Alignment: + score: 23.3 bits + xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF + MADE1 43 taatggcaaaaaccgcaattacttttgcaccaacctaa 80 + taatg caaaaacc caattacttttgcac aacctaa + humanchr1/239220001-239550000 174493 TAATGACAAAAACCACAATTACTTTTGCACTAACCTAA 174456 + 689********************************985 PP + +>> MADE1 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon + score bias Evalue hmmfrom hmm to alifrom ali to envfrom env to mod len acc + ------ ----- --------- ------- ------- --------- --------- --------- --------- --------- ---- + ? 6.3 7.0 1.4 41 72 .. 304073 304104 .. 304053 304109 .. 80 0.85 + + Alignment: + score: 6.3 bits + xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF + MADE1 41 tttaatggcaaaaaccgcaattacttttgcac 72 + tt a tgg aaaaa ca tta ttttgca + humanchr1/239220001-239550000 304073 TTAAGTGGGAAAAAATACACTTATTTTTGCAT 304104 + 455779************************86 PP + + + +Internal pipeline statistics summary: +------------------------------------- +Query sequence(s): 1 (660000 residues searched) +Target model(s): 1 (80 nodes) +Residues passing SSV filter: 61794 (0.0936); expected (0.02) +Residues passing bias filter: 46199 (0.07); expected (0.02) +Residues passing Vit filter: 2752 (0.00417); expected (0.001) +Residues passing Fwd filter: 2526 (0.00383); expected (1e-05) +Total number of hits: 5 (0.000405) +# CPU time: 0.02u 0.00s 00:00:00.02 Elapsed: 00:00:00.02 +# Mc/sec: 2640.00 +// +[ok] diff -r 25fa92aa34c0 -r 1a83249ddfff test-data/MADE1.nhmmscan_out.aliscoresout --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/MADE1.nhmmscan_out.aliscoresout Sat Apr 07 03:51:55 2018 -0400 @@ -0,0 +1,17 @@ +# target name accession query name accession hmmfrom hmm to alifrom ali to envfrom env to modlen strand E-value score bias description of target +#------------------- ---------- -------------------- ---------- ------- ------- ------- ------- ------- ------- ------- ------ --------- ------ ----- --------------------- +MADE1 DF0000629.2 humanchr1/239220001-239550000 - 4 80 302390 302466 302387 302466 80 + 1.2e-10 38.6 7.4 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +MADE1 DF0000629.2 humanchr1/239220001-239550000 - 1 43 174456 174498 174456 174518 80 + 7.8e-08 29.6 8.3 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +MADE1 DF0000629.2 humanchr1/239220001-239550000 - 1 77 302466 302390 302466 302387 80 - 1.2e-07 28.9 6.0 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +MADE1 DF0000629.2 humanchr1/239220001-239550000 - 43 80 174493 174456 174513 174456 80 - 7.2e-06 23.3 7.0 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +MADE1 DF0000629.2 humanchr1/239220001-239550000 - 41 72 304073 304104 304053 304109 80 + 1.4 6.3 7.0 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +# +# Program: hmmscan +# Version: 3.1b2 (February 2015) +# Pipeline mode: SCAN +# Query file: /tmp/tmpc_c3amjg/files/000/dataset_2.dat +# Target file: /tmp/tmpc_c3amjg/files/000/dataset_1.dat +# Option settings: nhmmscan --tblout /tmp/tmpc_c3amjg/files/000/dataset_4.dat --dfamtblout /tmp/tmpc_c3amjg/files/000/dataset_5.dat --notextw -E 10.0 --F1 0.02 --F2 0.001 --F3 1e-05 --seed 4 --B1 110 --B2 240 --B3 1000 --cpu 1 /tmp/tmpc_c3amjg/files/000/dataset_1.dat /tmp/tmpc_c3amjg/files/000/dataset_2.dat +# Current dir: /tmp/tmpc_c3amjg/job_working_directory/000/3/working +# Date: Wed Apr 4 16:58:45 2018 +# [ok] diff -r 25fa92aa34c0 -r 1a83249ddfff test-data/MADE1.nhmmscan_out.dfamtblout --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/MADE1.nhmmscan_out.dfamtblout Sat Apr 07 03:51:55 2018 -0400 @@ -0,0 +1,10 @@ +# hit scores +# ---------- +# +# target name acc query name bits e-value bias hmm-st hmm-en strand ali-st ali-en env-st env-en modlen description of target +# ------------------- ------------------- ------------------- ------ --------- ----- ------- ------- ------ ------- ------- ------- ------- ------- --------------------- +MADE1 DF0000629.2 humanchr1/239220001-239550000 38.6 1.2e-10 7.4 4 80 + 302390 302466 302387 302466 80 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +MADE1 DF0000629.2 humanchr1/239220001-239550000 29.6 7.8e-08 8.3 1 43 + 174456 174498 174456 174518 80 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +MADE1 DF0000629.2 humanchr1/239220001-239550000 28.9 1.2e-07 6.0 1 77 - 302466 302390 302466 302387 80 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +MADE1 DF0000629.2 humanchr1/239220001-239550000 23.3 7.2e-06 7.0 43 80 - 174493 174456 174513 174456 80 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +MADE1 DF0000629.2 humanchr1/239220001-239550000 6.3 1.4 7.0 41 72 + 304073 304104 304053 304109 80 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon diff -r 25fa92aa34c0 -r 1a83249ddfff test-data/MADE1.nhmmscan_out.tblout --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/MADE1.nhmmscan_out.tblout Sat Apr 07 03:51:55 2018 -0400 @@ -0,0 +1,17 @@ +# target name accession query name accession hmmfrom hmm to alifrom ali to envfrom env to modlen strand E-value score bias description of target +#------------------- ---------- -------------------- ---------- ------- ------- ------- ------- ------- ------- ------- ------ --------- ------ ----- --------------------- +MADE1 DF0000629.2 humanchr1/239220001-239550000 - 4 80 302390 302466 302387 302466 80 + 1.2e-10 38.6 7.4 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +MADE1 DF0000629.2 humanchr1/239220001-239550000 - 1 43 174456 174498 174456 174518 80 + 7.8e-08 29.6 8.3 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +MADE1 DF0000629.2 humanchr1/239220001-239550000 - 1 77 302466 302390 302466 302387 80 - 1.2e-07 28.9 6.0 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +MADE1 DF0000629.2 humanchr1/239220001-239550000 - 43 80 174493 174456 174513 174456 80 - 7.2e-06 23.3 7.0 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +MADE1 DF0000629.2 humanchr1/239220001-239550000 - 41 72 304073 304104 304053 304109 80 + 1.4 6.3 7.0 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon +# +# Program: hmmscan +# Version: 3.1b2 (February 2015) +# Pipeline mode: SCAN +# Query file: /tmp/tmpc_c3amjg/files/000/dataset_2.dat +# Target file: /tmp/tmpc_c3amjg/files/000/dataset_1.dat +# Option settings: nhmmscan --tblout /tmp/tmpc_c3amjg/files/000/dataset_4.dat --dfamtblout /tmp/tmpc_c3amjg/files/000/dataset_5.dat --notextw -E 10.0 --F1 0.02 --F2 0.001 --F3 1e-05 --seed 4 --B1 110 --B2 240 --B3 1000 --cpu 1 /tmp/tmpc_c3amjg/files/000/dataset_1.dat /tmp/tmpc_c3amjg/files/000/dataset_2.dat +# Current dir: /tmp/tmpc_c3amjg/job_working_directory/000/3/working +# Date: Wed Apr 4 16:58:45 2018 +# [ok] diff -r 25fa92aa34c0 -r 1a83249ddfff test-data/nhmmer.out.dfamtblout --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/nhmmer.out.dfamtblout Sat Apr 07 03:51:55 2018 -0400 @@ -0,0 +1,17 @@ +# target name accession query name accession hmmfrom hmm to alifrom ali to envfrom env to sq len strand E-value score bias description of target +# ------------------- ---------- -------------------- ---------- ------- ------- ------- ------- ------- ------- ------- ------ --------- ------ ----- --------------------- +humanchr1/239220001-239550000 - MADE1 DF0000629.2 4 80 302390 302466 302387 302466 330000 + 1.2e-10 38.6 7.4 - +humanchr1/239220001-239550000 - MADE1 DF0000629.2 1 43 174456 174498 174456 174518 330000 + 7.8e-08 29.6 8.3 - +humanchr1/239220001-239550000 - MADE1 DF0000629.2 1 77 302466 302390 302466 302387 330000 - 1.2e-07 28.9 6.0 - +humanchr1/239220001-239550000 - MADE1 DF0000629.2 43 80 174493 174456 174513 174456 330000 - 7.2e-06 23.3 7.0 - +humanchr1/239220001-239550000 - MADE1 DF0000629.2 41 72 304073 304104 304053 304109 330000 + 1.4 6.3 7.0 - +# +# Program: nhmmer +# Version: 3.1b2 (February 2015) +# Pipeline mode: SEARCH +# Query file: /tmp/tmpnpsbd_6p/files/000/dataset_1.dat +# Target file: /tmp/tmpnpsbd_6p/files/000/dataset_2.dat +# Option settings: nhmmer --tblout /tmp/tmpnpsbd_6p/files/000/dataset_4.dat --dfamtblout None --notextw -E 10.0 --F1 0.02 --F2 0.001 --F3 1e-05 --dna --seed 4 --cpu 1 /tmp/tmpnpsbd_6p/files/000/dataset_1.dat /tmp/tmpnpsbd_6p/files/000/dataset_2.dat +# Current dir: /tmp/tmpnpsbd_6p/job_working_directory/000/3/working +# Date: Fri Apr 6 16:39:46 2018 +# [ok] diff -r 25fa92aa34c0 -r 1a83249ddfff test-data/nhmmer.out.tblout --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/nhmmer.out.tblout Sat Apr 07 03:51:55 2018 -0400 @@ -0,0 +1,17 @@ +# target name accession query name accession hmmfrom hmm to alifrom ali to envfrom env to sq len strand E-value score bias description of target +# ------------------- ---------- -------------------- ---------- ------- ------- ------- ------- ------- ------- ------- ------ --------- ------ ----- --------------------- +humanchr1/239220001-239550000 - MADE1 DF0000629.2 4 80 302390 302466 302387 302466 330000 + 1.2e-10 38.6 7.4 - +humanchr1/239220001-239550000 - MADE1 DF0000629.2 1 43 174456 174498 174456 174518 330000 + 7.8e-08 29.6 8.3 - +humanchr1/239220001-239550000 - MADE1 DF0000629.2 1 77 302466 302390 302466 302387 330000 - 1.2e-07 28.9 6.0 - +humanchr1/239220001-239550000 - MADE1 DF0000629.2 43 80 174493 174456 174513 174456 330000 - 7.2e-06 23.3 7.0 - +humanchr1/239220001-239550000 - MADE1 DF0000629.2 41 72 304073 304104 304053 304109 330000 + 1.4 6.3 7.0 - +# +# Program: nhmmer +# Version: 3.1b2 (February 2015) +# Pipeline mode: SEARCH +# Query file: /tmp/tmpnpsbd_6p/files/000/dataset_1.dat +# Target file: /tmp/tmpnpsbd_6p/files/000/dataset_2.dat +# Option settings: nhmmer --tblout /tmp/tmpnpsbd_6p/files/000/dataset_4.dat --dfamtblout None --notextw -E 10.0 --F1 0.02 --F2 0.001 --F3 1e-05 --dna --seed 4 --cpu 1 /tmp/tmpnpsbd_6p/files/000/dataset_1.dat /tmp/tmpnpsbd_6p/files/000/dataset_2.dat +# Current dir: /tmp/tmpnpsbd_6p/job_working_directory/000/3/working +# Date: Fri Apr 6 16:39:46 2018 +# [ok]