Mercurial > repos > devteam > picard

comparison picard_EstimateLibraryComplexity.xml @ 33:3f254c5ced1d draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/picard commit 9ecbbb878d68a980ba35a90865e524c723ca3ed8
author iuc
date Sun, 03 Mar 2024 16:06:11 +0000
parents f9242e01365a
children
comparison
equal deleted inserted replaced
32:f9242e01365a 33:3f254c5ced1d
1 <tool name="EstimateLibraryComplexity" id="picard_EstimateLibraryComplexity" version="@TOOL_VERSION@.@WRAPPER_VERSION@"> 1 <tool name="EstimateLibraryComplexity" id="picard_EstimateLibraryComplexity" version="@TOOL_VERSION@.@WRAPPER_VERSION@" profile="@PROFILE@">
2 <description>assess sequence library complexity from read sequences</description> 2 <description>assess sequence library complexity from read sequences</description>
3 <macros> 3 <macros>
4 <import>picard_macros.xml</import> 4 <import>picard_macros.xml</import>
5 <token name="@WRAPPER_VERSION@">1</token> 5 <token name="@WRAPPER_VERSION@">0</token>
6 </macros> 6 </macros>
7 <expand macro="requirements" /> 7 <expand macro="requirements"/>
8 <command detect_errors="exit_code"><![CDATA[ 8 <command detect_errors="exit_code"><![CDATA[
9 @java_options@ 9 @java_options@
10 @symlink_element_identifier@ 10 @symlink_element_identifier@
11 picard 11
12 EstimateLibraryComplexity 12 picard EstimateLibraryComplexity
13 13
14 INPUT='$escaped_element_identifier' 14 --INPUT '$escaped_element_identifier'
15 OUTPUT="${outFile}" 15 --OUTPUT '${outFile}'
16 16
17 MIN_IDENTICAL_BASES="${min_identical_bases}" 17 --MIN_IDENTICAL_BASES '${min_identical_bases}'
18 MAX_DIFF_RATE="${max_diff_rate}" 18 --MAX_DIFF_RATE '${max_diff_rate}'
19 MIN_MEAN_QUALITY="${min_mean_quality}" 19 --MIN_MEAN_QUALITY '${min_mean_quality}'
20 MAX_GROUP_RATIO="${max_group_ratio}" 20 --MAX_GROUP_RATIO '${max_group_ratio}'
21 READ_NAME_REGEX='${ str( $read_name_regex ) }' 21 --READ_NAME_REGEX '${ str( $read_name_regex ) }'
22 OPTICAL_DUPLICATE_PIXEL_DISTANCE="${optical_duplicate_pixel_distance}" 22 --OPTICAL_DUPLICATE_PIXEL_DISTANCE '${optical_duplicate_pixel_distance}'
23 23
24 VALIDATION_STRINGENCY="${validation_stringency}" 24 --VALIDATION_STRINGENCY '${validation_stringency}'
25 QUIET=true 25 --QUIET true
26 VERBOSITY=ERROR 26 --VERBOSITY ERROR
27 27
28 ]]></command> 28 ]]></command>
29 <inputs> 29 <inputs>
30 <param format="bam" name="inputFile" type="data" label="Select SAM/BAM dataset or dataset collection" help="If empty, upload or import a SAM/BAM dataset" /> 30 <param format="bam" name="inputFile" type="data" label="Select SAM/BAM dataset or dataset collection" help="If empty, upload or import a SAM/BAM dataset"/>
31 <param name="min_identical_bases" type="integer" value="5" label="The minimum number of bases at the starts of reads that must be identical for reads to be grouped together for duplicate detection" help="MIN_IDENTICAL_BASES; In effect total_reads / 4^max_id_bases reads will be compared at a time, so lower numbers will produce more accurate results but consume exponentially more memory and CPU; default=5"/> 31 <param name="min_identical_bases" type="integer" value="5" label="The minimum number of bases at the starts of reads that must be identical for reads to be grouped together for duplicate detection" help="MIN_IDENTICAL_BASES; In effect total_reads / 4^max_id_bases reads will be compared at a time, so lower numbers will produce more accurate results but consume exponentially more memory and CPU; default=5"/>
32 <param name="max_diff_rate" type="float" value="0.03" label="The maximum rate of differences between two reads to call them identical" help="MAX_DIFF_RATE; default=0.03"/> 32 <param name="max_diff_rate" type="float" value="0.03" label="The maximum rate of differences between two reads to call them identical" help="MAX_DIFF_RATE; default=0.03"/>
33 <param name="min_mean_quality" type="integer" min="0" max="93" value="20" label="The minimum mean quality of the bases in a read pair for the read to be analyzed" help="MIN_MEAN_QUALITY; Reads with lower average quality are filtered out and not considered in any calculations; default=20"/> 33 <param name="min_mean_quality" type="integer" min="0" max="93" value="20" label="The minimum mean quality of the bases in a read pair for the read to be analyzed" help="MIN_MEAN_QUALITY; Reads with lower average quality are filtered out and not considered in any calculations; default=20"/>
34 <param name="max_group_ratio" type="integer" value="500" label="Do not process self-similar groups that are this many times over the mean expected group size" help="MAX_GROUP_RATIO; I.e. if the input contains 10m read pairs and MIN_IDENTICAL_BASES is set to 5, then the mean expected group size would be approximately 10 reads; default-500"/> 34 <param name="max_group_ratio" type="integer" value="500" label="Do not process self-similar groups that are this many times over the mean expected group size" help="MAX_GROUP_RATIO; I.e. if the input contains 10m read pairs and MIN_IDENTICAL_BASES is set to 5, then the mean expected group size would be approximately 10 reads; default-500"/>
35 35 <param name="read_name_regex" type="text" value="[a-zA-Z0-9]+:[0-9]:([0-9]+):([0-9]+):([0-9]+).*." label="Regular expression that can be used to parse read names in the incoming SAM/BAM dataset" help="READ_NAME_REGEX; Read names are parsed to extract three variables: tile/region, x coordinate and y coordinate. These values are used to estimate the rate of optical duplication in order to give a more accurate estimated library size. See help below for more info; default=[a-zA-Z0-9]+:[0-9]:([0-9]+):([0-9]+):([0-9]+).*.">
36 <param name="read_name_regex" type="text" value="[a-zA-Z0-9]+:[0-9]:([0-9]+):([0-9]+):([0-9]+).*." label="Regular expression that can be used to parse read names in the incoming SAM/BAM dataset" help="READ_NAME_REGEX; Read names are parsed to extract three variables: tile/region, x coordinate and y coordinate. These values are used to estimate the rate of optical duplication in order to give a more accurate estimated library size. See help below for more info; default=[a-zA-Z0-9]+:[0-9]:([0-9]+):([0-9]+):([0-9]+).*."> 36 <expand macro="sanitize_query"/>
37 <expand macro="sanitize_query" /> 37 </param>
38 </param> 38 <param name="optical_duplicate_pixel_distance" type="integer" value="100" min="0" max="500" label="The maximum offset between two duplicte clusters in order to consider them optical duplicates" help="OPTICAL_DUPLICATE_PIXEL_DISTANCE; default=100"/>
39 <param name="optical_duplicate_pixel_distance" type="integer" value="100" min="0" max="500" label="The maximum offset between two duplicte clusters in order to consider them optical duplicates" help="OPTICAL_DUPLICATE_PIXEL_DISTANCE; default=100"/> 39 <expand macro="VS"/>
40 40 </inputs>
41 <expand macro="VS" /> 41 <outputs>
42 42 <data format="tabular" name="outFile" label="${tool.name} on ${on_string}: Library complexity report"/>
43 </inputs> 43 </outputs>
44 44 <tests>
45 <outputs> 45 <test>
46 <data format="tabular" name="outFile" label="${tool.name} on ${on_string}: Library complexity report"/> 46 <param name="inputFile" value="picard_EstimateLibraryComplexity.bam" ftype="bam"/>
47 </outputs> 47 <param name="min_identical_bases" value="5"/>
48 48 <param name="max_diff_rate" value="0.03"/>
49 <tests> 49 <param name="min_mean_quality" value="20"/>
50 <test> 50 <param name="read_name_regex" value="[a-zA-Z0-9]+:[0-9]:([0-9]+):([0-9]+):([0-9]+).*."/>
51 <param name="inputFile" value="picard_EstimateLibraryComplexity.bam" ftype="bam"/> 51 <param name="optical_duplicate_pixel_distance" value="100"/>
52 <param name="min_identical_bases" value="5"/> 52 <param name="max_group_ratio" value="500"/>
53 <param name="max_diff_rate" value="0.03"/> 53 <param name="validation_stringency" value="LENIENT"/>
54 <param name="min_mean_quality" value="20"/> 54 <output name="outFile" file="picard_EstimateLibraryComplexity_test1.tab" ftype="tabular" lines_diff="4"/>
55 <param name="read_name_regex" value="[a-zA-Z0-9]+:[0-9]:([0-9]+):([0-9]+):([0-9]+).*."/> 55 </test>
56 <param name="optical_duplicate_pixel_distance" value="100"/> 56 </tests>
57 <param name="max_group_ratio" value="500"/> 57 <help>
58 <param name="validation_stringency" value="LENIENT"/>
59 <output name="outFile" file="picard_EstimateLibraryComplexity_test1.tab" ftype="tabular" lines_diff="4"/>
60 </test>
61 </tests>
62
63
64 <help>
65 58
66 **Purpose** 59 **Purpose**
67 60
68 Attempts to estimate library complexity from sequence of read pairs alone. Does so by sorting all reads by the first N bases (5 by default) 61 Attempts to estimate library complexity from sequence of read pairs alone. Does so by sorting all reads by the first N bases (5 by default)
69 of each read and then comparing reads with the first N bases identical to each other for duplicates. Reads are considered to be duplicates 62 of each read and then comparing reads with the first N bases identical to each other for duplicates. Reads are considered to be duplicates
74 67
75 Unpaired reads are ignored in this computation. 68 Unpaired reads are ignored in this computation.
76 The algorithm attempts to detect optical duplicates separately from PCR duplicates and excludes these in the calculation of library size. 69 The algorithm attempts to detect optical duplicates separately from PCR duplicates and excludes these in the calculation of library size.
77 70
78 Also, since there is no alignment to screen out technical reads one further filter is applied on the data. After examining all reads a Histogram 71 Also, since there is no alignment to screen out technical reads one further filter is applied on the data. After examining all reads a Histogram
79 is built of [#reads in duplicate set -> #of duplicate sets]; all bins that contain exactly one duplicate set are then removed from the Histogram 72 is built of [#reads in duplicate set -&gt; #of duplicate sets]; all bins that contain exactly one duplicate set are then removed from the Histogram
80 as outliers before library size is estimated. 73 as outliers before library size is estimated.
81 74
82 @dataset_collections@ 75 @dataset_collections@
83 76
84 @description@ 77 @description@
119 112
120 113
121 @more_info@ 114 @more_info@
122 115
123 </help> 116 </help>
124 <expand macro="citations" /> 117 <expand macro="citations"/>
125 </tool> 118 </tool>