Repository revision
8:8434f070e939

Repository 'openms_mapaligneridentification'
hg clone https://toolshed.g2.bx.psu.edu/repos/galaxyp/openms_mapaligneridentification

MapAlignerIdentification tool metadata
Miscellaneous
Corrects retention time distortions between maps based on common peptide identifications.
MapAlignerIdentification
toolshed.g2.bx.psu.edu/repos/galaxyp/openms_mapaligneridentification/MapAlignerIdentification/2.8+galaxy0
2.8+galaxy0
None
True
Version lineage of this tool (guids ordered most recent to oldest)
toolshed.g2.bx.psu.edu/repos/galaxyp/openms_mapaligneridentification/MapAlignerIdentification/2.8+galaxy0 (this tool)
toolshed.g2.bx.psu.edu/repos/galaxyp/openms_mapaligneridentification/MapAlignerIdentification/2.6+galaxy0
toolshed.g2.bx.psu.edu/repos/galaxyp/openms_mapaligneridentification/MapAlignerIdentification/2.5+galaxy0
toolshed.g2.bx.psu.edu/repos/galaxyp/openms_mapaligneridentification/MapAlignerIdentification/2.3.0
MapAlignerIdentification
Requirements (dependencies defined in the <requirements> tag set)
name version type
openms 2.8 package
openms-thirdparty 2.8 package
omssa 2.1.9 package
blast 2.13.0 package
ctdopts 1.4 package
Additional information about this tool
#def quote(s):
    #set $s = [ _ for _ in $s.split(" ") if _ != "" ]
    #set $q = False
    #for $i, $p in enumerate($s):
        #if $p == "":
            #continue
        #end if
        #if $p.startswith('"'):
            #set $q = True
        #end if
##        #if p.startswith('-'):
##            #set p = "\\" + p
##        #elif p.startswith('"-'):
##            #set p = "\\" + p[1:]
##        #end if
        #if not $q:
            #set $s[i] = '"%s"' % p
        #end if
        #if $p.endswith('"'):
            #set $q = False
        #end if
    #end for
    #return " ".join($s)
#end def
  
#def oms2gxyext(o)
    #set m={'txt': 'txt', 'tsv': 'tabular', 'bioml': 'xml', 'consensusXML': 'consensusxml', 'csv': 'csv', 'dta': 'dta', 'dta2d': 'dta2d', 'edta': 'edta', 'fa': 'fasta', 'fas': 'fasta', 'fasta': 'fasta', 'FASTA': 'fasta', 'featureXML': 'featurexml', 'featurexml': 'featurexml', 'html': 'html', 'HTML': 'html', 'idXML': 'idxml', 'json': 'json', 'kroenik': 'kroenik', 'mascotXML': 'mascotxml', 'mgf': 'mgf', 'mrm': 'mrm', 'ms': 'sirius.ms', 'ms2': 'ms2', 'msp': 'msp', 'mzData': 'mzdata', 'mzid': 'mzid', 'mzML': 'mzml', 'mzml': 'mzml', 'mzq': 'mzq', 'mzQC': 'mzqc', 'mzTab': 'mztab', 'mzXML': 'mzxml', 'novor': 'txt', 'obo': 'obo', 'oms': 'sqlite', 'omssaXML': 'idxml', 'osw': 'osw', 'OSW': 'osw', 'params': 'txt', 'paramXML': 'paramxml', 'peplist': 'peplist', 'pep.xml': 'pepxml', 'pepXML': 'pepxml', 'png': 'png', 'PNG': 'png', 'protXML': 'protxml', 'psms': 'psms', 'pqp': 'pqp', 'qcML': 'qcml', 'spec.xml': 'spec.xml', 'splib': 'splib', 'sqMass': 'sqmass', 'tandem.xml': 'tandem', 'trafoXML': 'trafoxml', 'traML': 'traml', 'TraML': 'traml', 'tab': 'tabular', 'raw': 'thermo.raw', 'xls': 'tsv', 'XML': 'xml', 'xml': 'xml', 'xquest.xml': 'xquest.xml', 'xsd': 'xml'}
    #return m[o]
#end def
#def gxy2omsext(g)
    #set m={'txt': 'txt', 'tabular': 'tsv', 'xml': 'bioml', 'consensusxml': 'consensusXML', 'csv': 'csv', 'dta': 'dta', 'dta2d': 'dta2d', 'edta': 'edta', 'fasta': 'fa', 'featurexml': 'featureXML', 'html': 'html', 'idxml': 'idXML', 'json': 'json', 'kroenik': 'kroenik', 'mascotxml': 'mascotXML', 'mgf': 'mgf', 'mrm': 'mrm', 'sirius.ms': 'ms', 'ms2': 'ms2', 'msp': 'msp', 'mzdata': 'mzData', 'mzid': 'mzid', 'mzml': 'mzML', 'mzq': 'mzq', 'mzqc': 'mzQC', 'mztab': 'mzTab', 'mzxml': 'mzXML', 'obo': 'obo', 'sqlite': 'oms', 'osw': 'osw', 'paramxml': 'paramXML', 'peff': 'fasta', 'peplist': 'peplist', 'pepxml': 'pep.xml', 'png': 'png', 'protxml': 'protXML', 'psms': 'psms', 'pqp': 'pqp', 'qcml': 'qcML', 'spec.xml': 'spec.xml', 'splib': 'splib', 'sqmass': 'sqMass', 'tandem': 'tandem.xml', 'trafoxml': 'trafoXML', 'traml': 'traML', 'thermo.raw': 'raw', 'tsv': 'xls', 'xquest.xml': 'xquest.xml'}
    #return m[g]
#end def

#import re

## Preprocessing
mkdir in_cond.in &&
#if $in_cond.in_select == "no"
mkdir ${' '.join(["'in_cond.in/%s'" % (i) for i, f in enumerate($in_cond.in) if f])} && 
${' '.join(["ln -s '%s' 'in_cond.in/%s/%s.%s' && " % (f, i, re.sub('[^\w\-_]', '_', f.element_identifier), $gxy2omsext(f.ext)) for i, f in enumerate($in_cond.in) if f])}
#else
ln -s '$in_cond.in' 'in_cond.in/${re.sub("[^\w\-_]", "_", $in_cond.in.element_identifier)}.$gxy2omsext($in_cond.in.ext)' &&
#end if
#if "out_FLAG" in str($OPTIONAL_OUTPUTS).split(',')
  mkdir out &&
  mkdir ${' '.join(["'out/%s'" % (i) for i, f in enumerate($in_cond.in) if f])} && 
#end if
#if "trafo_out_FLAG" in str($OPTIONAL_OUTPUTS).split(',')
  mkdir trafo_out &&
  mkdir ${' '.join(["'trafo_out/%s'" % (i) for i, f in enumerate($in_cond.in) if f])} && 
#end if
#if $design:
  mkdir design &&
  ln -s '$design' 'design/${re.sub("[^\w\-_]", "_", $design.element_identifier)}.$gxy2omsext($design.ext)' &&
#end if
#if $reference.file:
  mkdir reference.file &&
  ln -s '$reference.file' 'reference.file/${re.sub("[^\w\-_]", "_", $reference.file.element_identifier)}.$gxy2omsext($reference.file.ext)' &&
#end if

## Main program call

set -o pipefail &&
MapAlignerIdentification -write_ctd ./ &&
python3 '$__tool_directory__/fill_ctd.py' 'MapAlignerIdentification.ctd' '$args_json' '$hardcoded_json' &&
MapAlignerIdentification -ini MapAlignerIdentification.ctd
-in
#if $in_cond.in_select == "no"
${' '.join(["'in_cond.in/%s/%s.%s'"%(i, re.sub('[^\w\-_]', '_', f.element_identifier), $gxy2omsext(f.ext)) for i, f in enumerate($in_cond.in) if f])}
#else
'in_cond.in/${re.sub("[^\w\-_]", "_", $in_cond.in.element_identifier)}.$gxy2omsext($in_cond.in.ext)'
#end if
#if "out_FLAG" in str($OPTIONAL_OUTPUTS).split(',')
  -out
  ${' '.join(["'out/%s/%s.%s'"%(i, re.sub('[^\w\-_]', '_', f.element_identifier), f.ext) for i, f in enumerate($in_cond.in) if f])}
#end if
#if "trafo_out_FLAG" in str($OPTIONAL_OUTPUTS).split(',')
  -trafo_out
  ${' '.join(["'trafo_out/%s/%s.%s'"%(i, re.sub('[^\w\-_]', '_', f.element_identifier), $gxy2omsext("trafoxml")) for i, f in enumerate($in_cond.in) if f])}
#end if
#if $design:
  -design
  'design/${re.sub("[^\w\-_]", "_", $design.element_identifier)}.$gxy2omsext($design.ext)'
#end if
#if $reference.file:
  -reference:file
  'reference.file/${re.sub("[^\w\-_]", "_", $reference.file.element_identifier)}.$gxy2omsext($reference.file.ext)'
#end if
#if len(str($OPTIONAL_OUTPUTS).split(',')) == 0
  | tee '$stdout'
#end if

## Postprocessing
#if "trafo_out_FLAG" in str($OPTIONAL_OUTPUTS).split(',')
  ${' '.join(["&& mv -n 'trafo_out/%(bn)s/%(id)s.%(gext)s' 'trafo_out/%(bn)s/%(id)s'"%{"bn": i, "id": re.sub('[^\w\-_]', '_', f.element_identifier), "gext": $gxy2omsext("trafoxml")} for i, f in enumerate($in_cond.in) if f])}
#end if
#if "ctd_out_FLAG" in $OPTIONAL_OUTPUTS
  && mv 'MapAlignerIdentification.ctd' '$ctd_out'
#end if
None
False
Functional tests
name inputs outputs required files
Test-1 in_cond|in: ['MapAlignerIdentification_1_input1.featureXML', 'MapAlignerIdentification_1_input2.featureXML']
store_original_rt: False
reference|index: 0
algorithm|score_type:
algorithm|score_cutoff: False
algorithm|min_score: 0.05
algorithm|min_run_occur: 2
algorithm|max_rt_shift: 0.0
algorithm|use_unassigned_peptides: True
algorithm|use_feature_rt: False
algorithm|use_adducts: True
model|type: b_spline
model|linear|symmetric_regression: False
model|linear|x_weight: None
model|linear|y_weight: None
model|linear|x_datum_min: 1e-15
model|linear|x_datum_max: 1000000000000000.0
model|linear|y_datum_min: 1e-15
model|linear|y_datum_max: 1000000000000000.0
model|b_spline|wavelength: 0.0
model|b_spline|num_nodes: 5
model|b_spline|extrapolate: linear
model|b_spline|boundary_condition: 2
model|lowess|span: 0.666666666666667
model|lowess|num_iterations: 3
model|lowess|delta: -1.0
model|lowess|interpolation_type: cspline
model|lowess|extrapolation_type: four-point-linear
model|interpolated|interpolation_type: cspline
model|interpolated|extrapolation_type: two-point-linear
adv_opts|force: False
adv_opts|test: true
OPTIONAL_OUTPUTS: ['ctd_out_FLAG', 'out_FLAG']
name: value
MapAlignerIdentification_1_input1.featureXML
MapAlignerIdentification_1_input2.featureXML
value
Test-2 in_cond|in: MapAlignerIdentification_1_input1.featureXML
store_original_rt: False
reference|file: MapAlignerIdentification_1_input2.featureXML
reference|index: 0
algorithm|score_type:
algorithm|score_cutoff: False
algorithm|min_score: 0.05
algorithm|min_run_occur: 2
algorithm|max_rt_shift: 0.0
algorithm|use_unassigned_peptides: True
algorithm|use_feature_rt: False
algorithm|use_adducts: True
model|type: b_spline
model|linear|symmetric_regression: False
model|linear|x_weight: None
model|linear|y_weight: None
model|linear|x_datum_min: 1e-15
model|linear|x_datum_max: 1000000000000000.0
model|linear|y_datum_min: 1e-15
model|linear|y_datum_max: 1000000000000000.0
model|b_spline|wavelength: 0.0
model|b_spline|num_nodes: 5
model|b_spline|extrapolate: linear
model|b_spline|boundary_condition: 2
model|lowess|span: 0.666666666666667
model|lowess|num_iterations: 3
model|lowess|delta: -1.0
model|lowess|interpolation_type: cspline
model|lowess|extrapolation_type: four-point-linear
model|interpolated|interpolation_type: cspline
model|interpolated|extrapolation_type: two-point-linear
adv_opts|force: False
adv_opts|test: true
OPTIONAL_OUTPUTS: ['ctd_out_FLAG', 'out_FLAG']
name: value
MapAlignerIdentification_1_input1.featureXML
MapAlignerIdentification_1_input2.featureXML
value
Test-3 in_cond|in: ['MapAlignerIdentification_1_input2.featureXML', 'MapAlignerIdentification_1_input1.featureXML']
store_original_rt: False
reference|index: 1
algorithm|score_type:
algorithm|score_cutoff: False
algorithm|min_score: 0.05
algorithm|min_run_occur: 2
algorithm|max_rt_shift: 0.0
algorithm|use_unassigned_peptides: True
algorithm|use_feature_rt: False
algorithm|use_adducts: True
model|type: b_spline
model|linear|symmetric_regression: False
model|linear|x_weight: None
model|linear|y_weight: None
model|linear|x_datum_min: 1e-15
model|linear|x_datum_max: 1000000000000000.0
model|linear|y_datum_min: 1e-15
model|linear|y_datum_max: 1000000000000000.0
model|b_spline|wavelength: 0.0
model|b_spline|num_nodes: 5
model|b_spline|extrapolate: linear
model|b_spline|boundary_condition: 2
model|lowess|span: 0.666666666666667
model|lowess|num_iterations: 3
model|lowess|delta: -1.0
model|lowess|interpolation_type: cspline
model|lowess|extrapolation_type: four-point-linear
model|interpolated|interpolation_type: cspline
model|interpolated|extrapolation_type: two-point-linear
adv_opts|force: False
adv_opts|test: true
OPTIONAL_OUTPUTS: ['ctd_out_FLAG', 'out_FLAG']
name: value
MapAlignerIdentification_1_input2.featureXML
MapAlignerIdentification_1_input1.featureXML
value
Test-4 in_cond|in: ['MapAlignerIdentification_1_input1.featureXML', 'MapAlignerIdentification_1_input2.featureXML']
store_original_rt: False
reference|index: 2
algorithm|score_type:
algorithm|score_cutoff: False
algorithm|min_score: 0.05
algorithm|min_run_occur: 2
algorithm|max_rt_shift: 0.0
algorithm|use_unassigned_peptides: True
algorithm|use_feature_rt: False
algorithm|use_adducts: True
model|type: b_spline
model|linear|symmetric_regression: False
model|linear|x_weight: None
model|linear|y_weight: None
model|linear|x_datum_min: 1e-15
model|linear|x_datum_max: 1000000000000000.0
model|linear|y_datum_min: 1e-15
model|linear|y_datum_max: 1000000000000000.0
model|b_spline|wavelength: 0.0
model|b_spline|num_nodes: 5
model|b_spline|extrapolate: linear
model|b_spline|boundary_condition: 2
model|lowess|span: 0.666666666666667
model|lowess|num_iterations: 3
model|lowess|delta: -1.0
model|lowess|interpolation_type: cspline
model|lowess|extrapolation_type: four-point-linear
model|interpolated|interpolation_type: cspline
model|interpolated|extrapolation_type: two-point-linear
adv_opts|force: False
adv_opts|test: true
OPTIONAL_OUTPUTS: ['ctd_out_FLAG', 'out_FLAG']
name: value
MapAlignerIdentification_1_input1.featureXML
MapAlignerIdentification_1_input2.featureXML
value
Test-5 in_cond|in: ['MapAlignerIdentification_5_input1.consensusXML', 'MapAlignerIdentification_5_input2.consensusXML']
store_original_rt: False
reference|index: 0
algorithm|score_type:
algorithm|score_cutoff: False
algorithm|min_score: 0.05
algorithm|min_run_occur: 2
algorithm|max_rt_shift: 0.0
algorithm|use_unassigned_peptides: True
algorithm|use_feature_rt: False
algorithm|use_adducts: True
model|type: b_spline
model|linear|symmetric_regression: False
model|linear|x_weight: None
model|linear|y_weight: None
model|linear|x_datum_min: 1e-15
model|linear|x_datum_max: 1000000000000000.0
model|linear|y_datum_min: 1e-15
model|linear|y_datum_max: 1000000000000000.0
model|b_spline|wavelength: 0.0
model|b_spline|num_nodes: 5
model|b_spline|extrapolate: linear
model|b_spline|boundary_condition: 2
model|lowess|span: 0.666666666666667
model|lowess|num_iterations: 3
model|lowess|delta: -1.0
model|lowess|interpolation_type: cspline
model|lowess|extrapolation_type: four-point-linear
model|interpolated|interpolation_type: cspline
model|interpolated|extrapolation_type: two-point-linear
adv_opts|force: False
adv_opts|test: true
OPTIONAL_OUTPUTS: ['ctd_out_FLAG', 'out_FLAG']
name: value
MapAlignerIdentification_5_input1.consensusXML
MapAlignerIdentification_5_input2.consensusXML
value
Test-6 in_cond|in: MapAlignerIdentification_1_input1.featureXML
store_original_rt: False
reference|file: MapAlignerIdentification_1_input2.featureXML
reference|index: 0
algorithm|score_type:
algorithm|score_cutoff: False
algorithm|min_score: 0.05
algorithm|min_run_occur: 2
algorithm|max_rt_shift: 0.0
algorithm|use_unassigned_peptides: True
algorithm|use_feature_rt: False
algorithm|use_adducts: True
model|type: b_spline
model|linear|symmetric_regression: False
model|linear|x_weight: None
model|linear|y_weight: None
model|linear|x_datum_min: 1e-15
model|linear|x_datum_max: 1000000000000000.0
model|linear|y_datum_min: 1e-15
model|linear|y_datum_max: 1000000000000000.0
model|b_spline|wavelength: 0.0
model|b_spline|num_nodes: 5
model|b_spline|extrapolate: linear
model|b_spline|boundary_condition: 2
model|lowess|span: 0.666666666666667
model|lowess|num_iterations: 3
model|lowess|delta: -1.0
model|lowess|interpolation_type: cspline
model|lowess|extrapolation_type: four-point-linear
model|interpolated|interpolation_type: cspline
model|interpolated|extrapolation_type: two-point-linear
adv_opts|force: False
adv_opts|test: true
OPTIONAL_OUTPUTS: ['ctd_out_FLAG', 'trafo_out_FLAG']
name: value
MapAlignerIdentification_1_input1.featureXML
MapAlignerIdentification_1_input2.featureXML
value
Test-7 in_cond|in: MapAlignerIdentification_7_input1.idXML
store_original_rt: False
reference|file: MapAlignerIdentification_7_input2.idXML
reference|index: 0
algorithm|score_type:
algorithm|score_cutoff: False
algorithm|min_score: 0.05
algorithm|min_run_occur: 2
algorithm|max_rt_shift: 0.5
algorithm|use_unassigned_peptides: True
algorithm|use_feature_rt: False
algorithm|use_adducts: True
model|type: b_spline
model|linear|symmetric_regression: False
model|linear|x_weight: None
model|linear|y_weight: None
model|linear|x_datum_min: 1e-15
model|linear|x_datum_max: 1000000000000000.0
model|linear|y_datum_min: 1e-15
model|linear|y_datum_max: 1000000000000000.0
model|b_spline|wavelength: 0.0
model|b_spline|num_nodes: 5
model|b_spline|extrapolate: linear
model|b_spline|boundary_condition: 2
model|lowess|span: 0.666666666666667
model|lowess|num_iterations: 3
model|lowess|delta: -1.0
model|lowess|interpolation_type: cspline
model|lowess|extrapolation_type: four-point-linear
model|interpolated|interpolation_type: cspline
model|interpolated|extrapolation_type: two-point-linear
adv_opts|force: True
adv_opts|test: true
OPTIONAL_OUTPUTS: ['ctd_out_FLAG', 'out_FLAG', 'trafo_out_FLAG']
name: value
MapAlignerIdentification_7_input1.idXML
MapAlignerIdentification_7_input2.idXML
value
Test-8 in_cond|in: MapAlignerIdentification_8_input1.oms
store_original_rt: True
reference|file: MapAlignerIdentification_8_input2.oms
reference|index: 0
algorithm|score_type:
algorithm|score_cutoff: False
algorithm|min_score: 0.05
algorithm|min_run_occur: 2
algorithm|max_rt_shift: 0.5
algorithm|use_unassigned_peptides: True
algorithm|use_feature_rt: False
algorithm|use_adducts: True
model|type: b_spline
model|linear|symmetric_regression: False
model|linear|x_weight: None
model|linear|y_weight: None
model|linear|x_datum_min: 1e-15
model|linear|x_datum_max: 1000000000000000.0
model|linear|y_datum_min: 1e-15
model|linear|y_datum_max: 1000000000000000.0
model|b_spline|wavelength: 0.0
model|b_spline|num_nodes: 5
model|b_spline|extrapolate: linear
model|b_spline|boundary_condition: 2
model|lowess|span: 0.666666666666667
model|lowess|num_iterations: 3
model|lowess|delta: -1.0
model|lowess|interpolation_type: cspline
model|lowess|extrapolation_type: four-point-linear
model|interpolated|interpolation_type: cspline
model|interpolated|extrapolation_type: two-point-linear
adv_opts|force: False
adv_opts|test: true
OPTIONAL_OUTPUTS: ['ctd_out_FLAG', 'out_FLAG', 'trafo_out_FLAG']
name: value
MapAlignerIdentification_8_input1.oms
MapAlignerIdentification_8_input2.oms
value