Mercurial > repos > bgruening > hicup2juicer
diff hicup2juicer.xml @ 0:396bdc215399 draft
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/hicup commit 50173a1994a72d9774fd46777de94dd02d35bd42
author | bgruening |
---|---|
date | Mon, 24 Oct 2022 15:58:52 +0000 |
parents | |
children | 82641e9a86cb |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/hicup2juicer.xml Mon Oct 24 15:58:52 2022 +0000 @@ -0,0 +1,83 @@ +<tool id="hicup2juicer" name="Hicup to juicer converter" version="@VERSION@+galaxy0"> + <description></description> + <macros> + <import>hicup_macros.xml</import> + </macros> + <expand macro="requirements_hicup" /> + <command detect_errors="exit_code"><![CDATA[ + #if $input_file.ext != 'sam': + #set ext='bam' + #else: + #set ext='sam' + #end if + ln -s '$input_file' input.$ext && + hicup2juicer input.$ext + ]]> + </command> + <inputs> + <param name="input_file" type="data" format="qname_sorted.bam,sam" label="Output of HiCUP to convert" /> + </inputs> + <outputs> + <data name="output" format="tabular" label="${tool.name} on ${on_string}: Pairs in juicebox format" from_work_dir="*.prejuicer" /> + </outputs> + <tests> + <test expect_num_outputs="1"> + <param name="input_file" value="dataset1_2.hicup.bam"/> + <output name="output" value="dataset1_2.hicup.bam.prejuicer"/> + </test> + </tests> + <help><![CDATA[ + +HiCUP homepage: www.bioinformatics.babraham.ac.uk/projects/hicup + +The hicup2juicer script converts HiCUP BAM/SAM files to a format compatible +with Juicer and JuiceBox( https://github.com/aidenlab/juicer ). +Outputfiles generated by this script may be converted to Juicer ".hic" files +using the "pre" command as described at: https://github.com/aidenlab/juicer/wiki/Pre + +The script does not use restriction site coordinates when generating output. + +FUNCTION + +HiCUP generates SAM/BAM files of mapped, filtered paired-end reads +constituting the sequenced valid Hi-C di-tags. These may then be analysed by a +variety of specialised tools, but before this is possible the datasets will +need parsing into an appropriate format. + +The hicup2juicer script converts HiCUP BAM/SAM files to a tab-delimited format +comprising 7 columns, with read pairs on the same line: + +<readname> <str1> <chr1> <pos1> <frag1> <str2> <chr2> <pos2> <frag2> <mapq1> <mapq2> +str = strand (0 for forward, anything else for reverse) +chr = chromosome (must be a chromosome in the genome) +pos = position +frag = restriction site fragment +mapq = mapping quality score + +Column1: Readpair index number (starting at 1) +Column2: forward read strand (0 = positive strand, 1 = negative strand) +Column3: forward read chromosome name +Column4: forward read position +Column5: forward read fragment id (set to the dummy value 0) +Column6: reverse read strand (0 = positive strand, 1 = negative strand) +Column7: reverse read chromosome name +Column8: reverse read position +Column9: reverse read fragment id (set to the dummy value 1) +Column10: forward read MAPQ score +Column11: reverse read MAPQ score + + +COMMAND LINE OPTIONS + +--help Print help message and exit +--version Print the program version and exit +--zip Write output to a gzip file + +Full instructions on running the pipeline can be found at: +www.bioinformatics.babraham.ac.uk/projects/hicup + +Steven Wingett, Babraham Institute, Cambridge, UK +]]> + </help> + <expand macro="citation_hicup" /> +</tool> \ No newline at end of file