Mercurial > repos > fubar > mashmap
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 |