comparison CoverageReport.xml @ 0:49788db492f7 draft default tip

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/coverage_report commit 1b647bb088f62c1369d963f030caecaa9ee003c7

0:49788db492f7

<tool id="CoverageReport2" name="Panel Coverage Report" version="0.0.4">
    <description>Provide mapping statistics</description>
    <requirements>
      <requirement type="package" version="1.75">perl-number-format</requirement>    
      <requirement type="package" version="3.4.1">R</requirement>
      <requirement type="package" version="2.17.0">bedtools</requirement>
      <requirement type="package" version="0.1.18">samtools</requirement>
      <requirement type="package" version="0.1.6">tectonic</requirement>      
    </requirements>
    <command detect_errors="exit_code"><![CDATA[      
    perl '$__tool_directory__/CoverageReport.pl'
    ## input files
    -b '$input1'
    -t '$input2'
  
    ## output files
    -o '$output1'
  
    ## run parameters
    $perGene
    $PositionLevel
    -m $threshold
    -f $frac 
    ## sample name
    -n '$input1.name'
    ]]></command>
    <inputs>
      <param name="input1" type="data" format="bam" label="BAM file" help="BAM file of mapped reads" />
      <param name="input2" type="data" format="bed" label="Target Regions BED" help="BED file containing regions of interest. See below for format" />
      <param name="threshold" type="integer" value="40" label="Minimal Coverage Threshold" help="Default: 40" />   
      <param name="frac" type="float" value="0.2" label="Fraction of Average Coverage for usage in plot" help="Default: 0.2" />
      <param name="perGene" type="boolean" truevalue="-r" falsevalue="" checked="True" label="Plot exon coverages for all genes in targets"/>
      <param name="PositionLevel" type="select" label="Perform Per Exon Analysis" help="Only Failed: Only those exons not reaching global coverage above threshold, or 100%. All Exons: This can take a very long time for large panels! Select all failed to check all exons for local failures." >
	<option value='' selected="true">None</option>
	<option value='-s'>Plot Only Globally Failed</option>
	<option value='-S'>Plot All Failed Exons</option>
	<option value='-A'>Plot All Exons</option>
	<option value='-L'>List All Failed Exons</option>
      </param>
    </inputs>
    <outputs>
      <data format="pdf" name="output1" label="${tool.name} on ${on_string}: PDF Report"/>
    </outputs>
    <tests>
      <test>
	<param name="input1" value="ecoli.bam" />
	<param name="input2" value="ecoli_6.bed" />
	<output name="output1" value="answer.pdf" compare="sim_size"/>
      </test>
    </tests>
    <help><![CDATA[
**What it does**

This tool creates a coverage report for QC purposes. By default, average coverage statistics are provided, taken from samtools flagstats. If specified, it can also create overviews per gene in the BED file, and sub-exon plots for failed exons. 

------

**BED format**

The BED file containing targets of interest has very specific format requirements. You **must** use the following format::

  Column 1: Chromosome : Use the same syntax as the references used by Galaxy. Check your sam-headers for the correct format. ('chr1' vs '1')
  Column 2: Start Position
  Column 3: End Position
  Column 4: Target Name. Use : "GENE-NAME&lt;space&gt;Exon_number" : This is split on the space after 'GeneName' for correct grouping.
  Column 5: Score : ignored, use '0'
  Column 6: Strand: ignored,'+' or '-'

.. class:: infomark 

Note: The exons for the plots will be ordered in the same way as the exons in the BED file. 

------

**Input formats**

BAM file for reads, BED file for targets.

------

**Outputs**

The output files are a PDF report. 
The output tables are (tab seperated txt files): 

**Targets.Global.Coverage** : Original BED file + following columns::
  - Total coverage in target
  - Bases in target with coverage
  - Length of target
  - Percent of target covered

**Targets.Position.Coverage** : Original BED file + following columns::
  - Position in target region
  - Coverage at position

    ]]></help>  
    <citations>
      <citation type="bibtex">@ARTICLE{a1,
      title = {This tool creates a coverage report for QC purposes.},
      author = {Philip Mabon},
      url = {https://github.com/galaxyproject/tools-iuc/}
      }
      }</citation>
    </citations>
</tool>