|
140
|
1 <tool name="CollectRnaSeqMetrics" id="picard_CollectRnaSeqMetrics" version="1.106.0">
|
|
|
2 <description>Collect RNA-Seq Metrics</description>
|
|
|
3 <requirements><requirement type="package" version="1.106.0">picard</requirement></requirements>
|
|
|
4 <command interpreter="python">
|
|
|
5 picard_wrapper.py -i "${input_file}" -d "${html_file.files_path}" -t "${html_file}"
|
|
|
6 -n "${out_prefix}" --tmpdir "${__new_file_path__}" --assumesorted ${ASSUME_SORTED}
|
|
|
7 --refflat ${REF_FLAT}
|
|
|
8 #if $identify_ribosomal.opt == "yes"
|
|
|
9 --ribosomalintervals ${identify_ribosomal.RIBOSOMAL_INTERVALS}
|
|
|
10 #end if
|
|
|
11 --malevel "${malevel}"
|
|
|
12 --minlength ${MINIMUM_LENGTH}
|
|
|
13 --strandspecificity ${STRAND_SPECIFICITY}
|
|
|
14 --rrnafragmentpercentage ${RRNA_FRAGMENT_PERCENTAGE}
|
|
|
15 #for $i in $IGNORE_SEQUENCES
|
|
|
16 --ignoreseq "${i.IGNORE_SEQUENCE}"
|
|
|
17 #end for
|
|
|
18 -j "\$JAVA_JAR_PATH/CollectRnaSeqMetrics.jar"
|
|
|
19 </command>
|
|
|
20
|
|
|
21 <stdio>
|
|
|
22 <exit_code range="0" level="warning" description="Tool finished correctly" />
|
|
|
23 </stdio>
|
|
|
24
|
|
|
25 <inputs>
|
|
|
26 <param format="sam" name="input_file" type="data" label="Input SAM file." help="" />
|
|
|
27 <param name="out_prefix" value="RNA-Seq Metrics" type="text"
|
|
|
28 label="Title for the output file" help="Use this remind you what the job was for." size="80" />
|
|
|
29
|
|
|
30 <param format="data" name="REF_FLAT" type="data" label="Gene annotations in refFlat form. Format described here: http://genome.ucsc.edu/goldenPath/gbdDescriptionsOld.html#RefFlat Required." help="" />
|
|
|
31
|
|
|
32 <conditional name="identify_ribosomal">
|
|
|
33 <param name="opt" type="select" label="Identify ribosomal bases" help="If 'no' is selected, no bases will be identified as being ribosomal.">
|
|
|
34 <option value="no">no</option>
|
|
|
35 <option value="yes">yes</option>
|
|
|
36 </param>
|
|
|
37 <when value="no" />
|
|
|
38 <when value="yes">
|
|
|
39 <param format="data" name="RIBOSOMAL_INTERVALS" type="data" label="Location of rRNA sequences in genome, in interval_list format. If not specified no bases will be identified as being ribosomal. Format described here: http://picard.sourceforge.net/javadoc/net/sf/picard/util/IntervalList.html Default value: null." help="" />
|
|
|
40 </when>
|
|
|
41 </conditional>
|
|
|
42
|
|
|
43 <param name="STRAND_SPECIFICITY" type="select" label="For strand-specific library prep." help="For unpaired reads, use FIRST_READ_TRANSCRIPTION_STRAND if the reads are expected to be on the transcription strand.">
|
|
|
44 <option value="NONE" selected="True">None</option>
|
|
|
45 <option value="FIRST_READ_TRANSCRIPTION_STRAND">FIRST_READ_TRANSCRIPTION_STRAND</option>
|
|
|
46 <option value="SECOND_READ_TRANSCRIPTION_STRAND">SECOND_READ_TRANSCRIPTION_STRAND</option>
|
|
|
47 </param>
|
|
|
48
|
|
|
49 <param name="MINIMUM_LENGTH" type="text" value="500" label="When calculating coverage based values (e.g. CV of coverage) only use transcripts of this length or greater." help="" />
|
|
|
50
|
|
|
51 <repeat name="IGNORE_SEQUENCES" title="Ignore Sequences">
|
|
|
52 <param name="IGNORE_SEQUENCE" label="Ignore Sequence" type="text" help="If a read maps to a sequence specified with this option, all the bases in the read are counted as ignored bases." />
|
|
|
53 </repeat>
|
|
|
54
|
|
|
55 <param name="RRNA_FRAGMENT_PERCENTAGE" type="text" value="0.8" label="This percentage of the length of a fragment must overlap one of the ribosomal intervals for a read or read pair by this must in order to be considered rRNA." help="" />
|
|
|
56
|
|
|
57 <param name="malevel" value="0" type="select" multiple="true" label="Metric Accumulation Level"
|
|
|
58 help="Level(s) at which metrics will be accumulated">
|
|
|
59 <option value="ALL_READS" selected="true">All reads (default)</option>
|
|
|
60 <option value="SAMPLE" default="true">Sample</option>
|
|
|
61 <option value="LIBRARY" default="true">Library</option>
|
|
|
62 <option value="READ_GROUP" default="true">Read group</option>
|
|
|
63 </param>
|
|
|
64
|
|
|
65 <param checked="True" truevalue="true" falsevalue="false" name="ASSUME_SORTED" type="boolean" label="If true (default), then the sort order in the header file will be ignored." />
|
|
|
66 </inputs>
|
|
|
67 <outputs>
|
|
|
68 <data format="html" name="html_file" label="${out_prefix}.html"/>
|
|
|
69 </outputs>
|
|
|
70 <tests>
|
|
|
71 <test>
|
|
|
72 <!-- python picard_wrapper.py
|
|
|
73 -i "/home/~/PICARD-in.sam"
|
|
|
74 -d "/home/~/outputrnaseqùetrics_files"
|
|
|
75 -t "/home/~/outputrnaseqmetrics.htm"
|
|
|
76 -n "RNA-Seq Metrics"
|
|
|
77 --tmpdir "/home/dorine/galaxypicard/galaxy-central/database/tmp"
|
|
|
78 --assumesorted true
|
|
|
79 --refflat /home/~/refFlat.txt --malevel "ALL_READS" --minlength 500 --strandspecificity NONE
|
|
|
80 --rrnafragmentpercentage 0.8 -j "/home/~/CollectRnaSeqMetrics.jar" -->
|
|
|
81 <param name="inputFile" value="PICARD-in.sam" />
|
|
|
82 <param name="out_prefix" value="RNA-Seq Metrics" />
|
|
|
83 <param name="ASSUME_SORTED" value="true" />
|
|
|
84 <param name="REF_FLAT" value="refFlat.txt" />
|
|
|
85 <param name="malevel" value="ALL_READS" />
|
|
|
86 <param name="MINIMUM_LENGTH" value="500" />
|
|
|
87 <param name="STRAND_SPECIFICITY" value="NONE" />
|
|
|
88 <param name="RRNA_FRAGMENT_PERCENTAGE" value="O.8" />
|
|
|
89 <output name="html_file" file="outputrnaseqmetrics.html" ftype="html" lines_diff="30"/>
|
|
|
90 </test>
|
|
|
91 </tests>
|
|
|
92 <help>
|
|
|
93 Picard documentation says:
|
|
|
94
|
|
|
95
|
|
|
96 CollectRnaSeqMetrics
|
|
|
97
|
|
|
98 Documentation: http://picard.sourceforge.net/command-line-overview.shtml#CollectRnaSeqMetrics
|
|
|
99
|
|
|
100 Program to collect metrics about the alignment of RNA to various functional classes of loci in the genome: coding, intronic, UTR, intergenic, ribosomal. Also determines strand-specificity for strand-specific libraries.
|
|
|
101
|
|
|
102 </help>
|
|
|
103 </tool>
|
