/tools/filters/axt_to_lav.xml

https://bitbucket.org/cistrome/cistrome-harvard/ · XML · 94 lines · 71 code · 23 blank · 0 comment · 0 complexity · 3669f5cd22f008b55af89bd1a3017ffb MD5 · raw file

  1. <tool id="axt_to_lav_1" name="AXT to LAV">
  2. <description>Converts an AXT formatted file to LAV format</description>
  3. <command interpreter="python">axt_to_lav.py /galaxy/data/$dbkey_1/seq/%s.nib:$dbkey_1:${GALAXY_DATA_INDEX_DIR}/shared/ucsc/chrom/${dbkey_1}.len /galaxy/data/$dbkey_2/seq/%s.nib:$dbkey_2:${GALAXY_DATA_INDEX_DIR}/shared/ucsc/chrom/${dbkey_2}.len $align_input $lav_file $seq_file1 $seq_file2</command>
  4. <inputs>
  5. <param name="align_input" type="data" format="axt" label="Alignment File" optional="False"/>
  6. <param name="dbkey_1" type="genomebuild" label="Genome"/>
  7. <param name="dbkey_2" type="genomebuild" label="Genome"/>
  8. </inputs>
  9. <outputs>
  10. <data name="lav_file" format="lav"/>
  11. <data name="seq_file1" format="fasta" parent="lav_file"/>
  12. <data name="seq_file2" format="fasta" parent="lav_file"/>
  13. </outputs>
  14. <help>
  15. .. class:: warningmark
  16. **IMPORTANT**: AXT formatted alignments will be phased out from Galaxy in the coming weeks. They will be replaced with pairwise MAF alignments, which are already available. To try pairwise MAF alignments use "Extract Pairwise MAF blocks" tool in *Fetch Sequences and Alignments* section.
  17. --------
  18. **Syntax**
  19. This tool converts an AXT formatted file to the LAV format.
  20. - **AXT format** The alignments are produced from Blastz, an alignment tool available from Webb Miller's lab at Penn State University. The lav format Blastz output, which does not include the sequence, was converted to AXT format with lavToAxt. Each alignment block in an AXT file contains three lines: a summary line and 2 sequence lines. Blocks are separated from one another by blank lines.
  21. - **LAV format** LAV is an alignment format developed by Webb Miller's group. It is the primary output format for BLASTZ.
  22. - **FASTA format** a text-based format for representing both nucleic and protein sequences, in which base pairs or proteins are represented using a single-letter code.
  23. - This format contains an one line header. It starts with a ">" symbol. The first word on this line is the name of the sequence. The rest of the line is a description of the sequence.
  24. - The remaining lines contain the sequence itself.
  25. - Blank lines in a FASTA file are ignored, and so are spaces or other gap symbols (dashes, underscores, periods) in a sequence.
  26. - Fasta files containing multiple sequences are just the same, with one sequence listed right after another. This format is accepted for many multiple sequence alignment programs.
  27. -----
  28. **Example**
  29. - AXT format::
  30. 0 chr19 3001012 3001075 chr11 70568380 70568443 - 3500
  31. TCAGCTCATAAATCACCTCCTGCCACAAGCCTGGCCTGGTCCCAGGAGAGTGTCCAGGCTCAGA
  32. TCTGTTCATAAACCACCTGCCATGACAAGCCTGGCCTGTTCCCAAGACAATGTCCAGGCTCAGA
  33. 1 chr19 3008279 3008357 chr11 70573976 70574054 - 3900
  34. CACAATCTTCACATTGAGATCCTGAGTTGCTGATCAGAATGGAAGGCTGAGCTAAGATGAGCGACGAGGCAATGTCACA
  35. CACAGTCTTCACATTGAGGTACCAAGTTGTGGATCAGAATGGAAAGCTAGGCTATGATGAGGGACAGTGCGCTGTCACA
  36. - Convert the above file to LAV format::
  37. #:lav
  38. s {
  39. &quot;/galaxy/data/hg16/seq/chr19.nib&quot; 1 63811651 0 1
  40. &quot;/galaxy/data/mm5/seq/chr11.nib-&quot; 1 121648857 0 1
  41. }
  42. h {
  43. &quot;> hg16.chr19&quot;
  44. &quot;> mm5.chr11 (reverse complement)&quot;
  45. }
  46. a {
  47. s 3500
  48. b 3001012 70568380
  49. e 3001075 70568443
  50. l 3001012 70568380 3001075 70568443 81
  51. }
  52. a {
  53. s 3900
  54. b 3008279 70573976
  55. e 3008357 70574054
  56. l 3008279 70573976 3008357 70574054 78
  57. }
  58. #:eof
  59. - With two files in the FASTA format::
  60. &gt;hg16.chr19_-_3001011_3001075
  61. TCAGCTCATAAATCACCTCCTGCCACAAGCCTGGCCTGGTCCCAGGAGAGTGTCCAGGCTCAGA
  62. &gt;hg16.chr19_-_3008278_3008357
  63. CACAATCTTCACATTGAGATCCTGAGTTGCTGATCAGAATGGAAGGCTGAGCTAAGATGAGCGACGAGGCAATGTCACA
  64. **and**::
  65. &gt;mm5.chr11_-_70568379_70568443
  66. TCTGTTCATAAACCACCTGCCATGACAAGCCTGGCCTGTTCCCAAGACAATGTCCAGGCTCAGA
  67. &gt;mm5.chr11_-_70573975_70574054
  68. CACAGTCTTCACATTGAGGTACCAAGTTGTGGATCAGAATGGAAAGCTAGGCTATGATGAGGGACAGTGCGCTGTCACA
  69. </help>
  70. <code file="axt_to_lav_code.py"/>
  71. </tool>