comparison mashmap.xml @ 6:27df186d5446 draft

planemo upload
author fubar
date Sat, 24 Feb 2024 04:11:27 +0000
parents 10e4181a6443
children 53f601fb8664
comparison
equal deleted inserted replaced
5:10e4181a6443 6:27df186d5446
1 <tool name="mashmap" id="mashmap" version="3.1.3" profile="22.05"> 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--> 2 <!--Source in git at: https://github.com/fubar2/galaxy_tf_overlay-->
3 <!--Created by toolfactory@galaxy.org at 23/02/2024 20:50:55 using the Galaxy Tool Factory.--> 3 <!--Created by toolfactory@galaxy.org at 23/02/2024 21:34:16 using the Galaxy Tool Factory.-->
4 <description>Fast local alignment boundaries</description> 4 <description>Fast local alignment boundaries</description>
5 <requirements> 5 <requirements>
6 <requirement version="3.1.3" type="package">mashmap</requirement> 6 <requirement version="3.1.3" type="package">mashmap</requirement>
7 </requirements> 7 </requirements>
8 <version_command><![CDATA[echo "3.1.3"]]></version_command> 8 <version_command><![CDATA[echo "3.1.3"]]></version_command>
9 <command><![CDATA[bash '$runme']]></command> 9 <command><![CDATA[bash '$runme']]></command>
10 <configfiles> 10 <configfiles>
11 <configfile name="runme"><![CDATA[#if len($reflist) > 1: 11 <configfile name="runme"><![CDATA[#if len($reflist) > 1:
12 #for i, mash in enumerate($reflist): 12 #for i, mash in enumerate($reflist):
13 #if i == 0: 13 #if i == 0:
14 echo $mash > 'reflist' && 14 echo '$mash' > 'reflist' &&
15 #else: 15 #else:
16 echo $mash >> 'reflist' && 16 echo '$mash' >> 'reflist' &&
17 #end if 17 #end if
18 #end for 18 #end for
19 #end if 19 #end if
20 mashmap --pi '$perc_identity' -s '$seqLength' -f '$filtermode' $dense \ 20 mashmap --pi '$perc_identity' -s '$seqLength' -f '$filtermode' $dense \
21 #if int($sketchSize) > 0: 21 #if int($sketchSize) > 0:
22 -J '$sketchSize' \ 22 -J '$sketchSize' \
23 #end if 23 #end if
24 #if len($reflist) == 1: 24 #if len($reflist) == 1:
25 -r $reflist -q '$query' && 25 -r '$reflist' -q '$query' &&
26 #else 26 #else
27 --rl 'reflist' -q '$query' && 27 --rl 'reflist' -q '$query' &&
28 #end if 28 #end if
29 cp 'mashmap.out' $mashout]]></configfile> 29 cp 'mashmap.out' '$mashout']]></configfile>
30 </configfiles> 30 </configfiles>
31 <inputs> 31 <inputs>
32 <param name="query" type="data" optional="false" label="Query sequences (as fasta) to mash against the references supplied below" help="" format="fasta" multiple="false"/> 32 <param name="query" type="data" optional="false" label="Query sequences (as fasta) to mash against the references supplied below" help="" format="fasta" multiple="false"/>
33 <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"/> 33 <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"/>
34 <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."/> 34 <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."/>
68 68
69 Map set of query sequences against a reference genome: 69 Map set of query sequences against a reference genome:
70 70
71 mashmap -r reference.fna -q query.fa 71 mashmap -r reference.fna -q query.fa
72 72
73 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. 73 The output is a paf format file (https://github.com/lh3/miniasm/blob/master/PAF.md).
74 Thi 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.
74 75
75 Map set of query seqences against a list of reference genomes: 76 Map set of query seqences against a list of reference genomes:
76 77
77 mashmap --rl referenceList.txt -q query.fa 78 mashmap --rl referenceList.txt -q query.fa
78 79