Mercurial > repos > iuc > mummer_dnadiff

comparison dnadiff.xml @ 0:deacc6b40cc0 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mummer4 commit 8133565adbfc012fa54b96449c2a18d044049107
author iuc
date Wed, 05 Dec 2018 02:38:01 -0500
parents
children 315c678562d2
comparison
equal deleted inserted replaced
-1:000000000000 0:deacc6b40cc0
1 <tool id="mummer_dnadiff" name="DNAdiff" version="@MUMMER_VERSION@">
2 <description>Evaluate similarities/differences between two sequences</description>
3 <macros>
4 <import>macros.xml</import>
5 </macros>
6 <expand macro="requirements" />
7 <command detect_errors="exit_code">
8 <![CDATA[
9 ln -s $reference_sequence reference.fa &&
10 ln -s $query_sequence query.fa &&
11 dnadiff
12 'reference.fa' 'query.fa'
13
14 ]]>
15 </command>
16 <inputs>
17 <param name="reference_sequence" type="data" format="fasta" label="Reference Sequence" help="FastA or multi-FastA" />
18 <param name="query_sequence" type="data" format="fasta" label="Query Sequence(s)" help="FastA or multi-FastA" />
19 <param name="report_only" type="select" label="Output only the general report file?" help="Select no to output all output files">
20 <option value="yes">YES</option>
21 <option value="no">NO</option>
22 </param>
23 </inputs>
24 <outputs>
25 <data name="report" format="txt" from_work_dir="out.report" label="${tool.name} on ${on_string}: report" />
26 <data name="delta" format="tabular" from_work_dir="out.delta" label="${tool.name} on ${on_string}: delta">
27 <filter> report_only == 'no'</filter>
28 </data>
29 <data name="1delta" format="tabular" from_work_dir="out.1delta" label="${tool.name} on ${on_string}: 1delta" >
30 <filter> report_only == 'no' </filter>
31 </data>
32 <data name="mdelta" format="tabular" from_work_dir="out.mdelta" label="${tool.name} on ${on_string}: mdelta" >
33 <filter> report_only == 'no' </filter>
34 </data>
35 <data name="1coords" format="tabular" from_work_dir="out.1coords" label="${tool.name} on ${on_string}: 1coords" >
36 <filter> report_only == 'no' </filter>
37 <actions>
38 <action name="column_names" type="metadata" default="[S1], [E1], [S2], [E2], [LEN 1], [LEN 2], [% IDY], [LEN R], [LEN Q], [COV R], [COV Q], [REF TAG], [QUERY TAG]" />
39 </actions>
40 </data>
41 <data name="mcoords" format="tabular" from_work_dir="out.mcoords" label="${tool.name} on ${on_string}: mcoords" >
42 <filter> report_only == 'no' </filter>
43 <actions>
44 <action name="column_names" type="metadata" default="[S1], [E1], [S2], [E2], [LEN 1], [LEN 2], [% IDY], [LEN R], [LEN Q], [COV R], [COV Q], [REF TAG], [QUERY TAG]" />
45 </actions>
46 </data>
47 <data name="snps" format="tabular" from_work_dir="out.snps" label="${tool.name} on ${on_string}: snps" >
48 <filter> report_only == 'no' </filter>
49 <actions>
50 <action name="column_names" type="metadata" default="[P1], [REF SUB], [QUERY SUB], [P2], [BUFF], [DIST], [LEN REF], [LEN QUERY], [REF FRAME], [QUERY FRAME], [REF TAG], [QUERY TAG]" />
51 </actions>
52 </data>
53 <data name="rdiff" format="tabular" from_work_dir="out.rdiff" label="${tool.name} on ${on_string}: rdiff" >
54 <filter> report_only == 'no' </filter>
55 <actions>
56 <action name="column_names" type="metadata" default="Seq ID, Feature Type, Feature Start, Feature End, Feature Length" />
57 </actions>
58 </data>
59 <data name="qdiff" format="tabular" from_work_dir="out.qdiff" label="${tool.name} on ${on_string}: qdiff" >
60 <filter> report_only == 'no' </filter>
61 <actions>
62 <action name="column_names" type="metadata" default="Seq ID, Feature Type, Feature Start, Feature End, Feature Length" />
63 </actions>
64 </data>
65 </outputs>
66 <tests>
67 <test>
68 <param name="reference_sequence" ftype="fasta" value="human_aqp3.fasta" />
69 <param name="query_sequence" ftype="fasta" value="mouse_aqp3.fasta"/>
70 <param name="input_type" value="sequence" />
71 <param name="report_only" value="no" />
72 <output name="report" ftype="txt" compare="diff" sort="true" lines_diff="2" value="report.txt"/>
73 <output name="delta" ftype="tabular" compare="diff" lines_diff="2" value="delta.txt"/>
74 <output name="1delta" ftype="tabular" compare="diff" lines_diff="2" value="1delta.txt"/>
75 <output name="mdelta" ftype="tabular" compare="diff" lines_diff="2" value="mdelta.txt"/>
76 <output name="1coords" ftype="tabular" compare="diff" value="1coords.txt"/>
77 <output name="mcoords" ftype="tabular" compare="diff" value="mcoords.txt"/>
78 <output name="snps" ftype="tabular" compare="diff" value="snps.txt"/>
79 <output name="rdiff" ftype="tabular" compare="diff" value="rdiff.txt"/>
80 <output name="qdiff" ftype="tabular" compare="diff" value="qdiff.txt"/>
81 </test>
82 </tests>
83 <help><![CDATA[
84 This script is a wrapper around nucmer that builds an alignment using default parameters, and runs many of nucmer's helper scripts to process the output and report alignment statistics, SNPs, breakpoints, etc. It is designed for evaluating the sequence and structural similarity of two highly similar sequence sets. E.g. comparing two different assemblies of the same organism, or comparing two strains of the same species.
85
86 **Output files:**
87 * report: Summary of alignments, differences and SNPs
88 * delta: Standard nucmer alignment output
89 * 1delta: 1-to-1 alignment from delta-filter -1
90 * mdelta: M-to-M alignment from delta-filter -m
91 * 1coords: 1-to-1 coordinates from show-coords -THrcl .1delta
92 * mcoords: M-to-M coordinates from show-coords -THrcl .mdelta
93 * snps: SNPs from show-snps -rlTHC .1delta
94 * rdiff: Classified ref breakpoints from show-diff -rH .mdelta
95 * qdiff: Classified qry breakpoints from show-diff -qH .mdelta
96 ]]></help>
97 <expand macro="citation" />
98 </tool>