Mercurial > repos > iuc > umi_tools_count

annotate umi-tools_counts.xml @ 1:3c932ad4a174 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
author iuc
date Sat, 14 Jul 2018 06:14:24 -0400
parents 8db56d2f8b72
children b557acca0b56
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
1 <tool id="umi_tools_count" name="UMI-tools count" version="@VERSION@.1">
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
2 <description>performs quantification of UMIs from BAM files</description>
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
3 <macros>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
4 <import>macros.xml</import>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
5 <xml name="sanitize_tag" >
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
6 <sanitizer invalid_char="">
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
7 <valid initial="string.letters,string.digits" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
8 </sanitizer>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
9 </xml>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
10 </macros>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
11 <expand macro="requirements" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
12 <command detect_errors="exit_code"><![CDATA[
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
13 ln -s '${input_bam}' 'input.bam' &&
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
14 ln -s '${input_bam.metadata.bam_index}' 'input.bam.bai' &&
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
15
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
16 umi_tools count
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
17 -I input.bam
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
18 '$paired'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
19 --extract-umi-method='$barcodes.extract_umi_method.value'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
20 #if str($barcodes.extract_umi_method) == 'read_id':
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
21 --umi-separator='$barcodes.umi_separator.value'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
22 #else if str($barcodes.extract_umi_method) == 'tag':
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
23 --umi-tag='$barcodes.umi_tag.value'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
24 --cell-tag='$barcodes.cell_tag.value'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
25 #end if
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
26 --method='$method.value'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
27 --edit-distance-threshold='$edit_distance_threshold'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
28 --mapping-quality='$advanced.mapping_quality'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
29 --per-gene
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
30 '$wide_format_cell_counts'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
31 '$advanced.per_contig'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
32 '$advanced.per_cell'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
33 #if str($advanced.gene_tag) != "":
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
34 --gene-tag='$advanced.gene_tag.value'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
35 #end if
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
36 #if str($advanced.skip_tags_regex) != "":
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
37 --skip-tags-regex='$advanced.skip_tags_regex.value'
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
38 #end if
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
39 #if '$advanced.random_seed' != 0:
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
40 --random-seed='$advanced.random_seed'
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
41 #end if
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
42 -S '$out_counts'
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
43 ]]></command>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
44 <inputs>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
45 <param name="input_bam" type="data" format="bam" label="Sorted BAM file" help="Please use the samtools sort tool to ensure a correct BAM input" />
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
46 <param argument="--paired" type="boolean" truevalue="--paired" falsevalue="" checked="false" label="Bam is paired-end" help="both read pairs will be output. This will also force the use of the template length to determine reads with the same mapping coordinates." />
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
47 <conditional name="barcodes" >
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
48 <param argument="--extract-umi-method" name="extract_umi_method" type="select" label="Umi Extract Method" help="How are the barcodes encoded in the read?" >
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
49 <option value="read_id" selected="true">Barcodes are contained at the end of the read seperated by a delimiter</option>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
50 <option value="tag" >Barcodes are contained in tags</option>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
51 <option value="umis" >Barcodes were extracted using umis</option>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
52 </param>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
53 <when value="read_id" >
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
54 <param argument="--umi-separator" name="umi_separator" type="text" label="Delimiter between read id and the UMI" value="_" >
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
55 <sanitizer invalid_char="" >
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
56 <valid initial="string.punctuation" />
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
57 </sanitizer>
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
58 </param>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
59 </when>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
60 <when value="tag" >
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
61 <param argument="--umi-tag" name="umi_tag" type="text" label="Tag which contains the UMI" >
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
62 <expand macro="sanitize_tag" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
63 </param>
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
64 <param argument="--cell-tag" name="cell_tag" type="text" label="Tag which contains the cell barcode" >
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
65 <expand macro="sanitize_tag" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
66 </param>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
67 </when>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
68 <when value="umis"></when>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
69 </conditional>
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
70 <param argument="--method" type="select" label="Method to identify group of reads" help="UMIs with the same (or similar) codes can be grouped together. The simplest methods 'unique' and 'percentile' group identical
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
71 UMIs, however 'cluster', 'adjacency', and 'directional' can group similar umis with edit distances less than some threshold. Unique: Reads group share the exact same UMI. Percentile: Reads group share the same UMI, and UMIs with
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
72 counts &lt; 1% of the median counts for UMIs at the same position are ignored. Cluster: Identify clusters of connected UMIs (based on hamming distance threshold). Adjacency: Same as cluster, but considers only directly ajacent UMIs in the cluster. Directional: Identify cluster of connected UMIs based on hamming distance and umi." >
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
73 <option value="unique" >Unique</option>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
74 <option value="percentile">Percentile</option>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
75 <option value="cluster">Cluster</option>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
76 <option value="adjacency">Adjacency</option>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
77 <option value="directional" selected="true" >Directional</option>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
78 </param>
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
79 <param argument="--edit-distance-threshold" name="edit_distance_threshold" type="integer" label="Edit distance threshold" min="0" value="1" />
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
80 <param argument="--wide-format-cell-counts" name="wide_format_cell_counts" type="boolean" truevalue="--wide-format-cell-counts" falsevalue="" checked="true" label="Output a matrix of genes and cells, instead of a flat file" />
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
81 <section name="advanced" title="Extra parameters" >
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
82 <param argument="--mapping-quality" name="mapping_quality" type="integer" min="0" value="0" label="Minimum mapping quality" />
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
83 <!-- Currently hard-coded parameter. Leave here if useful to future wrapper -->
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
84 <!-- <param argument="-\-per-gene" name="per_gene" type="text" label="Group reads together if they have the same gene" help="Reads will be grouped together if they have the same gene. This is useful if your library
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
85 prep generates PCR duplicates with non-identical alignment positions such as CEL-Seq. Note this option is hardcoded to be on with the count command. I.e counting is always performed per-gene. Must be combined with either
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
86 -\-gene-tag or -\-per-contig option" /> -->
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
87 <param argument="--gene-tag" name="gene_tag" type="text" label="Deduplicate per gene." help="The gene information is encoded in the bam read tag." value="XT" >
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
88 <expand macro="sanitize_tag" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
89 </param>
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
90 <param argument="--skip-tags-regex" name="skip_tags_regex" type="text" label="Skip any reads where the gene matches this tag" value="" >
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
91 <sanitizer invalid_char="">
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
92 <valid initial="string.letters,string.digits">
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
93 <add value="!="/>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
94 <add value="-"/>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
95 <add value="_"/>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
96 <add value="."/>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
97 <add value="?"/>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
98 <add value="&lt;"/><!-- left triangle bracket -->
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
99 <add value="&gt;"/><!-- right triangle bracket -->
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
100 <add value="&#91;"/> <!-- left square bracket -->
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
101 <add value="&#93;"/> <!-- right square bracket -->
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
102 <add value="&#94;"/> <!-- caret -->
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
103 <add value="&#123;"/> <!-- left curly -->
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
104 <add value="&#125;"/> <!-- right curly -->
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
105 <add value="&#40;"/> <!-- left parenthesis -->
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
106 <add value="&#41;"/> <!-- right parenthesis -->
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
107 </valid>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
108 </sanitizer>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
109 </param>
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
110 <param argument="--per-contig" name="per_contig" type="boolean" truevalue="--per-contig" falsevalue="" checked="false" label="Deduplicate per contig (field 3 in BAM; RNAME)" help="All reads with the same contig will be considered to have the same alignment position. This is useful if you have aligned to a reference transcriptome with one transcript per gene." />
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
111 <param argument="--per-cell" name="per_cell" type="boolean" truevalue="--per-cell" falsevalue="" checked="true" label="Group reads only if they have the same cell barcode." />
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
112 <param argument="--random-seed" name="random_seed" type="integer" min="0" value="0" label="Random Seed" />
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
113 </section>
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
114 </inputs>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
115 <outputs>
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
116 <data name="out_counts" format="tabular" />
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
117 </outputs>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
118 <tests>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
119 <test><!--count_single_gene_tag:-->
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
120 <param name="input_bam" value="chr19_gene_tags.bam" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
121 <param name="random_seed" value="123456789" />
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
122 <param name="method" value="directional" />
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
123 <param name="gene_tag" value="XF" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
124 <param name="skip_tags_regex" value="^[__|Unassigned]" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
125 <param name="extract_umi_method" value="umis" />
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
126 <param name="wide_format_cell_counts" value="false" />
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
127 <param name="per_cell" value="false" />
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
128 <output name="out_counts" value="count_single_gene_tag.tsv" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
129 </test>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
130 <test><!--count_single_cells_gene_tag:-->
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
131 <param name="input_bam" value="chr19_gene_tags.bam" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
132 <param name="random_seed" value="123456789" />
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
133 <param name="method" value="directional" />
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
134 <param name="gene_tag" value="XF" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
135 <param name="skip_tags_regex" value="^[__|Unassigned]" />
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
136 <param name="per_cell" value="true" />
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
137 <param name="extract_umi_method" value="umis" />
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
138 <param name="wide_format_cell_counts" value="false" />
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
139 <output name="out_counts" value="count_single_cells_gene_tag.tsv" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
140 </test>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
141 <test><!--count_single_cells_wide_gene_tag:-->
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
142 <param name="input_bam" value="chr19_gene_tags.bam" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
143 <param name="random_seed" value="123456789" />
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
144 <param name="method" value="directional" />
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
145 <param name="gene_tag" value="XF" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
146 <param name="skip_tags_regex" value="^[__|Unassigned]" />
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
147 <param name="per_cell" value="true" />
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
148 <param name="extract_umi_method" value="umis" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
149 <param name="wide_format_cell_counts" value="true" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
150 <output name="out_counts" value="count_single_cells_gene_tag_wide.tsv" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
151 </test>
1
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
152 <test><!-- count ENSDARG00000019692, with defaults -->
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
153 <param name="input_bam" value="fc.ENSDARG00000019692.bam" />
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
154 <param name="method" value="unique" />
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
155 <output name="out_counts" value="fc.ENSDARG00000019692.counts" />
3c932ad4a174 planemo upload commit 9a3aeb2c588f9f67824ea5568923ce70b048499a
iuc
parents: 0
diff changeset
156 </test>
0
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
157 </tests>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
158 <help><![CDATA[
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
159
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
160 UMI Tools count - Count reads per gene from BAM using UMIs
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
161 ----------------------------------------------------------
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
162
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
163 Purpose
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
164 -------
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
165
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
166 The purpose of this command is to count the number of reads per gene based
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
167 on the mapping co-ordinate and the UMI attached to the read.
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
168
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
169
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
170 It is assumed that the FASTQ files were processed with extract_umi.py
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
171 before mapping and thus the UMI is the last word of the read name. e.g:
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
172
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
173 @HISEQ:87:00000000_AATT
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
174
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
175 where AATT is the UMI sequeuence.
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
176
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
177 If you have used an alternative method which does not separate the
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
178 read id and UMI with a "_", such as bcl2fastq which uses ":", you can
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
179 specify the separator, or if your UMIs are encoded in a tag you can also specify this.
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
180
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
181 ]]></help>
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
182 <expand macro="citations" />
8db56d2f8b72 planemo upload commit c79a5f4a05156bb2a6035a844aa9ad8f0e59ecb5
iuc
parents:
diff changeset
183 </tool>