<tool id="cshl_fastx_trimmer" version="1.0.0" name="Trim sequences">
	<description></description>
    <requirements>
        <requirement type="package" version="0.0.13">fastx_toolkit</requirement>
    </requirements>
	<command>zcat -f '$input' | fastx_trimmer -v -f $first -l $last -o $output
#if $input.ext == "fastqsanger":
-Q 33
#end if
	</command>

	<inputs>
		<param format="fasta,fastqsolexa,fastqsanger" version="1.0.0" name="input" type="data" label="Library to clip" />

		<param version="1.0.0" name="first" size="4" type="integer" value="1">
			<label>First base to keep</label>
		</param>

		<param version="1.0.0" name="last" size="4" type="integer" value="21">
			<label>Last base to keep</label>
		</param>
	</inputs>

	<tests>
		<test>
			<!-- Trim a FASTA file - remove first four bases (e.g. a barcode) -->
			<param version="1.0.0" name="input" value="fastx_trimmer1.fasta" />
			<param version="1.0.0" name="first" value="5"/>
			<param version="1.0.0" name="last" value="36"/>
			<output version="1.0.0" name="output" file="fastx_trimmer1.out" />
		</test>
		<test>
			<!-- Trim a FASTQ file - remove last 9 bases (e.g. keep only miRNA length sequences) -->
			<param version="1.0.0" name="input" value="fastx_trimmer2.fastq" ftype="fastqsolexa"/>
			<param version="1.0.0" name="first" value="1"/>
			<param version="1.0.0" name="last" value="27"/>
			<output version="1.0.0" name="output" file="fastx_trimmer2.out" />
		</test>
	</tests>

	<outputs>
		<data format="input" version="1.0.0" name="output" metadata_source="input" />
	</outputs>
	<help>
**What it does**

This tool trims (cut bases from) sequences in a FASTA/Q file.
  
--------

**Example**

Input Fasta file (with 36 bases in each sequences)::

    >1-1
    TATGGTCAGAAACCATATGCAGAGCCTGTAGGCACC
    >2-1
    CAGCGAGGCTTTAATGCCATTTGGCTGTAGGCACCA
    

Trimming with First=1 and Last=21, we get a FASTA file with 21 bases in each sequences (starting from the first base)::

    >1-1
    TATGGTCAGAAACCATATGCA
    >2-1
    CAGCGAGGCTTTAATGCCATT

Trimming with First=6 and Last=10, will generate a FASTA file with 5 bases (bases 6,7,8,9,10) in each sequences::

    >1-1
    TCAGA
    >2-1
    AGGCT
    
    ------

This tool is based on `FASTX-toolkit`__ by Assaf Gordon.

 .. __: http://hannonlab.cshl.edu/fastx_toolkit/
    
</help>
<!-- FASTX-Trimmer is part of the FASTX-toolkit, by A.Gordon (gordon@cshl.edu) -->
</tool>