Mercurial > repos > abossers > mummer_toolsuite

annotate MUMmer/mummer_clustering.xml @ 0:59f302448cf6

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
author abossers
date Tue, 07 Jun 2011 17:22:27 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
1 <tool id="mummer_clustering" name="MUMmer Clustering" version="0.9.alx" force_history_refresh="True">
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
2 <description>: order sequence matches in clusters</description>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
3 <command>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
4 <!-- update this path to the installed location -->
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
5 /opt/MUMmer/MUMmer/$tool.cmd
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
6 #if $tool.cmd=="gaps":
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
7 $in_reference
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
8 #if $tool.gaps_r=="yes":
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
9 -r
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
10 #end if
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
11 #end if
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
12 #if $tool.cmd=="mgaps":
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
13 #if $tool.cmd_C=="yes":
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
14 -C
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
15 #end if
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
16 -d $tool.cmd_d
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
17 #if $tool.cmd_e=="yes":
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
18 -e
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
19 #end if
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
20 -f $tool.cmd_f
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
21 -l $tool.cmd_l
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
22 -s $tool.cmd_s
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
23 #end if
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
24 &lt; $tool.in_match_list
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
25 &gt; $out_tool
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
26
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
27 </command>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
28 <inputs>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
29 <conditional name="tool">
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
30 <param name="cmd" type="select" label="MUMmer maximal matching" help="Algorithms are run with default parameters (none). For specific args see help below" >
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
31 <option value="gaps" selected="true">gaps</option>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
32 <option value="mgaps">mgaps</option>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
33 </param>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
34 <when value="gaps">
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
35 <param name="in_reference" type="data" format="fasta" label="Reference FastA file" />
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
36 <param name="gaps_r" type="select" label="Use reversed [-r]" >
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
37 <option value="no" selected="true">No</option>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
38 <option value="yes">Yes</option>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
39 </param>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
40 <param name="in_match_list" type="data" format="text" label="MUMmer match list" help="See help for more details" />
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
41 </when>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
42 <when value="mgaps">
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
43 <param name="in_match_list" type="data" format="text" label="MUMmer match list" help="See help for more details" />
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
44 <param name="cmd_C" type="select" label="Check input header labels have reversed keyword [-C]" >
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
45 <option value="no" selected="true">No</option>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
46 <option value="yes">Yes</option>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
47 </param>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
48 <param name="cmd_d" type="integer" size="5" value="5" label="Max fixed diagonal difference [-d]" />
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
49 <param name="cmd_e" type="select" label="Use extent of cluster [-e]" >
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
50 <option value="no" selected="true">No</option>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
51 <option value="yes">Yes</option>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
52 </param>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
53 <param name="cmd_f" type="float" size="5" value="0.05" label="Max fraction separation for diagonal difference [-f]" />
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
54 <param name="cmd_l" type="integer" size="5" value="200" label="Min cluster length [-l]" />
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
55 <param name="cmd_s" type="integer" size="5" value="1000" label="Max separation adjecent matches in cluster [-s]" />
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
56 </when>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
57 </conditional>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
58 </inputs>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
59 <outputs>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
60 <data name="out_tool" format="text" label="Clustering output" />
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
61 </outputs>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
62 <requirements>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
63 <requirement type="binary">gaps</requirement>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
64 <requirement type="binary">mgaps</requirement>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
65 </requirements>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
66 <tests>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
67 <test>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
68 </test>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
69 </tests>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
70 <help>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
71 |
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
72
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
73
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
74 **Reference**
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
75 =============
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
76
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
77 - **MUMmer clustering Galaxy tool wrapper:** Alex Bossers, CVI of Wageningen UR, The Netherlands.
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
78
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
79 - **MUMmer suite v3.22:** http://mummer.sourceforge.net
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
80
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
81 - **MUMmer tutorials:** http://mummer.sourceforge.net/examples/
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
82
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
83 If you found these tools/wrappers usefull in your research, please acknowledge our work. If you improve
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
84 or modify the wrappers please add instead of substitute yourself into the acknowlegement section :)
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
85
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
86
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
87 **MUMmer Clustering**
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
88 =====================
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
89
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
90 MUMmer's clustering algorithms attempt to order small individual matches into larger match clusters
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
91 in order to make the output of mummer more intelligible. A dot plot makes it easy to spot alignment
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
92 regions from a match list, however when examining the data without graphic aids, it is very difficult
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
93 to draw any reasonable conclusions from the simple flat file list of matches. Clustering the matches
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
94 together into larger groups of neighboring matches makes this process much easier by ordering the
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
95 data and removing spurious matches.
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
96
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
97
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
98 Gaps
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
99 ----
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
100
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
101 *gaps* is the primary clustering algorithm for run-mummer1, and although classified as a "clustering"
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
102 step, gaps is more of a sorting routine. It implements the LIS (longest increasing subset) algorithm
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
103 to extract the longest consistent set of matches between two sequences, and generates a single
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
104 cluster that represents the best "straight-line" arrangement of matches between the sequences. By
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
105 straight-line, we mean no rearrangements or inversions, just a simple path of agreeing matches
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
106 between the two sequences. This limits the usability of this program to the alignment of genomes
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
107 that are very similar and with no large scale mutations. *gaps* is best suited for the comparison of
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
108 near identical sequences with the goal of finding minor mutations like SNPs and small indels.
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
109
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
110 Input can be filtered mummer output. The strange syntax is a result of a legacy issue described in
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
111 the Known problems (manual) section, and requires the header be stripped from the mummer output. In
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
112 addition, gaps is only designed to handle a single reference and a single query sequence, thus the
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
113 preceding mummer run must also follow this constraint. The -r is optional and designates the incoming
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
114 matches as reverse complement matches which must reference the reverse complement of the sequence,
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
115 therefore forcing mummer to be run without the -c option.
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
116
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
117 Reference: http://mummer.sourceforge.net/manual/#gaps
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
118
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
119 **Output:**
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
120 ::
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
121
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
122 > /home/aphillip/data/GHP.1con Consistent matches
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
123 183 17 22 none - -
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
124 238 72 108 none 33 33
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
125 347 181 92 none 1 1
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
126 458 292 50 none 19 19
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
127 705 539 44 none 1 1
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
128 750 584 38 none 1 1
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
129 807 641 23 -16 0 4
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
130 (output continues ...)
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
131 > Wrap around
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
132 334398 329917 47 none - 225
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
133 334446 329965 62 none 1 1
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
134 334539 330058 20 none 31 31
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
135 334560 330079 92 none 1 1
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
136 334653 330172 77 none 1 1
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
137 334740 330259 41 none 10 10
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
138 (output continues ...)
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
139 > /home/aphillip/data/GHP.1con Other matches
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
140 1317231 4891 21 none - -
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
141 1317275 4927 21 none - -
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
142 1317804 5399 25 none 508 451
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
143 947580 5436 36 none - -
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
144 23406 5518 34 none - -
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
145 333079 6592 32 none - -
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
146 (output continues ...)
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
147
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
148 Where the first line is the location of the reference file, and the first three columns are the same
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
149 as the three column match format described in the mummer section. The final three columns are the
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
150 overlap between this match and the previous match, the gap between the start of this match and the
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
151 end of the previous match in the reference, and the gap between the start of this match and the end
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
152 of the previous match in the query respectively.
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
153
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
154
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
155 mgaps
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
156 -----
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
157
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
158 *mgaps* was introduced into the MUMmer pipeline in an effort to better handle large-scale
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
159 rearrangements and duplications. Unlike gaps, mgaps is a full clustering algorithm that is capable
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
160 of generating multiple groups of consistently ordered matches. Clustering is controlled by a set of
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
161 command-line parameters that adjust the minimum cluster size, maximum gap between matches, etc. Only
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
162 matches that were included in clusters will appear in the output, so by adjusting the command-line
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
163 parameters it is possible to filter out many of the spurious matches, thus leaving only the larger
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
164 areas of conservation between the input sequences. The major advantage of mgaps is its ability to
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
165 identify these "islands" of conservation. This frees the user from the single LIS restraints of the
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
166 gaps program and allows for the identification of large-scale rearrangements, duplications, gene
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
167 families and so on.
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
168
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
169 Gaps can fail to identify clusters because they were not consistent with the LIS. However, by using
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
170 mgaps, all regions of conservation can now been identified. The only fallback being the increased
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
171 complexity of the output, where you once had only one cluster for the whole comparison, you usually
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
172 now get more. Because of this, it can sometimes be difficult separating the repetitive clusters from
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
173 "correct" clusters, *making mgaps more suited for global alignments instead of localized error detection*.
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
174
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
175 Input can be raw mummer output. *mgaps* is only designed to handle a single reference and one or
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
176 more query sequences, thus the preceding mummer run must also follow this constraint. Please refer
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
177 to the run-mummer3 script (see online manual) for an example of how to use this program in an
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
178 alignment pipeline. Note that in order to cluster reverse complement matches, the reverse complement
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
179 matches must reference the reverse complement strand of the query sequence, therefore forcing mummer
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
180 to be run without the -c option. A rewrite of this algorithm to handle multiple reference sequences
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
181 and a better coordinate system (forward coordinates for reverse complement matches) is doubtful but
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
182 may eventually appear.
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
183
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
184 The -d option can be interpreted as the number of insertions allowed between two matches in the same
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
185 cluster, while the -f option is a fraction equal to (diagonal difference / match separation) where
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
186 a higher value will increase the indel tolerance. Minimum cluster length is the sum of the contained
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
187 matches unless the -e option is used. The best way to get a feel for what each parameter controls
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
188 is to cluster the same data set numerous times with different values and observe the resulting
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
189 differences. It can also be helpful to set these parameters to the size of the element you wish to
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
190 capture, i.e. set the minimum cluster size to say the smallest exon you expect and set the max gap
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
191 to the smallest intron you expect to obtain clusters that could represent single exons (depending
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
192 of course of the similarity of the two sequences).
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
193
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
194 Reference: http://mummer.sourceforge.net/manual/#mgaps
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
195
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
196 **Output format**
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
197
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
198 Output of *mgaps* shares much in common with the output of mummer and gaps, with a slightly different
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
199 header formatting than gaps to allow for multiple query sequences and multiple clusters. The output
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
200 of mgaps run on both forward and reverse complement matches is as follows:
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
201 ::
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
202
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
203 > ID41
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
204 > ID41 Reverse
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
205 5177399 1 232 none - -
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
206 5177632 234 6794 none 1 1
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
207 5184433 7035 24 none 7 7
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
208 5184468 7069 23 none 11 10
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
209 > ID42
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
210 10181 43 1521 none - -
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
211 > ID42 Reverse
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
212 4654536 17 36 none - -
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
213 4654578 57 298 none 6 4
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
214 4654877 356 226 none 1 1
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
215 #
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
216 4655139 845 28 none - -
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
217 4655178 884 694 none 11 11
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
218 4655873 1579 20 none 1 1
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
219 #
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
220 4850044 17 1492 none - -
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
221 4851537 1510 711 none 1 1
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
222 4852249 2222 42 none 1 1
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
223 (output continues ...)
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
224
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
225
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
226 Headers containing the ID for each query sequence are listed after the '>' characters, and a
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
227 following Reverse keyword identifies the reverse matches for that query sequence. Individual clusters
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
228 for each sequence are separated by a '#' character, and the six columns are exactly the same as the
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
229 gaps output (see the gaps section for more details).
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
230
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
231
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
232 |
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
233 |
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
234
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
235 </help>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
236 </tool>
59f302448cf6 Migrated tool version 1.0.0 from old tool shed archive to new tool shed repository
abossers
parents:
diff changeset
237