Mercurial > repos > artbio > probecoverage
diff probecoverage.xml @ 0:dbeb4a0abfc6 draft
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/probecoverage commit 102edc0c210d94e9d72f913e2d18c19220c4167c
author | artbio |
---|---|
date | Sat, 23 Sep 2017 05:11:49 -0400 |
parents | |
children | ebe5ec2e244d |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/probecoverage.xml Sat Sep 23 05:11:49 2017 -0400 @@ -0,0 +1,84 @@ +<tool id="probecoverage" name="Probe Coverage" version="0.1.0"> + <description></description> + <requirements> + <requirement type="package" version="1.4.1">samtools</requirement> + <requirement type="package" version="2.26.0">bedtools</requirement> + <requirement type="package" version="1.3.2=r3.3.1_0">r-optparse</requirement> + <requirement type="package" version="2.2.1=r3.3.1_0">r-ggplot2</requirement> + <requirement type="package" version="1.4.2=r3.3.1_0">r-reshape2</requirement> + </requirements> + <stdio> + <exit_code range="1:" level="fatal" description="Tool exception" /> + </stdio> + <command detect_errors="exit_code"><![CDATA[ + #for $file in $inputs + samtools index '$file' && + #end for + bedtools multicov + -bams + #for $file in $inputs + '$file' + #end for + -bed '$bed' > $bedtools_table && + Rscript '$__tool_directory__'/probecoverage.r + --input '$bedtools_table' + --title 'Probe coverage depth (cumulative distribution)' + --xlab 'Depth' + --ylab 'Fraction of covered regions with coverage >= Depth' + --output '$distribution' + --sample " + #for $file in $inputs + $file.element_identifier, + #end for + " + + ]]></command> + <inputs> + <param name="inputs" type="data" format="bam" label="Select multiple Bam alignments to parse" multiple="True"/> + <param name="bed" type="data" format="bed" label="Select a bed file describing the genomic regions to analyze" /> + </inputs> + + <outputs> + <data format="tabular" name="bedtools_table" label="bedtools multicov output" /> + <data format="pdf" name="distribution" label="Cumulative distribution of region coverages" /> +</outputs> + + <tests> + <test> + <param name="inputs" value="sample1,sample2,sample3" ftype="bam"/> + <param name="bed" value="probes.bed" ftype="bed"/> + <output file="coverage.tab" name="bedtools_table" /> + <output file="graph.pdf" name="distribution" /> + </test> + </tests> + + +<help> + +**What it does** + +Generates counts of alignments from multiple position-sorted BAM files +that overlap intervals in a BED file. + +Shows data as a cumulative plot of numbers of regions with coverage > x. + +**Inputs** + +bam alignment files that must be + + - sorted + - mapped to the same reference + +**Output** + +A data frame of computed counts + +A pdf file generated by R + +</help> + +<citations> + <citation type="doi">10.1093/bioinformatics/btq033</citation> +</citations> +</tool> +