Mercurial > repos > nilesh > sickle
annotate sickle.xml @ 7:dd2ec1b37e96 draft default tip
Added support for the "pair"-data collection type.
| author | yhoogstrate |
|---|---|
| date | Fri, 20 Feb 2015 02:49:16 -0500 |
| parents | c6a790326d1e |
| children |
| rev | line source |
|---|---|
| 3 | 1 <tool id="sickle" name="Sickle"> |
| 2 <description>Windowed Adaptive Trimming of FastQ data</description> | |
| 3 <requirements> | |
| 4 | 4 <requirement type="package" version="1.21.0">sickle</requirement> |
| 3 | 5 </requirements> |
| 6 <command> | |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
7 sickle |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
8 |
| 3 | 9 #if str($readtype.single_or_paired) == "se": |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
10 se |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
11 #else |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
12 pe |
| 3 | 13 #end if |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
14 |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
15 --quiet |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
16 |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
17 #if str($readtype.single_or_paired) == "se": |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
18 -f $input_single -t $qual_type -o $output_single |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
19 #end if |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
20 |
| 3 | 21 #if str($readtype.single_or_paired) == "pe": |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
22 -f $input_paired1 -r $input_paired2 -o $output_paired1 -p $output_paired2 -s $output_paired_single -t $qual_type |
| 3 | 23 #end if |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
24 |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
25 #if str($readtype.single_or_paired) == "pe_collection": |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
26 -f "$readtype.input_paired.forward" -r "$readtype.input_paired.reverse" -o $output_paired1 -p $output_paired2 -s $output_paired_single -t $qual_type |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
27 #end if |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
28 |
| 3 | 29 #if str($qual_threshold) != "": |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
30 -q $qual_threshold |
| 3 | 31 #end if |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
32 |
| 3 | 33 #if str($length_threshold) != "": |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
34 -l $length_threshold |
| 3 | 35 #end if |
| 36 | |
| 37 #if $no_five_prime: | |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
38 -x |
| 3 | 39 #end if |
| 40 | |
| 41 #if $discard_n: | |
| 42 -n | |
| 43 #end if | |
| 44 </command> | |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
45 |
| 3 | 46 <inputs> |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
47 <!-- original wrapper code --> |
| 3 | 48 <conditional name="readtype"> |
| 49 <param name="single_or_paired" type="select" optional="false" label="Single-End or Paired-End reads?"> | |
| 50 <option value="se" selected="true">Single-End</option> | |
| 51 <option value="pe">Paired-End</option> | |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
52 <option value="pe_collection">Paired-End (as collection)</option> |
| 3 | 53 </param> |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
54 |
| 3 | 55 <when value="se"> |
| 56 <param format="fastq, fastqsanger" name="input_single" type="data" optional="false" label="Single-End FastQ Reads"/> | |
| 57 </when> | |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
58 |
| 3 | 59 <when value="pe"> |
| 60 <param format="fastq, fastqsanger" name="input_paired1" type="data" optional="false" label="Paired-End Forward Strand FastQ Reads"/> | |
| 61 <param format="fastq, fastqsanger" name="input_paired2" type="data" optional="false" label="Paired-End Reverse Strand FastQ Reads"/> | |
| 62 </when> | |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
63 |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
64 <when value="pe_collection"> |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
65 <param format="fastq, fastqsanger" name="input_paired" type="data_collection" collection_type="paired" optional="false" label="Paired-End FastQ Reads as paired collection" /> |
|
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
66 </when> |
| 3 | 67 </conditional> |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
68 |
| 3 | 69 <param name="qual_type" type="select" optional="false" label="Quality type"> |
| 70 <option value="illumina" selected="true">Illumina</option> | |
| 71 <option value="solexa">Solexa</option> | |
| 72 <option value="sanger">Sanger</option> | |
| 73 </param> | |
| 74 | |
| 75 <param name="qual_threshold" value="20" type="integer" optional="true" label="Quality Threshold"> | |
| 76 <validator type="in_range" min="0" message="Minimum value is 0"/> | |
| 77 </param> | |
| 78 | |
| 79 <param name="length_threshold" value="20" type="integer" optional="true" label="Length Threshold"> | |
| 80 <validator type="in_range" min="0" message="Minimum value is 0"/> | |
| 81 </param> | |
| 82 | |
| 83 <param name="no_five_prime" type="boolean" label="Don't do 5' trimming"/> | |
| 84 <param name="discard_n" type="boolean" label="Discard sequences with Ns"/> | |
| 85 </inputs> | |
| 86 | |
| 87 <outputs> | |
| 88 <data format_source="input_single" name="output_single" label="Single-End output of ${tool.name} on ${on_string}"> | |
| 89 <filter>(readtype['single_or_paired'] == 'se')</filter> | |
| 90 </data> | |
| 91 | |
| 92 <data format_source="input_paired1" name="output_paired1" label="Paired-End forward strand output of ${tool.name} on ${on_string}"> | |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
93 <filter>(readtype['single_or_paired'] == 'pe' or readtype['single_or_paired'] == 'pe_collection')</filter> |
| 3 | 94 </data> |
| 95 | |
| 96 <data format_source="input_paired2" name="output_paired2" label="Paired-End reverse strand output of ${tool.name} on ${on_string}"> | |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
97 <filter>(readtype['single_or_paired'] == 'pe' or readtype['single_or_paired'] == 'pe_collection')</filter> |
| 3 | 98 </data> |
| 99 | |
| 100 <data format_source="input_paired1" name="output_paired_single" label="Singletons from Paired-End output of ${tool.name} on ${on_string}"> | |
|
7
dd2ec1b37e96
Added support for the "pair"-data collection type.
yhoogstrate
parents:
6
diff
changeset
|
101 <filter>(readtype['single_or_paired'] == 'pe' or readtype['single_or_paired'] == 'pe_collection')</filter> |
| 3 | 102 </data> |
| 103 </outputs> | |
| 6 | 104 <tests> |
| 105 <param name="single_or_paired" value="pe"/> | |
| 106 <param name="input_paired1" value="test.f.fastq"/> | |
| 107 <param name="input_paired2" value="test.r.fastq"/> | |
| 108 <param name="qual_type" value="illumina"/> | |
| 109 <param name="qual_threshold" value="20"/> | |
| 110 <param name="length_threshold" value="20"/> | |
| 111 <param name="no_five_prime" value="false"/> | |
| 112 <param name="discard_n" value="false"/> | |
| 113 <output name="output_paired1" value="output.f.fastq"/> | |
| 114 <output name="output_paired2" value="output.r.fastq"/> | |
| 115 <output name="output_paired_single" value="output.fastq"/> | |
| 116 </tests> | |
| 117 | |
| 3 | 118 <help> |
| 119 Most modern sequencing technologies produce reads that have deteriorating quality towards the 3'-end and some towards the 5'-end as well. Incorrectly called bases in both regions negatively impact assembles, mapping, and downstream bioinformatics analyses. | |
| 120 | |
| 121 Sickle is a tool that uses sliding windows along with quality and length thresholds to determine when quality is sufficiently low to trim the 3'-end of reads and also determines when the quality is sufficiently high enough to trim the 5'-end of reads. It will also discard reads based upon the length threshold. It takes the quality values and slides a window across them whose length is 0.1 times the length of the read. If this length is less than 1, then the window is set to be equal to the length of the read. Otherwise, the window slides along the quality values until the average quality in the window rises above the threshold, at which point the algorithm determines where within the window the rise occurs and cuts the read and quality there for the 5'-end cut. Then when the average quality in the window drops below the threshold, the algorithm determines where in the window the drop occurs and cuts both the read and quality strings there for the 3'-end cut. However, if the length of the remaining sequence is less than the minimum length threshold, then the read is discarded entirely. 5'-end trimming can be disabled. | |
| 122 | |
| 123 Sickle also has an option to discard reads with any Ns in them. | |
| 124 | |
| 125 Sickle supports three types of quality values: Illumina, Solexa, and Sanger. Note that the Solexa quality setting is an approximation (the actual conversion is a non-linear transformation). The end approximation is close. Illumina quality refers to qualities encoded with the CASAVA pipeline between versions 1.3 and 1.7. Illumina quality using CASAVA >= 1.8 is Sanger encoded. | |
| 126 | |
| 127 Note that Sickle will remove the 2nd fastq record header (on the "+" line) and replace it with simply a "+". This is the default format for CASAVA >= 1.8. | |
| 128 | |
| 129 Sickle also supports gzipped file inputs. | |
| 130 </help> | |
| 131 | |
| 132 </tool> |