Mercurial > repos > jjohnson > shear
annotate shear_sv.xml @ 4:a82400332451
Use shear_wrapper.py to generate reference indexes when needed
| author | Jim Johnson <jj@umn.edu> |
|---|---|
| date | Tue, 16 Jul 2013 12:38:48 -0500 |
| parents | 630532975fa9 |
| children | aaaa5a071ff0 |
| rev | line source |
|---|---|
| 0 | 1 <?xml version="1.0"?> |
| 2 <tool id="shear_sv" name="SHEAR-SV" version="0.0.1"> | |
| 3 <description>find SVs and estimate the heterogeneity levels from a BAM</description> | |
| 4 <requirements> | |
| 5 <requirement type="package" version="0.1.2">shear</requirement> | |
| 6 <requirement type="package" version="0.1.18">samtools</requirement> | |
| 7 <requirement type="package" version="0.6.2">bwa</requirement> | |
| 8 <requirement type="package" version="1.0.1">crest</requirement> | |
| 9 <requirement type="package" version="latest">kent</requirement> | |
| 10 </requirements> | |
| 11 <!-- | |
| 12 <version_command></version_command> | |
| 13 --> | |
|
4
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
14 <command interpreter="python"> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
15 shear_wrapper.py -j \$SHEAR_JAR_PATH/SHEAR.jar --command sv |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
16 -p $prefix |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
17 -b $bamfile |
| 0 | 18 #if $genomeSource.refGenomeSource == 'indexed': |
| 19 -f $genomeSource.ref_fastas.fields.path | |
| 20 -i $genomeSource.bwa_indices.fields.path | |
| 3 | 21 -t $genomeSource.twobit_indices.fields.path |
|
4
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
22 #elif $genomeSource.refGenomeSource == 'svidx': |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
23 -f $genomeSource.svidx.extra_files_path/${genomeSource.svidx.metadata.base_name}.fa |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
24 -i $genomeSource.svidx.extra_files_path/${genomeSource.svidx.metadata.base_name}.fa |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
25 -t $genomeSource.svidx.extra_files_path/${genomeSource.svidx.metadata.base_name}.2bit |
| 0 | 26 #else: |
| 27 -f $genomeSource.ref_fasta | |
|
4
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
28 #if $genomeSource.save_svidx: |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
29 -D $sv_idx.extra_files_path |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
30 #end if |
| 0 | 31 #end if |
|
4
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
32 --report $report |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
33 --sdi $sdi |
| 0 | 34 </command> |
| 35 <inputs> | |
| 36 <param name="bamfile" type="data" format="bam" label="BAM alignment file containing the input sequences to the assembly."/> | |
|
4
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
37 <param name="prefix" type="text" value="shear_sv" label="Prefix for all generated files"> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
38 <validator type="regex" message="Prefix should start with a letter and contain only letter, digit, and '_' or '-' characters">[a-zA-Z0-9][_a-zA-Z0-9-]*</validator> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
39 </param> |
| 0 | 40 <!-- reference data --> |
| 41 <conditional name="genomeSource"> | |
| 42 <param name="refGenomeSource" type="select" label="Will you select a reference genome from your history or use a cached file?"> | |
| 43 <option value="indexed" selected="true">Use a cached reference genome</option> | |
|
4
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
44 <option value="history">Use fasta from the history</option> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
45 <option value="svidx">Use SHEAR sv index from the history</option> |
| 0 | 46 </param> |
| 47 <when value="indexed"> | |
|
1
a718ecb27e02
Comment out getting reference from history, need wrapper for that
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
48 <param name="ref_fastas" type="select" label="Select a reference genome fasta"> |
| 0 | 49 <options from_data_table="all_fasta"> |
| 50 <filter type="sort_by" column="2" /> | |
| 51 <validator type="no_options" message="No reference fasta files are available" /> | |
| 52 </options> | |
| 53 </param> | |
|
1
a718ecb27e02
Comment out getting reference from history, need wrapper for that
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
54 <param name="twobit_indices" type="select" label="Select the reference genome 2bit"> |
| 0 | 55 <options from_data_table="twobit"> |
| 56 <filter type="sort_by" column="0" /> | |
|
1
a718ecb27e02
Comment out getting reference from history, need wrapper for that
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
57 <!-- |
| 0 | 58 <filter type="data_meta" ref="ref_fastas" key="dbkey" column="0" /> |
|
1
a718ecb27e02
Comment out getting reference from history, need wrapper for that
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
59 --> |
| 0 | 60 <validator type="no_options" message="No indexes are available" /> |
| 61 </options> | |
| 62 </param> | |
|
1
a718ecb27e02
Comment out getting reference from history, need wrapper for that
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
63 <param name="bwa_indices" type="select" label="Select the reference genome bwa indices"> |
| 0 | 64 <options from_data_table="bwa_indexes"> |
| 65 <filter type="sort_by" column="2" /> | |
|
1
a718ecb27e02
Comment out getting reference from history, need wrapper for that
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
66 <!-- |
| 0 | 67 <filter type="data_meta" ref="ref_fastas" key="dbkey" column="0" /> |
|
1
a718ecb27e02
Comment out getting reference from history, need wrapper for that
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
68 --> |
| 0 | 69 <validator type="no_options" message="No indexes are available" /> |
| 70 </options> | |
| 71 </param> | |
| 72 </when> | |
| 73 <when value="history"> | |
|
4
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
74 <param name="ref_fasta" type="data" format="fasta" label="Select a reference from history" /> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
75 <param name="save_svidx" type="boolean" truevalue="yes" falsevalue="no" label="Save sv generated reference indexes"/> |
| 0 | 76 </when> |
|
4
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
77 <when value="svidx"> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
78 <param name="svidx" type="data" format="shear.svidx" label="Select a reference from history" /> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
79 </when> |
| 0 | 80 </conditional> |
| 81 </inputs> | |
| 82 <stdio> | |
| 83 <exit_code range="1:" level="fatal" description="Error" /> | |
| 84 </stdio> | |
| 85 <outputs> | |
|
4
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
86 <data format="shear.sdi" name="sdi" label="${tool.name} on ${on_string}: ${prefix}.sdi" /> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
87 <data format="txt" name="report" label="${tool.name} on ${on_string}: ${prefix}.report" /> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
88 <data format="shear.svidx" name="sv_idx" label="${tool.name} on ${on_string}: ${prefix}.svidx"> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
89 <filter>genomeSource['refGenomeSource'] == 'history' and genomeSource['save_svidx'] == True</filter> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
90 </data> |
| 0 | 91 </outputs> |
| 92 <tests> | |
| 93 <test> | |
|
4
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
94 <param name="prefix" value="shear_test" /> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
95 <param name="bamfile" value="simulated-data.bam" ftype="bam" /> |
| 0 | 96 <param name="refGenomeSource" value="history"/> |
|
4
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
97 <param name="ref_fasta" value="syn.fa" ftype="fasta" /> |
|
a82400332451
Use shear_wrapper.py to generate reference indexes when needed
Jim Johnson <jj@umn.edu>
parents:
3
diff
changeset
|
98 <output name="sdi" file="shear_sv.sdi"/> |
| 0 | 99 <output name="report" file="shear_sv.report"/> |
| 100 </test> | |
| 101 </tests> | |
| 102 <help> | |
| 103 SHEAR-SV will process a BAM alignment file to find SVs and estimate the | |
| 104 heterogeneity level of those SVs. It will output a *.report file reporing the | |
| 105 SVs and their heterogeneity levels, and it will output a *.sdi file that can be | |
| 106 used for SHEAR-Assemble to reconstruct the new genome. | |
| 107 </help> | |
| 108 </tool> |