Mercurial > repos > iuc > bam_to_scidx
changeset 0:741b8af6386c draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bam_to_scidx commit dfaab7cbfdee82c9fe4ff34ce02b42fc456b9db9
author | iuc |
---|---|
date | Sat, 12 Dec 2015 15:46:03 -0500 |
parents | |
children | ebbf2a6a4c55 |
files | BAMtoscIDX.jar bam_to_scidx.xml test-data/input.bam test-data/output.scidx |
diffstat | 4 files changed, 84 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/bam_to_scidx.xml Sat Dec 12 15:46:03 2015 -0500 @@ -0,0 +1,63 @@ +<tool id="bam_to_scidx" name="Convert BAM to ScIdx" version="1.0.0"> + <description></description> + <command> + <![CDATA[ + ln -f -s "${input_bam.metadata.bam_index}" "${input_bam}.bai" && + java -jar $__tool_directory__/BAMtoscIDX.jar -b "$input_bam" -i "${input_bam}.bai" -p $require_proper_mate_pairing -r $read -m $min_insert_size -M $max_insert_size -o "$output" 1>/dev/null + ]]> + </command> + <inputs> + <param name="input_bam" type="data" format="bam" label="BAM file" /> + <param name="require_proper_mate_pairing" type="select" label="Require proper mate-pairing?" help="Required if filtering by insert size for single-end Reqd 1 (below)." > + <option value="1" selected="True">Yes</option> + <option value="0">No</option> + </param> + <param name="read" type="select" label="Read to output"> + <option value="0" selected="True">Read1</option> + <option value="1">Read2</option> + <option value="2">Combined</option> + </param> + <param name="min_insert_size" type="integer" value="" optional="True" min="0" label="Minimum insert size to output" help="Will not filter out single-end Read 1 unless proper mate-pairing is required (above)."/> + <param name="max_insert_size" type="integer" value="" optional="True" min="0" label="Maximum insert size to output" help="Will not filter out single-end Read 1 unless proper mate-pairing is required (above)." /> + </inputs> + <outputs> + <data name="output" format="scidx" /> + </outputs> + <tests> + <test> + <param name="input_bam" value="input.bam" ftype="bam" /> + <param name="require_proper_mate_pairing" value="1" /> + <param name="read" value="0" /> + <output name="output" file="output.scidx" lines_diff="1" ftype="scidx" /> + </test> + </tests> + <help> + +**What it does** + +Converts BAM data to ScIdx, the Strand-specific coordinate count format, which is used by tools within +the Chip-exo Galaxy flavor. ScIdx files are 1-based. The format consists of 5 columns: the chromosome, +the position of the genomic coordinate, the number of tags on the forward strand, the number of tags on +the reverse strand and the number of total tags on the position. With pair-end reads, only the 5’ end of +READ1 will be used to create the ScIdx data file. Tools that use this format include GeneTrack and MultiGPS. + +----- + +**Options** + +* **Require proper mate-pairing?** - Select **Yes** to require proper mate paring. Filtering by insert size parameters will not filter out single-end Read 1 unless proper mate-pairing is required. +* **Minimum insert size to output** - Insert size below the minimum will be filtered from the results, but single-end Read 1 will not be filtered unless proper mate-pairing is required. +* **Maximum insert size to output** - Insert size above the maximum will be filtered from the results, but single-end Read 1 will not be filtered unless proper mate-pairing is required. + + </help> + <citations> + <citation type="bibtex"> + @unpublished{None, + author = {Lai, William}, + title = {None}, + year = {None}, + eprint = {None}, + url = {http://www.huck.psu.edu/content/research/independent-centers-excellence/center-for-eukaryotic-gene-regulation} + }</citation> + </citations> +</tool>
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output.scidx Sat Dec 12 15:46:03 2015 -0500 @@ -0,0 +1,21 @@ +#2015-11-23 20:18:56.51;input.bam;READ1 +chrom index forward reverse value +chr17 1820100 0 1 1 +chr17 1849175 1 0 1 +chr17 2778159 0 1 1 +chr17 3205518 0 1 1 +chr17 12974350 1 0 1 +chr17 20624912 0 1 1 +chr17 27071176 0 1 1 +chr17 34920819 1 0 1 +chr17 39810583 0 1 1 +chr17 43542379 0 1 1 +chr17 45800656 1 0 1 +chr17 47124830 0 1 1 +chr17 48697995 0 1 1 +chr17 54474619 0 1 1 +chr17 64901353 0 1 1 +chr17 69072332 1 0 1 +chr17 76184543 1 0 1 +chr17 76518878 0 1 1 +chr17 76941015 1 0 1