diff defuse.xml @ 22:68494d6aabeb

Update datamanager and defuse.xml config file generation
author Jim Johnson <jj@umn.edu>
date Thu, 27 Jun 2013 13:11:50 -0500
parents 1af6f32ff592
children e8fc5de0578b
line wrap: on
line diff
--- a/defuse.xml	Mon Jun 24 08:12:05 2013 -0500
+++ b/defuse.xml	Thu Jun 27 13:11:50 2013 -0500
@@ -20,7 +20,7 @@
       <param name="index" type="select" label="Select a Reference Dataset" help="if your genome of interest is not listed - contact Galaxy team">
         <options from_file="defuse_reference.loc">
           <column name="name" index="1"/>
-          <column name="value" index="2"/>
+          <column name="value" index="3"/>
           <filter type="sort_by" column="0" />
           <validator type="no_options" message="No indexes are available" />
         </options>
@@ -77,6 +77,9 @@
                but they require considerable diskspace, and should be deleted and purged when no longer needed."/>
   <param name="do_get_reads" type="boolean" checked="true" truevalue="yes" falsevalue="no" label="Run get_reads on each cluster"/>
  </inputs>
+ <stdio>
+   <exit_code range="1:"  level="fatal" description="Error Running Defuse" />
+ </stdio>
  <outputs>
   <data format="txt" name="config_txt" label="${tool.name} on ${on_string}: config.txt"/>
   <data format="txt" name="defuse_log" label="${tool.name} on ${on_string}: defuse.log" />
@@ -99,12 +102,11 @@
 #set $ds = chr(36)
 #if $refGenomeSource.genomeSource == "history":
 #set config_file = $refGenomeSource.config.__str__
-#set 
 #else 
 #set config_file = $refGenomeSource.index.value
 #end if
 #set pat = '^\s*([^#=][^=]*?)\s*=\s*(.*?)\s*$'
-#set fh = open()
+#set fh = open($config_file)
 #set keys = ['dataset_directory','ensembl_organism','ensembl_prefix','ensembl_version','ensembl_genome_version','ucsc_genome_version','ncbi_organism','ncbi_prefix','chromosomes','mt_chromosome','gene_sources','ig_gene_sources','rrna_gene_sources']
 #set kv = []
 #for $line in $fh:
@@ -146,18 +148,68 @@
 /project/db/genomes/Hsapiens/hg19/defuse
 #end try
 
+# Organism IDs
+ensembl_organism = #slurp
+#try
+$ref_dict['ensembl_organism']
+#except
+homo_sapiens
+#end try
+
+ensembl_prefix = #slurp
+#try
+$ref_dict['ensembl_prefix']
+#except
+Homo_sapiens
+#end try
+
+ensembl_version = #slurp
+#try
+$ref_dict['ensembl_version']
+#except
+71
+#end try
+
+ensembl_genome_version = #slurp
+#try
+$ref_dict['ensembl_genome_version']
+#except
+GRCh37
+#end try
+
+ucsc_genome_version = #slurp
+#try
+$ref_dict['ucsc_genome_version']
+#except
+hg19
+#end try
+
+ncbi_organism = #slurp
+#try
+$ref_dict['ncbi_organism']
+#except
+Homo_sapiens
+#end try
+
+ncbi_prefix = #slurp
+#try
+$ref_dict['ncbi_prefix']
+#except
+Hs
+#end try
+
 # Input genome and gene models
 gene_models = #slurp
 #try
 $ref_dict['gene_models']
 #except
-\$(dataset_directory)/Homo_sapiens.GRCh37.62.gtf
+\$(dataset_directory)/\$(ensembl_prefix).\$(ensembl_genome_version).\$(ensembl_version).gtf
 #end try
 genome_fasta = #slurp
 #try
 $ref_dict['genome_fasta']
 #except
-\$(dataset_directory)/Homo_sapiens.GRCh37.62.dna.chromosome.fa
+\$(dataset_directory)/\$(ensembl_prefix).\$(ensembl_genome_version).\$(ensembl_version).dna.chromosomes.fa
 #end try
 
 # Repeat table from ucsc genome browser