comparison tools/mira4/mira4_mapping.xml @ 13:7fcabeeca5df draft

Uploaded v0.0.2 preview 5, fixes for MIRA 4.0 (final), more verbose error if $MIRA4 path wrong
author peterjc
date Sun, 16 Feb 2014 16:30:48 -0500
parents 302d13490b23
children b0ffe0e7282b
comparison
equal deleted inserted replaced
12:a6a56440567c 13:7fcabeeca5df
78 </outputs> 78 </outputs>
79 <configfiles> 79 <configfiles>
80 <configfile name="manifest"> 80 <configfile name="manifest">
81 project = MIRA 81 project = MIRA
82 job = mapping,${job_type},${job_quality} 82 job = mapping,${job_type},${job_quality}
83 parameters = -GE:not=1 -NW:cmrnl -DI:trt=/tmp 83 parameters = -NW:cmrnl=no -DI:trt=/tmp
84 ## -GE:not is short for -GENERAL:number_of_threads and using one (1) 84 ## -GE:not is short for -GENERAL:number_of_threads and using one (1)
85 ## can be useful for repeatability of assemblies and bug hunting. 85 ## can be useful for repeatability of assemblies and bug hunting.
86 ## This is overriden by the command line -t switch which is easier
87 ## to set from within Galaxy.
86 ## 88 ##
87 ## -NW:cmrnl is short for -NAG_AND_WARN:check_maxreadnamelength 89 ## -NW:cmrnl is short for -NAG_AND_WARN:check_maxreadnamelength
88 ## and without this MIRA aborts with read names over 40 characters 90 ## and without this MIRA aborts with read names over 40 characters
89 ## due to limitations of some downstream tools. 91 ## due to limitations of some downstream tools.
90 ## 92 ##
126 ##This is perhaps redundant as MIRA defaults to StrainX for the reads: 128 ##This is perhaps redundant as MIRA defaults to StrainX for the reads:
127 strain = StrainX 129 strain = StrainX
128 #end if 130 #end if
129 ##Record the segment placement (if any) 131 ##Record the segment placement (if any)
130 #if str($rg.segments.type) == "paired" 132 #if str($rg.segments.type) == "paired"
131 segmentplacement = ${rg.segments.placement} 133 segment_placement = ${rg.segments.placement}
132 segmentnaming = ${rg.segments.naming} 134 segment_naming = ${rg.segments.naming}
133 #end if 135 #end if
134 #if str($rg.segments.type) == "none" 136 ##if str($rg.segments.type) == "none"
135 segmentplacement = ? 137 ##MIRA4 manual says use segment_placement = unknown or ? for unpaired data
136 #end if 138 ##but this stopped working in MIRA 4.0 RC5 and 4.0 (final). See:
139 ##http://www.freelists.org/post/mira_talk/Unpaired-reads-and-segment-placement--or-unknown
140 ##segment_placement = ?
141 ##end if
137 ##MIRA will accept multiple filenames on one data line, or multiple data lines 142 ##MIRA will accept multiple filenames on one data line, or multiple data lines
138 #for $f in $rg.filenames 143 #for $f in $rg.filenames
139 ##Must now map Galaxy datatypes to MIRA file types... 144 ##Must now map Galaxy datatypes to MIRA file types...
140 #if $f.ext.startswith("fastq") 145 #if $f.ext.startswith("fastq")
141 ##MIRA doesn't like fastqsanger etc, just plain old fastq: 146 ##MIRA doesn't like fastqsanger etc, just plain old fastq: