|
2
|
1 <tool name="mashmap" id="mashmap" version="3.1.3" profile="22.05">
|
|
|
2 <!--Source in git at: https://github.com/fubar2/galaxy_tf_overlay-->
|
|
4
|
3 <!--Created by toolfactory@galaxy.org at 23/02/2024 20:28:46 using the Galaxy Tool Factory.-->
|
|
2
|
4 <description>Fast local alignment boundaries</description>
|
|
|
5 <requirements>
|
|
|
6 <requirement version="3.1.3" type="package">mashmap</requirement>
|
|
|
7 </requirements>
|
|
|
8 <version_command><![CDATA[echo "3.1.3"]]></version_command>
|
|
|
9 <command><![CDATA[bash '$runme']]></command>
|
|
|
10 <configfiles>
|
|
4
|
11 <configfile name="runme"><![CDATA[mashmap --pi '$perc_identity' -s '$seqLength' -f '$filtermode' $dense \
|
|
2
|
12 #if int($sketchSize) > 0:
|
|
|
13 -J '$sketchSize' \
|
|
|
14 #end if
|
|
|
15 #if len($reflist) == 1:
|
|
3
|
16 -r $reflist -q '$query'
|
|
2
|
17 #else
|
|
|
18 #for i, mash in enumerate($reflist):
|
|
|
19 #if i == 0:
|
|
3
|
20 echo $mash > 'reflist'
|
|
2
|
21 #else:
|
|
3
|
22 echo $mash >> 'reflist'
|
|
2
|
23 #end if
|
|
|
24 #end for
|
|
|
25 --rl 'reflist' -q '$query'
|
|
|
26 #end if
|
|
3
|
27 cp 'mashmap.out' $mashout]]></configfile>
|
|
2
|
28 </configfiles>
|
|
|
29 <inputs>
|
|
|
30 <param name="query" type="data" optional="false" label="Query sequences (as fasta) to mash against the references supplied below" help="" format="fasta" multiple="false"/>
|
|
|
31 <param name="reflist" type="data" optional="false" label="Reference or references to mash the query sequences on" help="Choose one or more reference sequences to mash the query sequences against." format="fasta" multiple="true"/>
|
|
|
32 <param name="perc_identity" type="float" value="85.0" label="Identity threshold" help="By default, it is set to 85, implying mappings with 85 or more identity should be reported. For example, it can be set to 80to account for more noisy long-read datasets or 95 for mapping human genome assembly to human reference."/>
|
|
|
33 <param name="seqLength" type="integer" value="5000" label="Minimum segment length" help="Default is 5,000 bp. Sequences below this length are ignored. Mashmap provides guarantees on reporting local alignments of length twice this value."/>
|
|
|
34 <param name="sketchSize" type="integer" value="0" label="Sketch size - leave 0 for automatic setting based" help="This parameter sets the seed density of the winnowing scheme, gauranteeing that the minhash will be calculated from a sample of sketchSize k-mers for each segment. It is set automatically based on --pi but can be manually set as well."/>
|
|
4
|
35 <param name="dense" type="select" label="Dense sketching" help="This flag will increase the seed density substantially, resulting in a density of roughly 0.02 * (1 + (1 - pi) / .05) where pi is the perc_identity threshold. This leads to longer runtimes and higher RAM usage, but significantly more accurate estimates of ANI.">
|
|
|
36 <option value="">No dense sketching</option>
|
|
|
37 <option value="--dense">Dense sketching</option>
|
|
3
|
38 </param>
|
|
2
|
39 <param name="filtermode" type="select" label="Filter mode" help="Mashmap implements a plane-sweep based algorithm to perform the alignment filtering. Similar to delta-filter in nucmer, different filtering options are provided that are suitable for long read or assembly mapping. Option -f map is suitable for reporting the best mappings for long reads, whereas -f one-to-one is suitable for reporting orthologous mappings among all computed assembly to genome mappings.">
|
|
|
40 <option value="map">map - best mapping for long reads</option>
|
|
|
41 <option value="one-to-one">one-to-one - best for mapping orthologous reads</option>
|
|
|
42 <option value="none">None</option>
|
|
|
43 </param>
|
|
|
44 </inputs>
|
|
|
45 <outputs>
|
|
|
46 <data name="mashout" format="paf" label="mashmap on $query.element_identifier" hidden="false"/>
|
|
|
47 </outputs>
|
|
|
48 <tests>
|
|
|
49 <test>
|
|
|
50 <output name="mashout" value="mashout_sample" compare="diff" lines_diff="0"/>
|
|
|
51 <param name="query" value="query_sample"/>
|
|
|
52 <param name="reflist" value="reflist_sample"/>
|
|
|
53 <param name="perc_identity" value="85.0"/>
|
|
|
54 <param name="seqLength" value="5000"/>
|
|
|
55 <param name="sketchSize" value="0"/>
|
|
4
|
56 <param name="dense" value=""/>
|
|
2
|
57 <param name="filtermode" value="map"/>
|
|
|
58 </test>
|
|
|
59 </tests>
|
|
|
60 <help><![CDATA[
|
|
|
61 *MashMap* implements a fast and approximate algorithm for computing local alignment boundaries between long DNA sequences. It can be useful for mapping genome assembly or long reads (PacBio/ONT) to reference genome(s). Given a minimum alignment length and an identity threshold for the desired local alignments,
|
|
|
62
|
|
|
63 Mashmap computes alignment boundaries and identity estimates using k-mers. It does not compute the alignments explicitly, but rather estimates an unbiased k-mer based Jaccard similarity using a combination of minmers (a novel winnowing scheme) and MinHash. This is then converted to an estimate of sequence identity using the Mash distance. An appropriate k-mer sampling rate is automatically determined using the given minimum local alignment length and identity thresholds.
|
|
|
64
|
|
|
65 As an example, Mashmap can map a human genome assembly to the human reference genome in about one minute total execution time and < 4 GB memory using just 8 CPU threads, achieving more than an order of magnitude improvement in both runtime and memory over alternative methods. We describe the algorithms associated with Mashmap, and report on speed, scalability, and accuracy of the software in the publications listed below. Unlike traditional mappers, MashMap does not compute exact sequence alignments. In future, we plan to add an optional alignment support to generate base-to-base alignments.
|
|
|
66
|
|
|
67 Map set of query sequences against a reference genome:
|
|
|
68
|
|
|
69 mashmap -r reference.fna -q query.fa
|
|
|
70
|
|
|
71 The output is space-delimited with each line consisting of query name, length, 0-based start, end, strand, target name, length, start, end and mapping nucleotide identity.
|
|
|
72
|
|
|
73 Map set of query seqences against a list of reference genomes:
|
|
|
74
|
|
|
75 mashmap --rl referenceList.txt -q query.fa
|
|
|
76
|
|
|
77 File 'referenceList.txt' containing the list of reference genomes should contain path to the reference genomes, one per line.
|
|
|
78
|
|
|
79 Source code: https://github.com/marbl/MashMap
|
|
|
80 ]]></help>
|
|
|
81 <citations>
|
|
|
82 <citation type="doi">10.1093/bioinformatics/btad512</citation>
|
|
|
83 <citation type="doi">10.1093/bioinformatics/bts573</citation>
|
|
|
84 </citations>
|
|
|
85 </tool>
|
|
|
86
|
