annotate smudgeplot.xml @ 0:4c91967b3e6d draft default tip

planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
author iuc
date Thu, 30 Jun 2022 22:44:50 +0000
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
1 <tool id="smudgeplot" name="Smudgeplot" version="@TOOL_VERSION@+galaxy+@VERSION_SUFFIX@" profile="21.05">
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
2 <description>inference of ploidy and heterozygosity structure using whole genome sequencing</description>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
3
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
4 <macros>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
5 <token name="@TOOL_VERSION@">0.2.5</token>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
6 <token name="@VERSION_SUFFIX@">1</token>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
7 </macros>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
8
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
9 <xrefs>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
10 <xref type="bio.tools">smudgeplots</xref>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
11 </xrefs>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
12
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
13 <requirements>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
14 <requirement type="package" version="@TOOL_VERSION@">smudgeplot</requirement>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
15 <requirement type="package" version="2.3.0">kmer-jellyfish</requirement>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
16 </requirements>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
17
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
18 <command detect_errors="exit_code"><![CDATA[
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
19
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
20 #if $file.input.input_select == 'reads'
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
21
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
22 ## ~~~~~~~~~~~~~~~ Generate kmer-dump with presets ~~~~~~~~~~~~~~~~~~~~~
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
23
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
24 ## Jellyfish kmer count
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
25 ## ---------------------------------------------------------------------
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
26
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
27 #if $file.input.reads[0].is_of_type("fastqsanger.gz") or $file.input.reads[0].is_of_type("fasta.gz")
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
28 gunzip -c
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
29 #for $f in $file.input.reads
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
30 #if $f
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
31 '$f'
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
32 #end if
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
33 | jellyfish count -m 21 -t 4 -s 1M -o 1_counts.jf -C /dev/stdin
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
34 #end for
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
35
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
36 #else
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
37 jellyfish count -m 21 -t 4 -s 1M -o 1_counts.jf -C
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
38 #for $f in $file.input.reads
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
39 #if $f
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
40 '$f'
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
41 #end if
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
42 #end for
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
43 #end if
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
44
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
45 && jellyfish histo 1_counts.jf > 1_kmer_k21.hist
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
46
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
47 ## Calculate lower and upper kmer count cutoffs
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
48 ## ---------------------------------------------------------------------
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
49
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
50 #if $file.input.lower_cutoff is not None:
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
51 && L=$file.input.lower_cutoff
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
52 #else
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
53 && L=\$(smudgeplot.py cutoff 1_kmer_k21.hist L)
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
54 #end if
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
55
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
56 #if $file.input.upper_cutoff is not None:
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
57 && U=$file.input.upper_cutoff
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
58 #else
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
59 && U=\$(smudgeplot.py cutoff 1_kmer_k21.hist U)
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
60 #end if
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
61
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
62 ## ---------------------------------------------------------------------
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
63 ## Dump and extract coverage
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
64
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
65 && echo "Dump with cutoffs L=\$L, U=\$U"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
66 && jellyfish dump -c -L \$L -U \$U 1_counts.jf > 2_dump.jf
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
67 && smudgeplot.py hetkmers -o 2_kmer_pairs 2_dump.jf
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
68
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
69 #else
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
70
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
71 ## ~~~~~~~~~~~~~~~~~~~ Use provided kmer dump ~~~~~~~~~~~~~~~~~~~~~~~~~~
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
72
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
73 smudgeplot.py hetkmers -o 2_kmer_pairs '$file.input.dump'
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
74
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
75 #end if
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
76
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
77 ## ---------------------------------------------------------------------
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
78 ## Plot
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
79
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
80 && smudgeplot.py plot 2_kmer_pairs_coverages.tsv -o my_genome
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
81
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
82 ]]></command>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
83
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
84 <inputs>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
85 <section name="file" title="File inputs" expanded="true">
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
86 <conditional name="input">
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
87 <param
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
88 name="input_select" type="select" label="Select input type"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
89 help="For more control, create your own Kmer dump using Jellyfish.
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
90 See Smudgeplot on GitHub for more details: https://github.com/KamilSJaron/smudgeplot"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
91 >
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
92 <option value="reads" selected="true">Sequencing reads</option>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
93 <option value="dump">Kmer dump file</option>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
94 </param>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
95
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
96 <when value="reads">
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
97 <param
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
98 name="reads" type="data" format="fastqsanger,fastqsanger.gz,fasta.gz,fasta"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
99 label="Sequencing reads" multiple="true"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
100 help="Sequencing reads corresponding to your genome.
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
101 Don't worry about read pairing as it is not used in Kmer-counting.
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
102 If selecting multiple datasets, please do not mix datatypes!"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
103 />
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
104
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
105 <param
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
106 name="lower_cutoff"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
107 label="Lower kmer cutoff"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
108 type="integer"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
109 optional="true"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
110 help="Optionally set a manual lower limit for filtering kmers with
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
111 smudgeplot hetkmers. If no value is set, a cutoff will be
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
112 estimated with smudgeplot cutoff. Use the GenomeScope tool to
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
113 visualize your kmer histogram when choosing cutoff values."
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
114 />
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
115
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
116 <param
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
117 name="upper_cutoff"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
118 label="Upper kmer cutoff"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
119 type="integer"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
120 optional="true"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
121 help="Optionally set a manual upper limit for filtering kmers with
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
122 smudgeplot hetkmers. If no value is set, a cutoff will be
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
123 estimated with smudgeplot cutoff. Use the GenomeScope tool to
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
124 visualize your kmer histogram when choosing cutoff values."
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
125 />
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
126 </when>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
127
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
128 <when value="dump">
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
129 <param
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
130 name="dump" type="data" format="txt"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
131 label="Kmer dump"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
132 help="Upload your own Kmer dump file created with the Jellyfish or KMC tool.
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
133 This enables control over kmer-counting parameters."
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
134 />
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
135 </when>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
136 </conditional>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
137 </section>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
138
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
139 <param name="table_output" type="boolean" label="Output summary table"></param>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
140 <param name="verbose_output" type="boolean" label="Output verbose summary"></param>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
141 <param name="warnings_output" type="boolean" label="Output genome warnings"></param>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
142 </inputs>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
143
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
144 <outputs>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
145 <data
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
146 name="smudgeplot" format="png"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
147 from_work_dir="my_genome_smudgeplot.png"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
148 label="${tool.name} on ${on_string}: Smudgeplot"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
149 />
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
150 <data
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
151 name="smudgeplot_log" format="png"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
152 from_work_dir="my_genome_smudgeplot_log10.png"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
153 label="${tool.name} on ${on_string}: Smudgeplot (log10)"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
154 />
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
155 <data
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
156 name="genome_summary" format="tabular"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
157 from_work_dir="my_genome_summary_table.tsv"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
158 label="${tool.name} on ${on_string}: Genome summary table"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
159 >
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
160 <filter>table_output</filter>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
161 </data>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
162 <data
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
163 name="genome_summary_verbose" format="txt"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
164 from_work_dir="my_genome_verbose_summary.txt"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
165 label="${tool.name} on ${on_string}: Genome verbose summary"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
166 >
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
167 <filter>verbose_output</filter>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
168 </data>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
169 <data
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
170 name="genome_warnings" format="txt"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
171 from_work_dir="my_genome_warnings.txt"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
172 label="${tool.name} on ${on_string}: Genome warnings"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
173 >
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
174 <filter>warnings_output</filter>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
175 </data>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
176 </outputs>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
177
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
178 <tests>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
179 <!-- Standard run -->
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
180 <test expect_num_outputs="2">
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
181 <param name="input_select" value="reads"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
182 <param name="reads" value="test_reads.fasta" ftype="fasta"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
183 <param name="lower_cutoff" value="2"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
184 <param name="upper_cutoff" value="25"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
185 <output name="smudgeplot" ftype="png" file="my_genome_smudgeplot.png" compare="sim_size"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
186 </test>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
187 <!-- Standard run with gzipped input -->
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
188 <test expect_num_outputs="2">
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
189 <param name="input_select" value="reads"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
190 <param name="reads" value="test_reads.fasta.gz" ftype="fasta.gz"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
191 <param name="lower_cutoff" value="2"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
192 <param name="upper_cutoff" value="25"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
193 <output name="smudgeplot" ftype="png" file="my_genome_smudgeplot.png" compare="sim_size"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
194 </test>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
195 <!-- Multiple input read files -->
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
196 <test expect_num_outputs="2">
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
197 <param name="input_select" value="reads"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
198 <param name="lower_cutoff" value="2"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
199 <param name="upper_cutoff" value="80"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
200 <param
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
201 name="reads"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
202 value="test_reads.fasta,test_reads_2.fasta,test_reads_3.fasta"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
203 ftype="fasta"
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
204 />
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
205 <output name="smudgeplot" ftype="png" file="my_genome_smudgeplot.png" compare="sim_size"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
206 </test>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
207 <!-- With additional outputs-->
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
208 <test expect_num_outputs="5">
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
209 <param name="input_select" value="reads"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
210 <param name="reads" value="test_reads.fasta" ftype="fasta"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
211 <param name="lower_cutoff" value="2"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
212 <param name="upper_cutoff" value="25"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
213 <param name="table_output" value="true"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
214 <param name="verbose_output" value="true"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
215 <param name="warnings_output" value="true"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
216 <output name="smudgeplot" ftype="png" file="my_genome_smudgeplot.png" compare="sim_size"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
217 <output name="genome_summary" ftype="tabular" file="my_genome_summary_table.tsv"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
218 <output name="genome_summary_verbose" ftype="txt" file="my_genome_verbose_summary.txt"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
219 <output name="genome_warnings" ftype="txt" file="my_genome_warnings.txt"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
220 </test>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
221 <!-- K-mer dump input -->
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
222 <test expect_num_outputs="2">
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
223 <param name="input_select" value="dump"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
224 <param name="dump" value="dump.jf" ftype="txt"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
225 <output name="smudgeplot" ftype="png" file="my_genome_smudgeplot.png" compare="sim_size"/>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
226 </test>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
227 </tests>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
228
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
229 <help><![CDATA[
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
230
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
231 .. class:: infomark
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
232
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
233 **What it does**
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
234
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
235 This tool extracts heterozygous kmer pairs from kmer count databases and performs gymnastics with them. We are able to disentangle genome structure by comparing the sum of kmer pair coverages (CovA + CovB) to their relative coverage (CovB / (CovA + CovB)). Such an approach also allows us to analyze obscure genomes with duplications, various ploidy levels, etc.
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
236
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
237 Smudgeplots are computed from raw or even better from trimmed reads and show the haplotype structure using heterozygous kmer pairs. For example:
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
238
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
239 .. image:: $PATH_TO_IMAGES/smudge.png
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
240 :height: 520
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
241 :alt: Example smudgeplot graph
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
242
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
243 Every haplotype structure has a unique smudge on the graph and the heat of the smudge indicates how frequently the haplotype structure is represented in the genome compared to the other structures. The image above is an ideal case, where the sequencing coverage is sufficient to beautifully separate all the smudges, providing very strong and clear evidence of triploidy.
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
244
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
245 Please see `Smudgeplot on GitHub <https://github.com/KamilSJaron/smudgeplot>`_
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
246 for further documentation and tutorials.
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
247
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
248 **Inputs**
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
249
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
250 You have two choices when running Smudgeplot in Galaxy:
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
251
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
252 1. Input reads file(s) for default kmer-counting with Jellyfish
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
253
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
254 This should be at least one file which providing coverage of your genome of interest.
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
255 The tool accepts compressed (.gz) inputs. If choosing this option, you can
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
256 (optionally) specify manual cutoff values for the kmer dump step. The Smudgeplot
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
257 docs suggest that you can use GenomeScope on a kmer histogram in order to choose
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
258 reasonable lower and upper cutoff values.
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
259
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
260 2. Input your own kmer dump file for more control of kmer counting parameters
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
261
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
262 This file would be created by running ``jellyfish count`` and then ``jellyfish dump`` - the process is well described
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
263 `on GitHub <https://github.com/KamilSJaron/smudgeplot>`_.
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
264
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
265 **Outputs**
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
266
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
267 - ``smudgeplot.png`` smudgeplot image
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
268 - ``smudgeplot_log10.png`` smudgeplot with log scale
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
269 - ``my_genome_summary.tsv`` summarized genome statistics
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
270 - ``my_genome_verbose.txt`` detailed genome statistics
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
271 - ``my_genome_warnings.txt`` warnings emitted from the Smudgeplot tool
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
272
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
273 **Default operation**
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
274
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
275 If choosing reads as the input, a default kmer counting procedure will be used
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
276 to create a kmer dump. This default process is summarized as follows:
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
277
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
278 - ``jellyfish count -m 21 > counts.jf``
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
279 - ``jellyfish histo counts.jf > counts.hist``
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
280 - ``smudgeplot.py cutoff counts.hist`` to get kmer cutoff values (U & L)
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
281 - ``jellyfish dump -c -L <L> -U <U> counts.jf > dump.jf``
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
282
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
283 The kmer dump file is then used to create a smudgeplot:
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
284
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
285 - ``smudgeplot.py hetkmers -o kmer_pairs dump.jf``
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
286 - ``smudgeplot.py plot kmer_pairs_coverages.tsv -o my_genome``
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
287
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
288 ]]></help>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
289 <citations>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
290 <citation type="doi">10.1038/s41467-020-14998-3</citation>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
291 </citations>
4c91967b3e6d planemo upload for repository https://github.com/galaxyproject/tools-iuc commit cbe90253166d9908b68beb36b9488478178d225b
iuc
parents:
diff changeset
292 </tool>