Mercurial > repos > peterjc > clc_assembly_cell
comparison tools/clc_assembly_cell/clc_mapper.xml @ 5:31a9169b6222 draft default tip
"Update all the pico_galaxy tools on main Tool Shed"
author | peterjc |
---|---|
date | Fri, 16 Apr 2021 22:32:53 +0000 |
parents | 46a667b5e19e |
children |
comparison
equal
deleted
inserted
replaced
4:46a667b5e19e | 5:31a9169b6222 |
---|---|
1 <tool id="clc_mapper" name="CLC Mapper" version="0.0.6"> | 1 <tool id="clc_mapper" name="CLC Mapper" version="0.0.7"> |
2 <description>Maps reads giving a SAM/BAM file</description> | 2 <description>Maps reads giving a SAM/BAM file</description> |
3 <requirements> | 3 <requirements> |
4 <requirement type="binary">clc_mapper</requirement> | |
5 <requirement type="binary">clc_cas_to_sam</requirement> | |
6 <requirement type="binary">samtools</requirement> | |
7 <requirement type="package" version="0.1.19">samtools</requirement> | 4 <requirement type="package" version="0.1.19">samtools</requirement> |
8 </requirements> | 5 </requirements> |
9 <stdio> | 6 <version_command> |
10 <!-- Assume anything other than zero is an error --> | 7 ${CLC_ASSEMBLY_CELL:-/mnt/apps/clcBio/clc-assembly-cell-4.1.0-linux_64/}clc_mapper | grep -i version |
11 <exit_code range="1:" /> | 8 </version_command> |
12 <exit_code range=":-1" /> | 9 <command detect_errors="aggressive"> |
13 </stdio> | 10 echo Mapping reads with clc_mapper... |
14 <version_command>${CLC_ASSEMBLY_CELL:-/mnt/apps/clcBio/clc-assembly-cell-4.1.0-linux_64/}clc_mapper | grep -i version</version_command> | |
15 <command>echo Mapping reads with clc_mapper... | |
16 && \${CLC_ASSEMBLY_CELL:-/mnt/apps/clcBio/clc-assembly-cell-4.1.0-linux_64/}clc_mapper | 11 && \${CLC_ASSEMBLY_CELL:-/mnt/apps/clcBio/clc-assembly-cell-4.1.0-linux_64/}clc_mapper |
17 #for $ref in $references | 12 #for $ref in $references |
18 #if str($ref.ref_type)=="circular" | 13 #if str($ref.ref_type)=="circular" |
19 -d -z "$ref.ref_file" | 14 -d -z '$ref.ref_file' |
20 #else | 15 #else |
21 -d "$ref.ref_file" | 16 -d '$ref.ref_file' |
22 #end if | 17 #end if |
23 #end for | 18 #end for |
24 #for $rg in $read_group | 19 #for $rg in $read_group |
25 ##-------------------------------------- | 20 ##-------------------------------------- |
26 #if str($rg.segments.type) == "paired" | 21 #if str($rg.segments.type) == "paired" |
27 -p $rg.segments.placement $rg.segments.dist_mode $rg.segments.min_size $rg.segments.max_size -q -i "$rg.segments.filename1" "$rg.segments.filename2" | 22 -p $rg.segments.placement $rg.segments.dist_mode $rg.segments.min_size $rg.segments.max_size -q -i '$rg.segments.filename1' '$rg.segments.filename2' |
28 #end if | 23 #end if |
29 ##-------------------------------------- | 24 ##-------------------------------------- |
30 #if str($rg.segments.type) == "interleaved" | 25 #if str($rg.segments.type) == "interleaved" |
31 -p $rg.segments.placement $rg.segments.dist_mode $rg.segments.min_size $rg.segments.max_size -q "$rg.segments.filename" | 26 -p $rg.segments.placement $rg.segments.dist_mode $rg.segments.min_size $rg.segments.max_size -q '$rg.segments.filename' |
32 #end if | 27 #end if |
33 ##-------------------------------------- | 28 ##-------------------------------------- |
34 #if str($rg.segments.type) == "none" | 29 #if str($rg.segments.type) == "none" |
35 -p no -q | 30 -p no -q |
36 #for $f in $rg.segments.filenames | 31 #for $f in $rg.segments.filenames |
37 "$f" | 32 '$f' |
38 #end for | 33 #end for |
39 #end if | 34 #end if |
40 ##-------------------------------------- | 35 ##-------------------------------------- |
41 #end for | 36 #end for |
42 -o "temp_job.cas" | 37 -o "temp_job.cas" |
43 --cpus \${GALAXY_SLOTS:-4} | 38 --cpus "\${GALAXY_SLOTS:-4}" |
44 ## TODO - filtering out the progress lines seems to mess up the multiple commands | 39 ## TODO - filtering out the progress lines seems to mess up the multiple commands |
45 ## | grep -v "^Progress: " | 40 ## | grep -v "^Progress: " |
46 ##=========================================== | 41 ##=========================================== |
47 ## TODO - I've required all the input in Sanger FASTQ format (or FASTA) so can | 42 ## TODO - I've required all the input in Sanger FASTQ format (or FASTA) so can |
48 ## use the offset 33, rather then the CLCbio default of 64 which is only for | 43 ## use the offset 33, rather then the CLCbio default of 64 which is only for |
49 ## obsolete Illumina FASTQ files. Really need this option per input file... | 44 ## obsolete Illumina FASTQ files. Really need this option per input file... |
50 && echo Converting CAS file to BAM with clc_cas_to_sam... | 45 && echo Converting CAS file to BAM with clc_cas_to_sam... |
51 && \${CLC_ASSEMBLY_CELL:-/mnt/apps/clcBio/clc-assembly-cell-4.1.0-linux_64/}clc_cas_to_sam --cas "temp_job.cas" -o "temp_job.bam" --no-progress --qualityoffset 33 | 46 && \${CLC_ASSEMBLY_CELL:-/mnt/apps/clcBio/clc-assembly-cell-4.1.0-linux_64/}clc_cas_to_sam --cas 'temp_job.cas' -o 'temp_job.bam' --no-progress --qualityoffset 33 |
52 #if $discard_unmapped: | 47 #if $discard_unmapped: |
53 ## -u / --discardunmapped: Discard the unmapped reads | 48 ## -u / --discardunmapped: Discard the unmapped reads |
54 -u | 49 -u |
55 #end if | 50 #end if |
56 && rm "temp_job.cas" | 51 && rm 'temp_job.cas' |
57 ##=========================================== | 52 ##=========================================== |
58 && echo Sorting BAM file with samtools... | 53 && echo Sorting BAM file with samtools... |
59 && samtools sort "temp_job.bam" "temp_sorted" | 54 && samtools sort "temp_job.bam" "temp_sorted" |
60 && mv "temp_sorted.bam" "$out_bam" | 55 && mv "temp_sorted.bam" "$out_bam" |
61 && echo Indexing BAM file with samtools... | 56 && echo Indexing BAM file with samtools... |
145 <param name="read_group_0|segments|min_size" value="1" /> | 140 <param name="read_group_0|segments|min_size" value="1" /> |
146 <param name="read_group_0|segments|max_size" value="1000" /> | 141 <param name="read_group_0|segments|max_size" value="1000" /> |
147 <param name="read_group_0|segments|dist_mode" value="ss" /> | 142 <param name="read_group_0|segments|dist_mode" value="ss" /> |
148 <param name="read_group_0|segments|filename" value="SRR639755_mito_pairs.fastq.gz" ftype="fastqsanger" /> | 143 <param name="read_group_0|segments|filename" value="SRR639755_mito_pairs.fastq.gz" ftype="fastqsanger" /> |
149 <param name="discard_unmapped" value="false" /> | 144 <param name="discard_unmapped" value="false" /> |
150 <output name="out_fasta" file="SRR639755_mito_pairs_vs_NC_010642_clc.bam" ftype="bam" lines_diff="4"/> | 145 <output name="out_bam" file="SRR639755_mito_pairs_vs_NC_010642_clc.bam" ftype="bam" lines_diff="4"/> |
151 </test> | 146 </test> |
152 <test> | 147 <test> |
153 <param name="ref_file" value="NC_010642.fna" ftype="fasta" /> | 148 <param name="ref_file" value="NC_010642.fna" ftype="fasta" /> |
154 <param name="ref_type" value="circular" /> | 149 <param name="ref_type" value="circular" /> |
155 <param name="read_group_0|segments|type" value="interleaved" /> | 150 <param name="read_group_0|segments|type" value="interleaved" /> |
158 <param name="read_group_0|segments|min_size" value="1" /> | 153 <param name="read_group_0|segments|min_size" value="1" /> |
159 <param name="read_group_0|segments|max_size" value="1000" /> | 154 <param name="read_group_0|segments|max_size" value="1000" /> |
160 <param name="read_group_0|segments|dist_mode" value="ss" /> | 155 <param name="read_group_0|segments|dist_mode" value="ss" /> |
161 <param name="read_group_0|segments|filename" value="SRR639755_mito_pairs.fastq.gz" ftype="fastqsanger" /> | 156 <param name="read_group_0|segments|filename" value="SRR639755_mito_pairs.fastq.gz" ftype="fastqsanger" /> |
162 <param name="discard_unmapped" value="true" /> | 157 <param name="discard_unmapped" value="true" /> |
163 <output name="out_fasta" file="SRR639755_mito_pairs_vs_NC_010642_clc_u.bam" ftype="bam" lines_diff="4"/> | 158 <output name="out_bam" file="SRR639755_mito_pairs_vs_NC_010642_clc_u.bam" ftype="bam" lines_diff="4"/> |
164 </test> | 159 </test> |
165 </tests> | 160 </tests> |
166 <help> | 161 <help> |
167 | 162 |
168 **What it does** | 163 **What it does** |
188 Tool Shed at http://toolshed.g2.bx.psu.edu/view/peterjc/clc_assembly_cell | 183 Tool Shed at http://toolshed.g2.bx.psu.edu/view/peterjc/clc_assembly_cell |
189 </help> | 184 </help> |
190 <citations> | 185 <citations> |
191 <citation type="bibtex"> | 186 <citation type="bibtex"> |
192 @MISC{clcbio, | 187 @MISC{clcbio, |
193 AUTHOR = {CLC Bio} | 188 AUTHOR = {CLC Bio} |
194 title = {{CLC Assembly Cell}}, | 189 title = {{CLC Assembly Cell}}, |
195 url = {https://www.qiagenbioinformatics.com/products/clc-assembly-cell/}, | 190 url = {https://www.qiagenbioinformatics.com/products/clc-assembly-cell/}, |
196 year = {2008--2015} | 191 year = {2008--2015} |
197 } | 192 } |
198 </citation> | 193 </citation> |