|
8
|
1 <tool id="factera_1_4_3" name="FACTERA" version="1.4.3">
|
|
|
2
|
|
|
3 <description>Fusion And Chromosomal Translocation Enumeration and Recovery Algorithm</description>
|
|
|
4
|
|
|
5 <requirements>
|
|
|
6 <requirement type="package" version="0.1.19">package_samtools_0_1_19</requirement>
|
|
|
7 <requirement type="package" version="2.2.26">package_blast_plus_2_2_26</requirement>
|
|
|
8 <requirement type="package" version="1.0.0">package_ucsc_user_apps</requirement>
|
|
|
9 <requirement type="package" version="5.18.1">perl</requirement>
|
|
|
10 <requirement type="package" version="1.4.3">factera_1_4_3</requirement>
|
|
|
11 <requirement type="set_environment">FACTERA_DIR</requirement>
|
|
|
12 </requirements>
|
|
|
13
|
|
|
14 <command>
|
|
|
15 <!-- Run commands using Galaxy perl -->
|
|
|
16 #if $cond.additionalOpt == "no"
|
|
|
17 perl \$FACTERA_DIR/factera-v1.4.3.pl -o ./ $sortedBam $exonBed $genome2Bit
|
|
|
18 #else
|
|
|
19 perl \$FACTERA_DIR/factera-v1.4.3.pl -r $cond.r -m $cond.m -x $cond.x -s $cond.s -f $cond.f -S $cond.S
|
|
|
20 -k $cond.k -c $cond.c -b $cond.b -p $cond.p -a $cond.a -o ./ $sortedBam $exonBed $genome2Bit
|
|
|
21 #end if
|
|
|
22
|
|
|
23 <!-- Keep only the fusion file -->
|
|
|
24 cat *fusions.txt > $outfile;
|
|
|
25 </command>
|
|
|
26
|
|
|
27 <inputs>
|
|
|
28 <param format="bam" name="sortedBam" type="data" label="Alignment in BAM format" />
|
|
|
29 <param format="bed" name="exonBed" type="data" label="Genomic coordinates with gene/exon names in fourth column" />
|
|
|
30 <param format="twobit" name="genome2Bit" type="data" label="Two bit reference sequence" />
|
|
|
31 <conditional name="cond">
|
|
|
32 <param type="select" name="additionalOpt" label="Additional Options">
|
|
|
33 <option value="yes">Yes</option>
|
|
|
34 <option value="no" selected="True">No</option>
|
|
|
35 </param>
|
|
|
36 <when value="yes">
|
|
|
37 <!-- Names follow command line flags in factera.pl -->
|
|
|
38 <param type="integer" name="r" label="Minimum number of breakpoint-scanning reads required for output" value="5"/>
|
|
|
39 <param type="integer" name="m" label="Minimum number of discordant reads required for each candidate fusion" value="2"/>
|
|
|
40 <param type="integer" name="x" label="Maximum number of breakpoints to examine for each pair of genomic regions" value="5"/>
|
|
|
41 <param type="integer" name="s" label="Minimum number of reads with the same breakpoint" value="1"/>
|
|
|
42 <param type="float" name="f" min="0" max="1" label="Minimum fraction of read bases required for alignment to fusion template" value=".9"/>
|
|
|
43 <param type="float" name="S" min="0" max="1" label="Minimum similarity required for read to match fusion template" value=".95"/>
|
|
|
44 <param type="integer" name="k" label="K-mer size (number of bases) for fragment comparison" value="10"/>
|
|
|
45 <param type="integer" name="c" label="Minimum size of soft-clipped region to consider" value="16"/>
|
|
|
46 <param type="integer" name="b" label="Number of bases flanking breakpoint for fusion template" value="500"/>
|
|
|
47 <param type="integer" name="p" label="Number of threads for blastn search" value="4"/>
|
|
|
48 <param type="integer" name="a" label="Number of bases flanking breakpoint to provide in output" value="50"/>
|
|
|
49 </when>
|
|
|
50 </conditional>
|
|
|
51 </inputs>
|
|
|
52
|
|
|
53 <outputs>
|
|
|
54 <data format="tabular" name="outfile" />
|
|
|
55 </outputs>
|
|
|
56
|
|
|
57 <help>This tool does Factera</help>
|
|
|
58
|
|
|
59 </tool>
|
