annotate read_quality.xml @ 10:c63fbdcc0dbe draft

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 57f71aa633a43ab02bbf05acd0c6d7f406e01f1e"
author iuc
date Thu, 28 Nov 2019 20:46:57 +0000
parents caaa120457bc
children d7f6b3653d84
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
1 <tool id="rseqc_read_quality" name="Read Quality" version="@WRAPPER_VERSION@">
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
2 <description>determines Phred quality score</description>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
3
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
4 <macros>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
5 <import>rseqc_macros.xml</import>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
6 </macros>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
7
4
caaa120457bc planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 7f68686cac77df831f1a26a2126a238a2e480316
iuc
parents: 3
diff changeset
8 <expand macro="requirements">
caaa120457bc planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 7f68686cac77df831f1a26a2126a238a2e480316
iuc
parents: 3
diff changeset
9 <!--
caaa120457bc planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 7f68686cac77df831f1a26a2126a238a2e480316
iuc
parents: 3
diff changeset
10 Required due to conda solver bug: https://github.com/conda/conda/issues/6269
caaa120457bc planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 7f68686cac77df831f1a26a2126a238a2e480316
iuc
parents: 3
diff changeset
11 See: https://github.com/galaxyproject/tools-iuc/pull/1578 for more info
caaa120457bc planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 7f68686cac77df831f1a26a2126a238a2e480316
iuc
parents: 3
diff changeset
12 -->
caaa120457bc planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 7f68686cac77df831f1a26a2126a238a2e480316
iuc
parents: 3
diff changeset
13 <requirement type="package" version="3.4.1">r-base</requirement>
caaa120457bc planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 7f68686cac77df831f1a26a2126a238a2e480316
iuc
parents: 3
diff changeset
14 </expand>
3
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
15
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
16 <expand macro="stdio" />
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
17
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
18 <version_command><![CDATA[read_quality.py --version]]></version_command>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
19
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
20 <command><![CDATA[
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
21 read_quality.py
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
22 --input-file '${input}'
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
23 --out-prefix output
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
24 -r ${reduce}
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
25 --mapq ${mapq}
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
26 ]]>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
27 </command>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
28
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
29 <inputs>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
30 <expand macro="bam_sam_param" />
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
31 <param name="reduce" type="integer" label="Ignore Phred scores less than this amount (only applies to 'boxplot', default=1000)" value="1000" help="(--reduce)"/>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
32 <expand macro="mapq_param" />
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
33 <expand macro="rscript_output_param" />
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
34 </inputs>
1
f92b87abef3d just xmls
nilesh
parents:
diff changeset
35
3
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
36 <outputs>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
37 <data format="pdf" name="outputheatpdf" from_work_dir="output.qual.heatmap.pdf" label="${tool.name} on ${on_string} (Heatmap pdf)" />
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
38 <data format="pdf" name="outputboxpdf" from_work_dir="output.qual.boxplot.pdf" label="${tool.name} on ${on_string} (Boxplot pdf)" />
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
39 <expand macro="rscript_output_data" filename="output.qual.r" />
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
40 </outputs>
1
f92b87abef3d just xmls
nilesh
parents:
diff changeset
41
3
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
42 <tests>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
43 <test>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
44 <param name="input" value="pairend_strandspecific_51mer_hg19_random.bam"/>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
45 <param name="rscript_output" value="true" />
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
46 <output name="outputr" file="output.qual.r"/>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
47 <output name="outputheatpdf" file="output.qual.heatmap.pdf" compare="sim_size" />
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
48 <output name="outputboxpdf" file="output.qual.boxplot.pdf" compare="sim_size" />
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
49 </test>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
50 </tests>
1
f92b87abef3d just xmls
nilesh
parents:
diff changeset
51
3
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
52 <help><![CDATA[
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
53 read_quality.py
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
54 +++++++++++++++
1
f92b87abef3d just xmls
nilesh
parents:
diff changeset
55
3
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
56 According to SAM specification, if Q is the character to represent "base calling quality"
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
57 in SAM file, then Phred Quality Score = ord(Q) - 33. Here ord() is python function that
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
58 returns an integer representing the Unicode code point of the character when the argument
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
59 is a unicode object, for example, ord('a') returns 97. Phred quality score is widely used
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
60 to measure "reliability" of base-calling, for example, phred quality score of 20 means
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
61 there is 1/100 chance that the base-calling is wrong, phred quality score of 30 means there
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
62 is 1/1000 chance that the base-calling is wrong. In general: Phred quality score = -10xlog(10)P,
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
63 here P is probability that base-calling is wrong.
1
f92b87abef3d just xmls
nilesh
parents:
diff changeset
64
f92b87abef3d just xmls
nilesh
parents:
diff changeset
65 Inputs
f92b87abef3d just xmls
nilesh
parents:
diff changeset
66 ++++++++++++++
f92b87abef3d just xmls
nilesh
parents:
diff changeset
67
f92b87abef3d just xmls
nilesh
parents:
diff changeset
68 Input BAM/SAM file
3
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
69 Alignment file in BAM/SAM format.
1
f92b87abef3d just xmls
nilesh
parents:
diff changeset
70
f92b87abef3d just xmls
nilesh
parents:
diff changeset
71 Ignore phred scores less than this number (default=1000)
3
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
72 To avoid making huge vector in R, nucleotide with certain phred score represented less than this number will be ignored. Increase this number save more memory while reduce precision. This option only applies to the 'boxplot'.
1
f92b87abef3d just xmls
nilesh
parents:
diff changeset
73
f92b87abef3d just xmls
nilesh
parents:
diff changeset
74 Output
f92b87abef3d just xmls
nilesh
parents:
diff changeset
75 ++++++++++++++
f92b87abef3d just xmls
nilesh
parents:
diff changeset
76
f92b87abef3d just xmls
nilesh
parents:
diff changeset
77 1. output.qual.r
f92b87abef3d just xmls
nilesh
parents:
diff changeset
78 2. output.qual.boxplot.pdf
3
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
79 .. image:: $PATH_TO_IMAGES/36mer.qual.plot.png
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
80 :height: 600 px
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
81 :width: 600 px
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
82 :scale: 80 %
1
f92b87abef3d just xmls
nilesh
parents:
diff changeset
83 3. output.qual.heatmap.pdf
3
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
84 .. image:: $PATH_TO_IMAGES/36mer.qual.heatmap.png
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
85 :height: 600 px
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
86 :width: 600 px
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
87 :scale: 80 %
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
88
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
89 Heatmap: use different color to represent nucleotide density ("blue"=low density,"orange"=median density,"red"=high density")
1
f92b87abef3d just xmls
nilesh
parents:
diff changeset
90
3
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
91 @ABOUT@
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
92
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
93 ]]>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
94 </help>
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
95
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
96 <expand macro="citations" />
71ed55a3515a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 37fb1988971807c6a072e1afd98eeea02329ee83
iuc
parents: 1
diff changeset
97
1
f92b87abef3d just xmls
nilesh
parents:
diff changeset
98 </tool>