Mercurial > repos > iuc > multiqc
changeset 23:abfd8a6544d7 draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/multiqc commit 8a91236cee4d408ae2b53a3e9b6daebc332d631a
author | iuc |
---|---|
date | Sat, 10 Dec 2022 11:22:18 +0000 |
parents | 9a913cdee30e |
children | f7e2f1eb3a16 |
files | macros.xml multiqc.xml test-data/junction_annotation.txt test-data/junction_saturation.txt |
diffstat | 4 files changed, 84 insertions(+), 22 deletions(-) [+] |
line wrap: on
line diff
--- a/macros.xml Wed Nov 03 15:17:39 2021 +0000 +++ b/macros.xml Sat Dec 10 11:22:18 2022 +0000 @@ -1,6 +1,6 @@ <macros> <token name="@TOOL_VERSION@">1.11</token> - <token name="@VERSION_SUFFIX@">0</token> + <token name="@VERSION_SUFFIX@">1</token> <xml name="bio_tools"> <xrefs> <xref type="bio.tools">multiqc</xref> @@ -218,9 +218,11 @@ #set file_path = os.path.join($repeat_dir, 'file_' + str($k) + '.GC.xls') ln -s '$file' '$file_path' && #end for - #elif str($repeat2.type.type) == "junction_annotation" - #set $pattern = "Group Total_bases Tag_count Tags/Kb" - @LN_3_FILES@ + #elif str($repeat2.type.type) == "junction_saturation" + #for $k, $file in enumerate($repeat2.type.input) + #set file_path = os.path.join($repeat_dir, 'file_' + str($k) + '.junctionSaturation_plot.r') + ln -s '${file}' '$file_path' && + #end for #elif str($repeat2.type.type) == "read_distribution" #set $pattern = "Group Total_bases Tag_count Tags/Kb" @LN_3_FILES@
--- a/multiqc.xml Wed Nov 03 15:17:39 2021 +0000 +++ b/multiqc.xml Sat Dec 10 11:22:18 2022 +0000 @@ -413,42 +413,42 @@ <repeat name="output" title="RSeQC output" min="1"> <conditional name="type"> <param name="type" type="select" label="Type of RSeQC output?"> - <option value="bam_stat">bam_stat</option> - <option value="gene_body_coverage">gene_body_coverage</option> - <option value="infer_experiment">infer_experiment</option> - <option value="inner_distance">inner_distance</option> - <option value="junction_annotation">junction_annotation</option> - <option value="junction_saturation">junction_saturation</option> - <option value="read_distribution">read_distribution</option> - <option value="read_duplication_pos">read_duplication_pos</option> - <option value="read_gc">read_gc</option> + <option value="bam_stat">BAM/SAM mapping stats</option> + <option value="gene_body_coverage">Gene body coverage</option> + <option value="infer_experiment">Infer experiment</option> + <option value="inner_distance">Inner distance</option> + <option value="junction_annotation">Junction annotation</option> + <option value="junction_saturation">Junction saturation</option> + <option value="read_distribution">Read distribution</option> + <option value="read_duplication_pos">Read duplication</option> + <option value="read_gc">Read GC</option> </param> <when value="bam_stat"> - <param name="input" type="data" format="txt,tabular" multiple="true" label="RSeQC bam_stat output" help="It should contain 'Proper-paired reads map to different chrom:'"/> + <param name="input" type="data" format="txt,tabular" multiple="true" label="RSeQC BAM stat: output file" help="It should contain 'Proper-paired reads map to different chrom:'"/> </when> <when value="gene_body_coverage"> - <param name="input" type="data" format="txt" multiple="true" label="RSeQC gene_body_coverage output"/> + <param name="input" type="data" format="txt" multiple="true" label="RSeQC gene body coverage: stats file"/> </when> <when value="infer_experiment"> - <param name="input" type="data" format="txt" multiple="true" label="RSeQC infer_experiment output" help="It should contain 'Fraction of reads explained by'"/> + <param name="input" type="data" format="txt" multiple="true" label="RSeQC infer experiment: configuration output" help="It should contain 'Fraction of reads explained by'"/> </when> <when value="inner_distance"> - <param name="input" type="data" format="txt,xls,tabular" multiple="true" label="RSeQC inner_distance output"/> + <param name="input" type="data" format="txt" multiple="true" label="RSeQC inner distance: frequency file"/> </when> <when value="junction_annotation"> - <param name="input" type="data" format="txt,xls,tabular" multiple="true" label="RSeQC junction_annotation output" help="It should contain 'Partial Novel Splicing Junctions:'"/> + <param name="input" type="data" format="txt" multiple="true" label="RSeQC junction annotation: stats file" help="It should contain 'Partial Novel Splicing Junctions:'"/> </when> <when value="junction_saturation"> - <param name="input" type="data" format="txt" multiple="true" label="RSeQC junction_saturation output"/> + <param name="input" type="data" format="txt" multiple="true" label="RSeQC junction saturation: junction saturation plot Rscript file"/> </when> <when value="read_distribution"> - <param name="input" type="data" format="txt,xls,tabular" multiple="true" label="RSeQC read_distribution output" help="It should contain 'Group Total_bases Tag_count Tags/Kb'"/> + <param name="input" type="data" format="txt" multiple="true" label="RSeQC read distribution: stats output" help="It should contain 'Group Total_bases Tag_count Tags/Kb'"/> </when> <when value="read_duplication_pos"> - <param name="input" type="data" format="csv,xls,tabular" multiple="true" label="RSeQC read_duplication_pos output"/> + <param name="input" type="data" format="xls" multiple="true" label="RSeQC read duplication: positions XLS file"/> </when> <when value="read_gc"> - <param name="input" type="data" format="csv,xls,tabular" multiple="true" label="RSeQC read_gc output"/> + <param name="input" type="data" format="csv,xls,tabular" multiple="true" label="RSeQC read GC output"/> </when> </conditional> </repeat> @@ -1076,6 +1076,36 @@ </output> <output_collection name="stats" type="list" count="2"/> </test> + <!--Test 07--> + <test expect_num_outputs="2"> + <repeat name="results"> + <conditional name="software_cond"> + <param name="software" value="rseqc" /> + <repeat name="output"> + <conditional name="type"> + <param name="type" value="junction_annotation"/> + <param name="input" value="junction_annotation.txt"/> + </conditional> + </repeat> + <repeat name="output"> + <conditional name="type"> + <param name="type" value="junction_saturation"/> + <param name="input" value="junction_saturation.txt"/> + </conditional> + </repeat> + </conditional> + </repeat> + <param name="title" value="RSEQC report"/> + <!-- <param name="flat" value="true"/> --> + <output name="html_report"> + <assert_contents> + <has_text text="RSEQC report" /> + <has_text text="Junction Annotation" /> + <has_text text="Junction Saturation" /> + </assert_contents> + </output> + <output_collection name="stats" type="list" count="2"/> + </test> </tests> <help><![CDATA[ **What it does**
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/junction_annotation.txt Sat Dec 10 11:22:18 2022 +0000 @@ -0,0 +1,18 @@ +Reading reference bed file: /tmp/tmpwdbqfb2i/files/000/dataset_3.dat ... Done +Load BAM file ... Done + +=================================================================== +Total splicing Events: 4 +Known Splicing Events: 1 +Partial Novel Splicing Events: 1 +Novel Splicing Events: 1 +Filtered Splicing Events: 1 + +Total splicing Junctions: 3 +Known Splicing Junctions: 1 +Partial Novel Splicing Junctions: 1 +Novel Splicing Junctions: 1 + +=================================================================== +Create BED file ... +Create Interact file ...
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/junction_saturation.txt Sat Dec 10 11:22:18 2022 +0000 @@ -0,0 +1,12 @@ +pdf('output.junctionSaturation_plot.pdf') +x=c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95,100) +y=c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1) +z=c(0,0,0,0,0,0,1,1,1,1,1,1,1,2,2,2,2,2,2,3) +w=c(0,0,0,0,0,0,1,1,1,1,1,1,1,2,2,2,2,2,2,2) +m=max(0,0,0) +n=min(0,0,0) +plot(x,z/1000,xlab='percent of total reads',ylab='Number of splicing junctions (x1000)',type='o',col='blue',ylim=c(n,m)) +points(x,y/1000,type='o',col='red') +points(x,w/1000,type='o',col='green') +legend(5,0, legend=c("All junctions","known junctions", "novel junctions"),col=c("blue","red","green"),lwd=1,pch=1) +dev.off()