Mercurial > repos > devteam > picard
annotate picard_MergeBamAlignment.xml @ 12:05087b27692a draft
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 7491208ca0c917a053798a48c3e54c3e30e95d92
| author | devteam |
|---|---|
| date | Sun, 27 Nov 2016 15:11:50 -0500 |
| parents | 5eaa8a968300 |
| children | 7e6fd3d0f16e |
| rev | line source |
|---|---|
|
8
3a3234d7a2e8
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 00a7926c285bc4a339bd7deebf40b28f39c7d947-dirty
devteam
parents:
5
diff
changeset
|
1 <tool name="MergeBamAlignment" id="picard_MergeBamAlignment" version="@TOOL_VERSION@.0"> |
| 5 | 2 <description>merge alignment data with additional info stored in an unmapped BAM dataset</description> |
| 3 <macros> | |
| 4 <import>picard_macros.xml</import> | |
| 5 </macros> | |
|
8
3a3234d7a2e8
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 00a7926c285bc4a339bd7deebf40b28f39c7d947-dirty
devteam
parents:
5
diff
changeset
|
6 <expand macro="requirements" /> |
|
12
05087b27692a
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 7491208ca0c917a053798a48c3e54c3e30e95d92
devteam
parents:
9
diff
changeset
|
7 <command detect_errors="exit_code"><![CDATA[ |
| 5 | 8 @java_options@ |
| 9 #set $picard_dict = "localref.dict" | |
| 10 #set $ref_fasta = "localref.fa" ## This is done because picards "likes" .fa extension | |
| 11 | |
|
12
05087b27692a
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 7491208ca0c917a053798a48c3e54c3e30e95d92
devteam
parents:
9
diff
changeset
|
12 ln -s "${reference_source.ref_file}" "${ref_fasta}" && |
| 5 | 13 |
| 14 #if str( $reference_source.reference_source_selector ) == "history": | |
| 15 | |
|
12
05087b27692a
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 7491208ca0c917a053798a48c3e54c3e30e95d92
devteam
parents:
9
diff
changeset
|
16 picard CreateSequenceDictionary REFERENCE="${ref_fasta}" OUTPUT="${picard_dict}" |
| 5 | 17 QUIET=true |
| 18 VERBOSITY=ERROR | |
| 19 | |
|
12
05087b27692a
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 7491208ca0c917a053798a48c3e54c3e30e95d92
devteam
parents:
9
diff
changeset
|
20 && |
| 5 | 21 |
| 22 #else: | |
| 23 | |
| 24 #set $ref_fasta = str( $reference_source.ref_file.fields.path ) | |
| 25 | |
| 26 #end if | |
| 27 | |
|
12
05087b27692a
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 7491208ca0c917a053798a48c3e54c3e30e95d92
devteam
parents:
9
diff
changeset
|
28 picard |
| 5 | 29 MergeBamAlignment |
| 30 UNMAPPED_BAM="${unmapped_bam}" | |
| 31 | |
| 32 PAIRED_RUN=true ##This argument is ignored and will be removed. Required. Possible values: {true, false} | |
| 33 | |
| 34 #if str( $aligned_or_read1_and_read2.aligned_or_read1_and_read2_selector ) == "paired_one_file": | |
| 35 #for $dataset in $aligned_or_read1_and_read2.aligned_bams: | |
| 36 ALIGNED_BAM="${dataset.aligned_bam}" | |
| 37 #end for | |
| 38 #elif str( $aligned_or_read1_and_read2.aligned_or_read1_and_read2_selector ) == "paired_two_files": | |
| 39 #for $dataset in $aligned_or_read1_and_read2.read1_aligned_bams: | |
| 40 READ1_ALIGNED_BAM="${dataset.read1_aligned_bam}" | |
| 41 #end for | |
| 42 #for $dataset in $aligned_or_read1_and_read2.read2_aligned_bams: | |
|
8
3a3234d7a2e8
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 00a7926c285bc4a339bd7deebf40b28f39c7d947-dirty
devteam
parents:
5
diff
changeset
|
43 READ2_ALIGNED_BAM="${dataset.read2_aligned_bam}" |
| 5 | 44 #end for |
| 45 #else | |
| 46 #for $dataset in $aligned_or_read1_and_read2.read1_aligned_bams: | |
| 47 READ1_ALIGNED_BAM="${dataset.read1_aligned_bam}" | |
| 48 #end for | |
| 49 #end if | |
| 50 | |
| 51 OUTPUT="${outFile}" | |
| 52 REFERENCE_SEQUENCE="${ref_fasta}" | |
| 53 | |
| 54 CLIP_ADAPTERS="${clip_adapters}" | |
| 55 IS_BISULFITE_SEQUENCE="${is_bisulfite_sequence}" | |
| 56 ALIGNED_READS_ONLY="${aligned_reads_only}" | |
| 57 MAX_INSERTIONS_OR_DELETIONS="${max_insertions_or_deletions}" | |
| 58 | |
| 59 #for $attribute in $attributes_to_retain: | |
| 60 ATTRIBUTES_TO_RETAIN="${$attribute.attribute}" | |
| 61 #end for | |
| 62 | |
| 63 #for $attribute in $attributes_to_remove: | |
| 64 ATTRIBUTES_TO_REMOVE="${$attribute.attribute}" | |
| 65 #end for | |
| 66 | |
| 67 READ1_TRIM="${read1_trim}" | |
| 68 READ2_TRIM="${read2_trim}" | |
| 69 | |
| 70 #if str( $orientations ) != "None": | |
| 71 #for $orientation in str( $orientations ).split(','): ## See trello card https://trello.com/c/9nW02Zhd | |
| 72 EXPECTED_ORIENTATIONS="${orientation}" | |
| 73 #end for | |
| 74 #end if | |
| 75 | |
| 76 ALIGNER_PROPER_PAIR_FLAGS="${aligner_proper_pair_flags}" | |
| 77 PRIMARY_ALIGNMENT_STRATEGY="${primary_alignment_strategy}" | |
| 78 CLIP_OVERLAPPING_READS="${clip_overlapping_reads}" | |
| 79 INCLUDE_SECONDARY_ALIGNMENTS="${include_secondary_alignments}" | |
| 80 ADD_MATE_CIGAR="${add_mate_cigar}" | |
| 81 | |
| 82 VALIDATION_STRINGENCY="${validation_stringency}" | |
| 83 | |
| 84 SORT_ORDER=coordinate | |
| 85 QUIET=true | |
| 86 VERBOSITY=ERROR | |
| 87 | |
|
12
05087b27692a
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 7491208ca0c917a053798a48c3e54c3e30e95d92
devteam
parents:
9
diff
changeset
|
88 ]]></command> |
| 5 | 89 |
| 90 <inputs> | |
| 91 | |
| 92 <conditional name="reference_source"> | |
| 93 <param name="reference_source_selector" type="select" label="Load reference genome from"> | |
| 94 <option value="cached">Local cache</option> | |
| 95 <option value="history">History</option> | |
| 96 </param> | |
| 97 <when value="cached"> | |
| 98 <param name="ref_file" type="select" label="Use dictionary from the list" help="Select genome from the list"> | |
| 99 <options from_data_table="picard_indexes"> | |
| 100 <filter type="sort_by" column="2" /> | |
| 101 <validator type="no_options" message="No indexes are available" /> | |
| 102 </options> | |
| 103 <validator type="no_options" message="A built-in dictionary is not available for the build associated with the selected input file"/> | |
| 104 </param> | |
| 105 </when> | |
| 106 <when value="history"> | |
| 107 <param name="ref_file" type="data" format="fasta" label="Use the following dataset to create dictionary" help="You can upload a FASTA sequence to the history from which Picard will automatically generate dictionary using CreateSequenceDictionary command" /> | |
| 108 </when> | |
| 109 </conditional> | |
| 110 | |
| 111 <param format="sam,bam" name="unmapped_bam" type="data" label="Selected unaligned SAM or BAM with original reads" help="UNMAPPED_BAM; This dataset must be sorted in queryname order (use picard_SortSam to do this)" /> | |
| 112 <conditional name="aligned_or_read1_and_read2"> | |
| 113 <param name="aligned_or_read1_and_read2_selector" type="select" label="What type of aligned data do you have?"> | |
| 114 <option value="paired_one_file">Paired data in a single BAM file (ALIGNED_BAM)</option> | |
| 115 <option value="paired_two_files">Paired data in separate files (READ1_ALIGNED_BAM and READ2_ALIGNED_BAM)</option> | |
| 116 <option value="single_file">Singe end data (READ1_ALIGNED_BAM)</option> | |
| 117 </param> | |
| 118 <when value="paired_one_file"> | |
| 119 <repeat name="aligned_bams" title="Aligned SAM or BAM dataset" min="1" help="You can provide multiple datasets"> | |
| 120 <param name="aligned_bam" type="data" format="sam,bam" label="SAM or BAM dataset(s) with alignment data" help="ALIGNED_BAM"/> | |
| 121 </repeat> | |
| 122 </when> | |
| 123 <when value="paired_two_files"> | |
| 124 <repeat name="read1_aligned_bams" title="Aligned SAM or BAM dataset for Read 1" min="1" help="You can provide multiple datasets"> | |
| 125 <param name="read1_aligned_bam" type="data" format="sam,bam" label="SAM or BAM dataset(s) with alignment data for Read1" help="READ1_ALIGNED_BAM"/> | |
| 126 </repeat> | |
| 127 <repeat name="read2_aligned_bams" title="Aligned SAM or BAM dataset for Read 2" min="1" help="You can provide multiple datasets"> | |
| 128 <param name="read2_aligned_bam" type="data" format="sam,bam" label="SAM or BAM dataset(s) with alignment data for Read2" help="READ2_ALIGNED_BAM"/> | |
| 129 </repeat> | |
| 130 </when> | |
| 131 <when value="single_file"> | |
| 132 <repeat name="read1_aligned_bams" title="Aligned SAM or BAM dataset for Single Reads" min="1" help="You can provide multiple datasets"> | |
| 133 <param name="read1_aligned_bam" type="data" format="sam,bam" label="SAM or BAM dataset(s) with alignment data" help="READ1_ALIGNED_BAM"/> | |
| 134 </repeat> | |
| 135 </when> | |
| 136 </conditional> | |
| 137 | |
| 138 <param name="clip_adapters" type="boolean" checked="true" label="Whether to clip adapters where identified" help="CLIP_ADAPTERS; default=True"/> | |
| 139 <param name="is_bisulfite_sequence" type="boolean" label="Whether the data is from bisulfite sequencing (used when caculating the NM tag)" help="IS_BISULFITE_SEQUENCE; default=False"/> | |
| 140 <param name="aligned_reads_only" type="boolean" label="Whether to output only aligned reads" help="ALIGNED_READS_ONLY; default=False"/> | |
| 141 <param name="max_insertions_or_deletions" type="integer" value="1" label="The maximum number of insertions or deletions permitted for an alignment to be included" help="MAX_INSERTIONS_OR_DELETIONS; Alignments with more than this many insertions or deletions will be ignored. Set to -1 to allow any number of insertions or deletions. default=1"/> | |
| 142 | |
| 143 <repeat name="attributes_to_retain" title="Retain the following alignment attribute" min="0" help="You can provide multiple attributes"> | |
|
9
5eaa8a968300
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 1869970193a1878acbc0f8a79b81dd02b37f1dc1
devteam
parents:
8
diff
changeset
|
144 <param name="attribute" type="text" label="Reserved alignment attributes (tags starting with X, Y, or Z) that should be brought over from the alignment data when merging" help="ATTRIBUTES_TO_RETAIN; example: XA"/> |
| 5 | 145 </repeat> |
| 146 | |
| 147 <repeat name="attributes_to_remove" title="Remove the following alignment attribute" min="0" help="You can provide multiple attributes"> | |
|
9
5eaa8a968300
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 1869970193a1878acbc0f8a79b81dd02b37f1dc1
devteam
parents:
8
diff
changeset
|
148 <param name="attribute" type="text" label="Attributes from the alignment record that should be removed when merging." help="ATTRIBUTES_TO_REMOVE; This overrides ATTRIBUTES_TO_RETAIN if they share common tags"/> |
| 5 | 149 </repeat> |
| 150 | |
| 151 <param name="read1_trim" type="integer" value="0" label="The number of bases trimmed from the beginning of read 1 prior to alignment" help="READ1_TRIM; default=0"/> | |
| 152 <param name="read2_trim" type="integer" value="0" label="The number of bases trimmed from the beginning of read 2 prior to alignment" help="READ2_TRIM; default=0"/> | |
| 153 | |
| 154 <param name="orientations" type="select" multiple="True" display="checkboxes" label="The expected orientation of proper read pairs" help="EXPECTED_ORIENTATIONS; multiple orinetations can be selected"> | |
| 155 <option value="FR">Forward/Reverse (FR)</option> | |
| 156 <option value="RF">Reverse/Forward (RF)</option> | |
| 157 <option value="TANDEM">Tandem</option> | |
| 158 </param> | |
| 159 | |
| 160 | |
| 161 <param name="aligner_proper_pair_flags" type="boolean" label="Use the aligner's idea of what a proper pair is rather than computing in this program" help="ALIGNER_PROPER_PAIR_FLAGS; default=False"/> | |
| 162 | |
| 163 <param name="primary_alignment_strategy" type="select" label="Strategy for selecting primary alignment when the aligner has provided more than one alignment for a pair or fragments" help="PRIMARY_ALIGNMENT_STRATEGY; see help below for more info; default=BestMapq"> | |
| 164 <option value="BestMapq" selected="True">BestMapq</option> | |
| 165 <option value="EarliestFragment">EarliestFragment</option> | |
| 166 <option value="BestEndMapq">BestEndMapq</option> | |
| 167 <option value="MostDistant">MostDistant</option> | |
| 168 </param> | |
| 169 | |
| 170 <param name="clip_overlapping_reads" type="boolean" checked="True" label="For paired reads, soft clip the 3' end of each read if necessary so that it does not extend past the 5' end of its mate" help="CLIP_OVERLAPPING_READS; default=True"/> | |
| 171 <param name="include_secondary_alignments" type="boolean" checked="True" label="If false, do not write secondary alignments to output" help="INCLUDE_SECONDARY_ALIGNMENTS; default=True"/> | |
| 172 <param name="add_mate_cigar" type="boolean" checked="True" label="Adds the mate CIGAR tag (MC) if true, does not if false" help="ADD_MATE_CIGAR; default=True"/> | |
| 173 <expand macro="VS" /> | |
| 174 </inputs> | |
| 175 <outputs> | |
| 176 <data name="outFile" format="bam" label="${tool.name} on ${on_string}: BAM with merged alignments"/> | |
| 177 </outputs> | |
| 178 <tests> | |
| 179 <test> | |
| 180 <param name="reference_source_selector" value="history" /> | |
| 181 <param name="ref_file" value="picard_MergeBamAlignment_ref.fa" ftype="fasta" /> | |
| 182 <param name="unmapped_bam" value="picard_MergeBamAlignment_unaligned.bam" ftype="bam"/> | |
| 183 <param name="aligned_or_read1_and_read2_selector" value="paired_one_file"/> | |
| 184 <param name="aligned_bam" value="picard_MergeBamAlignment_aligned.bam" ftype="bam"/> | |
| 185 <param name="clip_adapters" value="True"/> | |
| 186 <param name="is_bisulfite_sequence" value="False"/> | |
| 187 <param name="aligned_reads_only" value="False"/> | |
| 188 <param name="max_insertions_or_deletions" value="1"/> | |
| 189 <param name="read1_trim" value="0"/> | |
| 190 <param name="read2_trim" value="0"/> | |
| 191 <param name="orientation" value="FR"/> | |
| 192 <param name="aligner_proper_pair_flags" value="False"/> | |
| 193 <param name="primary_alignment_strategy" value="BestMapq"/> | |
| 194 <param name="clip_overlapping_reads" value="True"/> | |
| 195 <param name="include_secondary_alignments" value="True"/> | |
| 196 <param name="add_mate_cigar" value="True"/> | |
|
12
05087b27692a
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/picard commit 7491208ca0c917a053798a48c3e54c3e30e95d92
devteam
parents:
9
diff
changeset
|
197 <output name="outFile" file="picard_MergeBamAlignment_test1.bam" ftype="bam" lines_diff="4"/> |
| 5 | 198 </test> |
| 199 </tests> | |
| 200 | |
| 201 | |
| 202 <help> | |
| 203 | |
| 204 .. class:: infomark | |
| 205 | |
| 206 **Purpose** | |
| 207 | |
| 208 Merges alignment data from a SAM or BAM dataset with additional data stored in an unmapped BAM dataset and produces a third SAM or BAM dataset of aligned and unaligned reads. | |
| 209 | |
| 210 @dataset_collections@ | |
| 211 | |
| 212 @description@ | |
| 213 | |
| 214 UNMAPPED_BAM=File | |
| 215 UNMAPPED=File Original SAM or BAM file of unmapped reads, which must be in queryname order. Required. | |
| 216 | |
| 217 ALIGNED_BAM=File | |
| 218 ALIGNED=File SAM or BAM file(s) with alignment data. This option may be specified 0 or more times. | |
| 219 Cannot be used in conjuction with option(s) READ1_ALIGNED_BAM (R1_ALIGNED) | |
| 220 READ2_ALIGNED_BAM (R2_ALIGNED) | |
| 221 | |
| 222 READ1_ALIGNED_BAM=File | |
| 223 R1_ALIGNED=File SAM or BAM file(s) with alignment data from the first read of a pair. This option may be | |
| 224 specified 0 or more times. Cannot be used in conjuction with option(s) ALIGNED_BAM | |
| 225 (ALIGNED) | |
| 226 | |
| 227 READ2_ALIGNED_BAM=File | |
| 228 R2_ALIGNED=File SAM or BAM file(s) with alignment data from the second read of a pair. This option may | |
| 229 be specified 0 or more times. Cannot be used in conjuction with option(s) ALIGNED_BAM | |
| 230 (ALIGNED) | |
| 231 | |
| 232 PAIRED_RUN=Boolean | |
| 233 PE=Boolean This argument is ignored and will be removed. Required. Possible values: {true, false} | |
| 234 | |
| 235 JUMP_SIZE=Integer | |
| 236 JUMP=Integer The expected jump size (required if this is a jumping library). Deprecated. Use | |
| 237 EXPECTED_ORIENTATIONS instead Default value: null. Cannot be used in conjuction with | |
| 238 option(s) EXPECTED_ORIENTATIONS (ORIENTATIONS) | |
| 239 | |
| 240 CLIP_ADAPTERS=Boolean Whether to clip adapters where identified. Default value: true. Possible values: {true, false} | |
| 241 | |
| 242 IS_BISULFITE_SEQUENCE=Boolean Whether the lane is bisulfite sequence (used when caculating the NM tag). Default value: | |
| 243 false. Possible values: {true, false} | |
| 244 | |
| 245 ALIGNED_READS_ONLY=Boolean Whether to output only aligned reads. Default value: false. Possible values: {true, false} | |
| 246 | |
| 247 MAX_INSERTIONS_OR_DELETIONS=Integer | |
| 248 MAX_GAPS=Integer The maximum number of insertions or deletions permitted for an alignment to be included. | |
| 249 Alignments with more than this many insertions or deletions will be ignored. Set to -1 to | |
| 250 allow any number of insertions or deletions. Default value: 1. | |
| 251 | |
| 252 ATTRIBUTES_TO_RETAIN=String Reserved alignment attributes (tags starting with X, Y, or Z) that should be brought over | |
| 253 from the alignment data when merging. This option may be specified 0 or more times. | |
| 254 | |
| 255 ATTRIBUTES_TO_REMOVE=String Attributes from the alignment record that should be removed when merging. This overrides | |
| 256 ATTRIBUTES_TO_RETAIN if they share common tags. This option may be specified 0 or more | |
| 257 times. | |
| 258 | |
| 259 READ1_TRIM=Integer | |
| 260 R1_TRIM=Integer The number of bases trimmed from the beginning of read 1 prior to alignment Default | |
| 261 value: 0. | |
| 262 | |
| 263 READ2_TRIM=Integer | |
| 264 R2_TRIM=Integer The number of bases trimmed from the beginning of read 2 prior to alignment Default | |
| 265 value: 0. | |
| 266 | |
| 267 EXPECTED_ORIENTATIONS=PairOrientation | |
| 268 ORIENTATIONS=PairOrientation The expected orientation of proper read pairs. Replaces JUMP_SIZE Possible values: {FR, | |
| 269 RF, TANDEM} This option may be specified 0 or more times. Cannot be used in conjuction | |
| 270 with option(s) JUMP_SIZE (JUMP) | |
| 271 | |
| 272 ALIGNER_PROPER_PAIR_FLAGS=Boolean | |
| 273 Use the aligner's idea of what a proper pair is rather than computing in this program. | |
| 274 Default value: false. Possible values: {true, false} | |
| 275 | |
| 276 SORT_ORDER=SortOrder | |
| 277 SO=SortOrder The order in which the merged reads should be output. Default value: coordinate. | |
| 278 Possible values: {unsorted, queryname, coordinate} | |
| 279 | |
| 280 PRIMARY_ALIGNMENT_STRATEGY=PrimaryAlignmentStrategy | |
| 281 Strategy for selecting primary alignment when the aligner has provided more than one | |
| 282 alignment for a pair or fragment, and none are marked as primary, more than one is marked | |
| 283 as primary, or the primary alignment is filtered out for some reason. BestMapq expects | |
| 284 that multiple alignments will be correlated with HI tag, and prefers the pair of | |
| 285 alignments with the largest MAPQ, in the absence of a primary selected by the aligner. | |
| 286 EarliestFragment prefers the alignment which maps the earliest base in the read. Note | |
| 287 that EarliestFragment may not be used for paired reads. BestEndMapq is appropriate for | |
| 288 cases in which the aligner is not pair-aware, and does not output the HI tag. It simply | |
| 289 picks the alignment for each end with the highest MAPQ, and makes those alignments | |
| 290 primary, regardless of whether the two alignments make sense together.MostDistant is also | |
| 291 for a non-pair-aware aligner, and picks the alignment pair with the largest insert size. | |
| 292 If all alignments would be chimeric, it picks the alignments for each end with the best | |
| 293 MAPQ. For all algorithms, ties are resolved arbitrarily. Default value: BestMapq. | |
| 294 Possible values: {BestMapq, EarliestFragment, BestEndMapq, MostDistant} | |
| 295 | |
| 296 CLIP_OVERLAPPING_READS=BooleanFor paired reads, soft clip the 3' end of each read if necessary so that it does not | |
| 297 extend past the 5' end of its mate. Default value: true. Possible values: {true, false} | |
| 298 | |
| 299 INCLUDE_SECONDARY_ALIGNMENTS=Boolean | |
| 300 If false, do not write secondary alignments to output. Default value: true. | |
| 301 Possible values: {true, false} | |
| 302 | |
| 303 ADD_MATE_CIGAR=Boolean | |
| 304 MC=Boolean Adds the mate CIGAR tag (MC) if true, does not if false. Possible values: {true, false} | |
| 305 | |
| 306 | |
| 307 | |
| 308 | |
| 309 @more_info@ | |
| 310 </help> | |
| 311 </tool> | |
| 312 | |
| 313 | |
| 314 | |
| 315 | |
| 316 | |
| 317 | |
| 318 | |
| 319 | |
| 320 | |
| 321 | |
| 322 | |
| 323 | |
| 324 |