jpayne@0
|
1 <tool id="cfsan-qaqc" name="CFSAN QA/QC process" version="0.1.0">
|
jpayne@0
|
2 <requirements>
|
jpayne@2
|
3 <requirement type="package" version="1.69">biopython</requirement>
|
jpayne@0
|
4 <requirement type="package" version="1.4.3">matplotlib</requirement>
|
jpayne@0
|
5 <requirement type="package" version="1.9.2">numpy</requirement>
|
jpayne@0
|
6 <requirement type="package" version="0.16.1">pandas</requirement>
|
jpayne@0
|
7 </requirements>
|
jpayne@0
|
8 <command detect_errors="exit_code"><![CDATA[
|
jpayne@0
|
9 #if $reads.reads_select == 'collection'
|
jpayne@0
|
10 #set name=$reads.coll.name
|
jpayne@0
|
11 #set forward=$reads.coll.forward
|
jpayne@0
|
12 #set reverse=$reads.coll.reverse
|
jpayne@0
|
13 #else
|
jpayne@0
|
14 #set name=$forward.name
|
jpayne@0
|
15 #end if
|
jpayne@0
|
16 mkdir plots &&
|
jpayne@1
|
17 $__tool_directory__/qa_core/gims_core.py --sequencingType "Illumina MiSeq"
|
jpayne@0
|
18 --approximateGenomeSize $genomeSize
|
jpayne@0
|
19 --coverageThreshold $covThreshold
|
jpayne@0
|
20 --minimumQScore $minQ
|
jpayne@0
|
21 --output_plots plots
|
jpayne@0
|
22 --parallelism \${GALAXY_SLOTS:-6}
|
jpayne@0
|
23 --annotate $name
|
jpayne@0
|
24 --log -
|
jpayne@0
|
25 suite
|
jpayne@0
|
26 $forward
|
jpayne@0
|
27 $reverse
|
jpayne@0
|
28 -o $output
|
jpayne@0
|
29 ]]></command>
|
jpayne@0
|
30 <inputs>
|
jpayne@0
|
31 <conditional name="reads">
|
jpayne@0
|
32 <param name="reads_select" type="select" label="Paired-end collection, or two datasets from your history">
|
jpayne@0
|
33 <option value="collection">Paired collection from your history</option>
|
jpayne@0
|
34 <option value="history">Two FASTQ datasets from your history</option>
|
jpayne@0
|
35 </param>
|
jpayne@0
|
36 <when value="collection">
|
jpayne@0
|
37 <param label="Paired reads" name="coll" type="data_collection" format="fastq,fastqsanger,fastq.gz,fastqsanger.gz,fastq.bz2,fastqsanger.bz2" collection_type="paired" />
|
jpayne@0
|
38 </when>
|
jpayne@0
|
39 <when value="history">
|
jpayne@0
|
40 <param label="Forward reads" type="data" name="forward" format="fastq,fastqsanger,fastq.gz,fastqsanger.gz,fastq.bz2,fastqsanger.bz2" />
|
jpayne@0
|
41 <param label="Reverse reads" type="data" name="reverse" format="fastq,fastqsanger,fastq.gz,fastqsanger.gz,fastq.bz2,fastqsanger.bz2" />
|
jpayne@0
|
42 </when>
|
jpayne@0
|
43 </conditional>
|
jpayne@0
|
44 <param name="genomeSize" type="text" label="Approximate genome size" value="4700000" />
|
jpayne@0
|
45 <param name="covThreshold" type="text" label="Coverage threshold" value="20" />
|
jpayne@0
|
46 <param name="minQ" type="text" label="Minimum average Q score" value="27.0" />
|
jpayne@0
|
47 </inputs>
|
jpayne@0
|
48 <outputs>
|
jpayne@0
|
49 <data label="QA/QC Results" name="output" format="tabular" />
|
jpayne@0
|
50 <data label="Nucleotide plot" name="nucleo_plot" format="png" from_work_dir="nucleotide_distribution.png" />
|
jpayne@0
|
51 <data label="QScore Distribution" name="qscore_dist" format="png" from_work_dir="qscore_distribution.png" />
|
jpayne@0
|
52 <data label="QScore Histogram" name="qscore_histro" format="png" from_work_dir="qscore_histogram.png" />
|
jpayne@0
|
53 </outputs>
|
jpayne@0
|
54 <tests>
|
jpayne@0
|
55 <test>
|
jpayne@0
|
56 <param name="reads_select" value="collection" />
|
jpayne@0
|
57 <param name="coll">
|
jpayne@0
|
58 <collection type="list:paired">
|
jpayne@0
|
59 <element name="TEST TEST TEST">
|
jpayne@0
|
60 <collection type="paired">
|
jpayne@0
|
61 <element name="forward" value="Input/forward_big.fastq.gz" ftype="fastqsanger.gz" />
|
jpayne@0
|
62 <element name="reverse" value="Input/reverse_big.fastq.gz" ftype="fastqsanger.gz" />
|
jpayne@0
|
63 </collection>
|
jpayne@0
|
64 </element>
|
jpayne@0
|
65 </collection>
|
jpayne@0
|
66 </param>
|
jpayne@0
|
67 <output name="output" value="Output/output.tsv" />
|
jpayne@0
|
68 </test>
|
jpayne@0
|
69 </tests>
|
jpayne@0
|
70 <help><![CDATA[
|
jpayne@0
|
71 TODO: Fill in help.
|
jpayne@0
|
72 ]]></help>
|
jpayne@0
|
73 </tool> |