/tools/hyphy/hyphy_branch_lengths_wrapper.xml

https://bitbucket.org/cistrome/cistrome-harvard/ · XML · 95 lines · 77 code · 12 blank · 6 comment · 0 complexity · 9ef7c03b3585b8745355a9bcf48f3ad4 MD5 · raw file

  1. <?xml version="1.2.1"?>
  2. <tool name="Branch Lengths" id="hyphy_branch_lengths_wrapper1">
  3. <description>Estimation</description>
  4. <command interpreter="python">hyphy_branch_lengths_wrapper.py $input1 $out_file1 "$tree" "$model" "$base_freq" "Global" ${GALAXY_DATA_INDEX_DIR}</command>
  5. <inputs>
  6. <page>
  7. <param format="fasta" name="input1" type="data" label="Fasta file"/>
  8. <param name="tree" type="text" label="Tree Definition" size="20" help="For example: ((hg17,panTro1),(mm5,rn3),canFam1)"/>
  9. <param name="model" type="select" label="Substitution Model">
  10. <option value="000000">F81</option>
  11. <option value="010010">HKY85</option>
  12. <option value="012345">REV</option>
  13. </param>
  14. <!-- <param name="model_options" type="select" label="Model Options">
  15. <option value="Local">All model parameters are estimated independently for each branch</option>
  16. <option value="Global">Model parameters are shared by all branches, branch lengths are estimated independently</option>
  17. <option value="Global w/variation">Model parameters are shared by all branches, branch lengths come from a user-chosen distribution, whose parameters are estimated</option>
  18. <option value="Global w/variation+HM">Model parameters are shared by all branches, branch lengths come from a user-chosen distribution, whose parameters is estimated; rates at adjacent sites are correlated via a simple Hidden Markov model with an autocorrelation parameter lambda</option>
  19. </param> -->
  20. <param name="base_freq" type="select" label="Base Frequencies">
  21. <option value="Observed">Nucleotide frequencies collected from the data file will be used as equilibrium frequencies</option>
  22. <option value="Equal">Equal (.25) frequencies are used as equilibrium frequencies</option>
  23. </param>
  24. </page>
  25. </inputs>
  26. <outputs>
  27. <data name="out_file1" format="tabular" />
  28. </outputs>
  29. <tests>
  30. <test>
  31. <param name="input1" value="branchlength_in.fasta"/>
  32. <param name="tree" value="((hg17,panTro1),(mm5,rn3),canFam1)"/>
  33. <param name="model" value="012345"/>
  34. <param name="base_freq" value="Observed"/>
  35. <output name="out_file1" file="branchlength_out.tabular"/>
  36. </test>
  37. </tests>
  38. <help>
  39. This tool takes a single or multiple FASTA alignment file and estimates branch lengths using HYPHY_, a maximum likelihood analyses package.
  40. For the tree definition, you only need to specify the species build names. For example, you could use the tree *((hg17,panTro1),(mm5,rn3),canFam1)*, if your FASTA file looks like this::
  41. &gt;hg17.chr7(+):26907301-26907310|hg17_0
  42. GTGGGAGGT
  43. &gt;panTro1.chr6(+):28037319-28037328|panTro1_0
  44. GTGGGAGGT
  45. &gt;mm5.chr6(+):52104022-52104031|mm5_0
  46. GTGGGAGGT
  47. &gt;rn3.chr4(+):80734395-80734404|rn3_0
  48. GTGGGAGGT
  49. &gt;canFam1.chr14(+):42826409-42826418|canFam1_0
  50. GTGGGAGGT
  51. &gt;hg17.chr7(+):26907310-26907326|hg17_1
  52. AGTCAGAGTGTCTGAG
  53. &gt;panTro1.chr6(+):28037328-28037344|panTro1_1
  54. AGTCAGAGTGTCTGAG
  55. &gt;mm5.chr6(+):52104031-52104047|mm5_1
  56. AGTCAGAGTGTCTGAG
  57. &gt;rn3.chr4(+):80734404-80734420|rn3_1
  58. AGTCAGAGTATCTGAG
  59. &gt;canFam1.chr14(+):42826418-42826434|canFam1_1
  60. AGTCAGAGTGTCTGAG
  61. &gt;hg17.chr7(+):26907326-26907338|hg17_2
  62. GTAGAAGACCCC
  63. &gt;panTro1.chr6(+):28037344-28037356|panTro1_2
  64. GTAGAAGACCCC
  65. &gt;mm5.chr6(+):52104047-52104059|mm5_2
  66. GTAGACGATGCC
  67. &gt;rn3.chr4(+):80734420-80734432|rn3_2
  68. GTAGATGATGCG
  69. &gt;canFam1.chr14(+):42826434-42826446|canFam1_2
  70. GTAGAAGACCCC
  71. &gt;hg17.chr7(+):26907338-26907654|hg17_3
  72. GGGGAAGGAACGCAGGGCGAAGAGCTGGACTTCTCTGAGGAT---TCCTCGGCCTTCTCGT-----CGTTTCCTGG----CGGGGTGGCCGGAGAGATGGGCAAGAGACCCTCCTTCTCACGTTTCTTTTGCTTCATTCGGCGGTTCTGGAACCAGATCTTCACTTGGGTCTCGTTGAGCTGCAGGGATGCAGCGATCTCCACCCTGCGGGCGCGCGTCAGGTACTTGTTGAAGTGGAACTCCTTCTCCAGTTCCGTGAGCTGCTTGGTAGTGAAGTTGGTGCGCACCGCGTTGGGTTGACCCAGGTAGCCGTACTCTCCAACTTTCC
  73. &gt;panTro1.chr6(+):28037356-28037672|panTro1_3
  74. GGGGAAGGAACGCAGGGCGAAGAGCTGGACTTCTCTGAGGAT---TCCTCGGCCTTCTCGT-----CGTTTCCTGG----CGGGGTGGCCGGAGAGATGGGCAAGAGACCCTCCTTCTCACGTTTCTTTTGCTTCATTCGGCGGTTCTGGAACCAGATCTTCACTTGGGTCTCGTTGAGCTGCAGGGATGCAGCGATCTCCACCCTGCGGGCGCGCGTCAGGTACTTGTTGAAGTGGAACTCCTTCTCCAGTTCCGTGAGCTGCTTGGTAGTGAAGTTGGTGCGCACCGCGTTGGGTTGACCCAGGTAGCCGTACTCTCCAACTTTCC
  75. &gt;mm5.chr6(+):52104059-52104375|mm5_3
  76. GGAGAAGGGGCACTGGGCGAGGGGCTAGATTTCTCAGATGAT---TCTTCCGTTTTCTCAT-----CGCTGCCAGG----AGGAGTGGCAGGGGAGATGGGCAGGAGCCCCTCCTTCTCACGCTTCTTCTGCTTCATGCGGCGATTCTGGAACCAGATCTTCACCTGGGTCTCATTGAGCTGTAGGGACGCGGCAATCTCCACCCTGCGCGCTCGTGTAAGGTACTTGTTGAAGTGGAACTCCTTCTCCAGCTCTGTGAGCTGCTTGGTGGTGAAATTGGTGCGCACTGCGTTGGGTTGACCCACGTAGCCGTACTCTCCAACTTTCC
  77. &gt;rn3.chr4(+):80734432-80734748|rn3_3
  78. GGAGAAGGGGCGCTGGGCGAGGAGCTGGATTTCTCAGATGAT---TCTTCAGTTTTCTCAT-----CGCTTCCAGG----AGGGGTGGCGGGTGAAATGGGCAAGAGCCCCTCTTTCTCGCGCTTCTTCTGCTTCATGCGGCGATTCTGGAACCAGATCTTCACCTGGGTCTCATTGAGTTGCAGGGACGCGGCTATCTCCACCCTGCGGGCTCTTGTTAGGTACTTGTTGAAGTGGAACTCCTTCTCCAGCTCTGTGAGCTGCTTGGTGGTGAAGTTGGTGCGCACTGCGTTGGGTTGACCCACGTAGCCATACTCTCCAACTTTCC
  79. &gt;canFam1.chr14(+):42826446-42826762|canFam1_3
  80. GGAGACGGAATGCAGGGCGAGGAGCTGGATTTCTCTGAAGAT---TCCTCCGCCTTCTCCT-----CACTTCCTGG----CGGGGTGGCAGGGGAGATGGGCAAAAGGCCCTCTTTCTCTCGTTTCTTCTGCTTCATCCGGCGGTTCTGGAACCAGATCTTCACCTGGGTCTCGTTGAGCTGCAGGGATGCTGCGATCTCCACCCTGCGGGCGCGGGTCAGATACTTATTGAAGTGGAACTCCTTTTCCAGCTCGGTGAGCTGCTTGGTGGTGAAGTTGGTACGCACTGCATTCGGTTGACCCACGTAGCCGTACTCTCCAACTTTCC
  81. .. _HYPHY: http://www.hyphy.org
  82. </help>
  83. </tool>