Mercurial > repos > rnateam > atactk_trim_adapters
diff trim_adapters.xml @ 0:9839a3fe72f7 draft default tip
planemo upload for repository https://github.com/bgruening/galaxytools/blob/master/tools/trim_adapters commit 3d0b670cda6522e5c442b144785b2f9f517f103d
author | rnateam |
---|---|
date | Wed, 20 Jun 2018 15:48:09 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/trim_adapters.xml Wed Jun 20 15:48:09 2018 -0400 @@ -0,0 +1,83 @@ +<tool id="atactk_trim_adapters" name="Trim putative adapter sequence" version="0.1.6"> + <requirements> + <requirement type="package" version="0.1.6">atactk</requirement> + </requirements> + <command detect_errors="exit_code"><![CDATA[ + ln -s '$input1' ./forward.${input1.ext} && + ln -s '$input2' ./reverse.${input2.ext} && + trim_adapters + #if '$settings.advanced' == 'advanced' + -d '$settings.edit_distance' + -f '$settings.fudge' + -s '$settings.trim_start' + -r '$settings.rc_length' + #end if + ./forward.${input1.ext} ./reverse.${input2.ext} + ]]></command> + <inputs> + <param type="data" name="input1" format="fastq,fastq.gz" /> + <param type="data" name="input2" format="fastq,fastq.gz" /> + <conditional name="settings"> + <param name="advanced" type="select" label="Specify advanced parameters"> + <option value="simple" selected="true">No, use program defaults.</option> + <option value="advanced">Yes, see full parameter list.</option> + </param> + <when value="simple"></when> + <when value="advanced"> + <param name="edit_distance" label="The maximum edit distance permitted when aligning the paired reads" type="integer" min="1" value="1" help="(-d)" /> + <param name="fudge" label="An arbitrary number of extra bases to trim from the ends of reads" type="integer" min="1" value="1" help="(-f)" /> + <param name="trim_start" label="Trim this number of bases from the start of each sequence" type="integer" min="0" value="0" help="(-s)"/> + <param name="rc_length" label="Use the reverse complement of this number of from the beginning of the reverse read to align reads" type="integer" min="1" value="20" help="(-r)"/> + </when> + </conditional> + </inputs> + <outputs> + <data name="output1" format="fastq.gz" from_work_dir="forward.trimmed.fastq.gz" /> + <data name="output2" format="fastq.gz" from_work_dir="reverse.trimmed.fastq.gz" /> + </outputs> + <tests> + <test> + <param name="input1" value="SP1_f.fastq" ftype="fastq" /> + <param name="input2" value="SP1_r.fastq" ftype="fastq" /> + <output name="output1" file="SP1_f.trimmed.fastq.gz" decompress="True"/> + <output name="output2" file="SP1_r.trimmed.fastq.gz" decompress="True"/> + </test> + <test> + <param name="input1" value="SP1_f.fastq.gz" ftype="fastq.gz" /> + <param name="input2" value="SP1_r.fastq.gz" ftype="fastq.gz" /> + <output name="output1" file="SP1_f.trimmed.fastq.gz" decompress="True"/> + <output name="output2" file="SP1_r.trimmed.fastq.gz" decompress="True"/> + </test> + </tests> + <help> +<![CDATA[ + +**What it does** + +The trim_adapters utility is based on a script by Jason Buenrostro. +Instead of looking for known adapter sequence, it aligns paired reads to each other +and trims off sequence outside the alignment. More precisely, it searches +the forward read for the reverse complement of a specified number of bases +(20 by default) at the beginning of the reverse read, then falls back to finding +the best alignment of the two reads, using the minimum Levenshtein distance between them. + +**Input** + +It requires 2 inputs: The (optionally gzipped) FASTQ file containing the +forward reads and the (optionally gzipped) FASTQ file containing the +reverse reads. + +**Output** + +Generates 2 gzipped fastq files with adapters trimmed. + + ]]></help> + <citations> + <citation type="bibtex">@unpublished{atactk: a toolkit for ATAC-seq data, + title = "atactk: a toolkit for ATAC-seq data", + author = "The Parker Lab at the University of Michigan", + url = "https://github.com/ParkerLab/atactk/", + year = "2018" +}</citation> + </citations> +</tool>