Mercurial > repos > simon-gladman > trimmomatic_032
annotate trimmomatic.xml @ 7:76d1aebc44bf draft
Uploaded: added dataset collection support for paired end trimmomatic
| author | simon-gladman |
|---|---|
| date | Thu, 10 Jul 2014 13:24:53 -0400 |
| parents | 9db0c1c78236 |
| children | 2c4e0eab2cae |
| rev | line source |
|---|---|
| 1 | 1 <tool id="trimmomatic" name="Trimmomatic" version="1.1.0"> |
| 2 | |
| 3 <description>Trim adapters and otherwise cleanup fastq files in a pair aware manner using trimmomatic 0.32</description> | |
| 4 <requirements> | |
| 5 <requirement type="set_environment">JAR_PATH</requirement> | |
| 6 </requirements> | |
| 7 | |
| 8 <command interpreter="perl"> | |
| 9 trimmomatic_wrapper.pl | |
| 10 paired ${paired.is_paired} | |
|
7
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
11 #if $paired.is_paired == "single": |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
12 fwdfile ${paired.forwards_file} |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
13 #end if |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
14 #if $paired.is_paired == "paired": |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
15 fwdfile ${paired.forwards_file} |
| 1 | 16 revfile ${paired.reverse_file} |
| 17 #end if | |
|
7
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
18 #if $paired.is_paired == "collection": |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
19 fwdfile ${collection_data.forward} |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
20 revfile ${collection_data.reverse} |
| 1 | 21 phred $phred |
| 22 cutadapt ${adapt.adapters} | |
| 23 #if $adapt.adapters: | |
| 24 adaptfile ${adapt.adapt_file} | |
| 25 adaptseed ${adapt.adapt_seed} | |
| 26 adaptpalindrome ${adapt.adapt_palindrome} | |
| 27 adaptsimple ${adapt.adapt_simple} | |
| 28 #end if | |
| 29 slidingwindow ${sliding.slidingwindow} | |
| 30 #if $sliding.slidingwindow: | |
| 31 slidingsize ${sliding.slid_window} | |
| 32 slidingqual ${sliding.slid_qual} | |
| 33 #end if | |
| 34 trimleading ${leading.trimleading} | |
| 35 #if $leading.trimleading: | |
| 36 leadingqual ${leading.lead_qual} | |
| 37 #end if | |
| 38 trimtrailing ${trailing.trimtrailing} | |
| 39 #if $trailing.trimtrailing: | |
| 40 trailingqual ${trailing.trail_qual} | |
| 41 #end if | |
| 42 crop ${crop.cropend} | |
| 43 #if $crop.cropend: | |
| 44 croplen ${crop.crop_len} | |
| 45 #end if | |
| 46 headcrop ${headcrop.headcrop_s} | |
| 47 #if $headcrop.headcrop_s: | |
| 48 headcroplen ${headcrop.headcrop_len} | |
| 49 #end if | |
| 50 minlen $minlen | |
| 51 log $log | |
| 52 #if $log == "True": | |
| 53 logfile $logfile | |
| 54 #end if | |
| 55 singles $singles | |
| 56 #if $paired.is_paired: | |
| 57 fwdpairs $fwdpairs | |
| 58 revpairs $revpairs | |
| 59 #end if | |
| 60 tool-dir \$JAR_PATH | |
| 4 | 61 threads \${GALAXY_SLOTS:-4} |
| 1 | 62 > $dummy_out |
| 63 </command> | |
| 64 | |
| 65 <inputs> | |
| 66 | |
| 67 <conditional name="paired"> | |
|
7
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
68 <!-- <param name="is_paired" type="boolean" checked="true" truevalue="paired" falsevalue="single" label="Paired end reads?" help="Do you have 2 separate paired end read files?"/> --> |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
69 <param name="is_paired" type="select" label="Read type selector" help="single file (not pair aware), two files (paired reads) or a dataset collection of pairs"> |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
70 <option value="single" selected="true">Single File</option> |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
71 <option value="paired">Two Files</option> |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
72 <option value="collection">Single dataset collection pair</option> |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
73 </param> |
| 1 | 74 <when value="paired"> |
| 75 <param name="forwards_file" type="data" format="fastq,fastqsanger,fastqsolexa,fastqillumina" label="Direction 1 fastq reads to trim"/> | |
| 76 <param name="reverse_file" type="data" format="fastq,fastqsanger,fastqsolexa,fastqillumina" label="Direction 2 fastq reads to trim"/> | |
| 77 </when> | |
| 78 <when value="single"> | |
| 79 <param name="forwards_file" type="data" format="fastq,fastqsanger,fastqsolexa,fastqillumina" label="Fastq reads to trim"/> | |
| 80 </when> | |
|
7
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
81 <when value="collection"> |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
82 <param name="collection_data" type="data_collection" format="fastq,fastqsanger,fastqsolexa,fastqillumina" label="Paired end dataset collection" collection_type="paired"/> |
|
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
83 </when> |
| 1 | 84 </conditional> |
| 85 <param name="phred" type="select" label="Quality encoding." help="Phred33 or Phred 64, probably Phred64"> | |
| 86 <option value="phred64">phred64</option> | |
| 87 <option value="phred33" selected="True">phred33</option> | |
| 88 </param> | |
| 89 | |
| 90 <conditional name="adapt"> | |
| 91 <param name="adapters" type="boolean" checked="false" falsevalue="False" truevalue="True" label="Clip Illumina adapters?" help="Cut adapter and other illumina-specific sequences from the read."/> | |
| 92 <when value="True"> | |
| 93 <param name="adapt_file" type="data" format="fasta" label="Fasta of adapters to clip" help="Please supply a list of adapters to clip"/> | |
| 94 <param name="adapt_seed" type="integer" value="2" label="Seed mismatches" help="Specifies the maximum mismatch count which will still allow a full match to be performed."/> | |
| 95 <param name="adapt_palindrome" type="integer" value="40" label="Palindrome clip threshold" help="Specifies how accurate the match between the two 'adapter ligated' reads must be for PE palindrome read alignment."/> | |
| 96 <param name="adapt_simple" type="integer" value="15" label="Simple clip threshold" help="Specifies how accurate the match between any adapter etc. sequence must be against a read."/> | |
| 97 </when> | |
| 98 <when value="False"> | |
| 99 </when> | |
| 100 </conditional> | |
| 101 | |
| 102 <conditional name="sliding"> | |
| 103 <param name="slidingwindow" type="boolean" checked="true" falsevalue="False" truevalue="True" label="Perform Sliding Window trimming?" help="Perform a sliding window trimming, cutting once the average quality within the window falls below a threshold."/> | |
| 104 <when value="True"> | |
| 105 <param name="slid_window" type="integer" value="4" label="Sliding window size" help="Specifies the number of bases to average across"/> | |
| 106 <param name="slid_qual" type="integer" value="15" label="Average quality required" help="Specifies the average quality required."/> | |
| 107 </when> | |
| 108 <when value="false"> | |
| 109 </when> | |
| 110 </conditional> | |
| 111 | |
| 112 <conditional name="leading"> | |
| 113 <param name="trimleading" type="boolean" checked="true" falsevalue="False" truevalue="True" label="Trim leading bases?" help="Cut bases off the start of a read, if below a threshold quality."/> | |
| 114 <when value="True"> | |
| 115 <param name="lead_qual" type="integer" value="3" label="Minimum quality" help="Specifies the minimum quality required to keep a base."/> | |
| 116 </when> | |
| 117 <when value="false"> | |
| 118 </when> | |
| 119 </conditional> | |
| 120 | |
| 121 <conditional name="trailing"> | |
| 122 <param name="trimtrailing" type="boolean" checked="true" falsevalue="False" truevalue="True" label="Trim trailing bases?" help="Cut bases off the end of a read, if below a threshold quality."/> | |
| 123 <when value="True"> | |
| 124 <param name="trail_qual" type="integer" value="3" label="Minimum quality" help="Specifies the minimum quality required to keep a base."/> | |
| 125 </when> | |
| 126 <when value="false"> | |
| 127 </when> | |
| 128 </conditional> | |
| 129 | |
| 130 <conditional name="crop"> | |
| 131 <param name="cropend" type="boolean" checked="false" falsevalue="False" truevalue="True" label="Crop reads?" help="Cut the read to a specified length."/> | |
| 132 <when value="True"> | |
| 133 <param name="crop_len" type="integer" value="0" label="Number of bases" help="The number of bases to keep, from the start of the read."/> | |
| 134 </when> | |
| 135 <when value="false"> | |
| 136 </when> | |
| 137 </conditional> | |
| 138 | |
| 139 <conditional name="headcrop"> | |
| 140 <param name="headcrop_s" type="boolean" checked="false" falsevalue="False" truevalue="True" label="Crop starts?" help="Cut the specified number of bases from the start of the read."/> | |
| 141 <when value="True"> | |
| 142 <param name="headcrop_len" type="integer" value="0" label="Number of bases" help="The number of bases to remove from the start of the read."/> | |
| 143 </when> | |
| 144 <when value="false"> | |
| 145 </when> | |
| 146 </conditional> | |
| 147 | |
| 148 <param name="minlen" type="integer" value="24" label="Minimum length read" help="Drop the read if it is below specified length"/> | |
| 149 | |
| 150 <param name="logfile" type="boolean" checked="false" falsevalue="False" truevalue="True" label="Turn on log?" help="Detailed log of trims on each read. (Very large file!)"/> | |
| 151 </inputs> | |
| 152 | |
| 153 <outputs> | |
| 154 <data name="log" format="tabular" label="${tool.name} on ${on_string}: Trimmomatic Log" hidden="True"> | |
| 155 </data> | |
| 6 | 156 <data name="fwdpairs" format="input" label="${tool.name} on ${on_string}: Dir1 trimmed pairs"> |
|
7
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
157 <filter>(paired['is_paired']=="single")</filter> |
| 1 | 158 </data> |
| 6 | 159 <data name="revpairs" format="input" label="${tool.name} on ${on_string}: Dir2 trimmed pairs"> |
|
7
76d1aebc44bf
Uploaded: added dataset collection support for paired end trimmomatic
simon-gladman
parents:
6
diff
changeset
|
160 <filter>(paired['is_paired'] == "paired" or paired['is_paired'] == "collection")</filter> |
| 1 | 161 </data> |
| 6 | 162 <data name="singles" format="input" label="${tool.name} on ${on_string}: trimmed reads"/> |
| 1 | 163 <data name="dummy_out" format="tabular" label="${tool.name} on ${on_string}: DEBUG OUTPUT" hidden="True"/> |
| 164 </outputs> | |
| 5 | 165 |
| 166 <stdio> | |
| 167 <exit_code range="1:" level="fatal" description="Trimmomatic error" /> | |
| 168 </stdio> | |
| 1 | 169 <help> |
| 170 ***Trimmomatic*** | |
| 171 | |
| 172 A flexible read trimming tool for Illumina NGS data | |
| 173 | |
| 174 | |
| 175 | |
| 176 Trimmomatic performs a variety of useful trimming tasks for illumina paired-end and single ended data.The selection of trimming steps and their associated parameters are supplied on the command line. | |
| 177 | |
| 178 The current trimming steps are: | |
| 179 | |
| 180 ILLUMINACLIP: Cut adapter and other illumina-specific sequences from the read. | |
| 181 | |
| 182 SLIDINGWINDOW: Perform a sliding window trimming, cutting once the average quality within the window falls below a threshold. | |
| 183 | |
| 184 LEADING: Cut bases off the start of a read, if below a threshold quality | |
| 185 | |
| 186 TRAILING: Cut bases off the end of a read, if below a threshold quality | |
| 187 | |
| 188 CROP: Cut the read to a specified length | |
| 189 | |
| 190 HEADCROP: Cut the specified number of bases from the start of the read | |
| 191 | |
| 192 MINLEN: Drop the read if it is below a specified length | |
| 193 | |
| 194 **The Adapter Fasta** | |
| 195 | |
| 196 Illumina adapter and other technical sequences are copyrighted by Illumina, which makes it illegal for us to make them available without a licence. As such, you need to get access to the illumina sequences from another source, and create an appropriate adapter FASTA. | |
| 197 | |
| 198 To make this fasta, you must first understand how it will be used. Trimmomatic uses two strategies for adapter trimming: Palindrome and Simple | |
| 199 | |
| 200 With 'simple' trimming, each adapter sequence is tested against the reads, and if a sufficiently accurate match is detected, the read is clipped appropriately. | |
| 201 | |
| 202 'Palindrome' trimming is specifically designed for the case of 'reading through' a short fragment into the adapter sequence on the other end. In this approach, the appropriate adapter sequences are 'in silico ligated' onto the start of the reads, and the combined adapter+read sequences, forward and reverse are aligned. If they align in a manner which indicates 'read-through', the forward read is clipped and the reverse read dropped (since it contains no new data). | |
| 203 | |
| 204 Naming of the sequences indicates how they should be used. For 'Palindrome' clipping, the sequence names should both start with 'Prefix', and end in '/1' for the forward adapter and '/2' for the reverse adapter.All other sequences are checked using 'simple' mode. Sequences with names ending in '/1' or '/2' will be checked only against the forward or reverse read. All other sequences will be checked against both the forward and reverse read. If you want to check for the reverse-complement of a specific sequence, you need to specifically include the reverse-complemented form of the sequence as well, with another name. | |
| 205 | |
| 206 The thresholds used are a simplified log-likelihood approach. Each matching base adds just over 0.6, while each mismatch reduces the alignment score by Q/10. Therefore, a perfect match of a 20 base sequence will score just over 12, while 25 bases are needed to score 15. As such we recommend values between 12 - 15 for this parameter. For palindromic matches, the entire read sequence plus (partial) adapter sequences can be used - therefore this threshold can be higher, in the range of 30-40. The 'seed mismatch' parameter is used to make alignments more efficient, specifying the maximum base mismatch count in the 'seed' (16 bases). Typical values here are 1 or 2. | |
| 207 | |
| 208 **Citation** | |
| 209 | |
| 210 Lohse M, Bolger AM, Nagel A, Fernie AR, Lunn JE, Stitt M, Usadel B. RobiNA: a | |
| 211 user-friendly, integrated software solution for RNA-Seq-based transcriptomics. | |
| 212 Nucleic Acids Res. 2012 Jul;40(Web Server issue):W622-7. | |
| 213 | |
| 214 | |
| 215 | |
| 216 **Contacts** | |
| 217 | |
| 218 Anthony Bolger bolger at mpimp-golm dot mpg dot de | |
| 219 | |
| 220 Federico M. Giorgi federico dot giorgi at gmail dot com | |
| 221 | |
| 222 Wrapped for Galaxy by: | |
| 223 | |
| 224 Simon Gladman (VBC/VLSCI) simon dot gladman at monash dot edu | |
| 225 </help> | |
| 226 </tool> |