# HG changeset patch
# User galaxyp
# Date 1484157966 18000
# Node ID ffed32e397cb34497bedd78958a1019a1aa6ad30
# Parent 2f8e9d0351a81ce19bd0c0f820018feb8e52496a
Uploaded
diff -r 2f8e9d0351a8 -r ffed32e397cb macros.xml
--- a/macros.xml Fri Jun 10 09:57:31 2016 -0400
+++ b/macros.xml Wed Jan 11 13:06:06 2017 -0500
@@ -4,6 +4,7 @@
+
@@ -25,8 +26,8 @@
-fi $forward_ion
-ri $reverse_ion
- 2
- 2.9
+ 3
+ 3.1.4
diff -r 2f8e9d0351a8 -r ffed32e397cb peptide_shaker.xml
--- a/peptide_shaker.xml Fri Jun 10 09:57:31 2016 -0400
+++ b/peptide_shaker.xml Wed Jan 11 13:06:06 2017 -0500
@@ -1,4 +1,4 @@
-
+
Perform protein identification using various search engines based on results from SearchGUI
@@ -6,7 +6,7 @@
macros.xml
- peptide_shaker
+ peptide-shaker
@@ -15,36 +15,30 @@
#set $exp_str = "Galaxy_Experiment_%s" % datetime.now().strftime("%Y%m%d%H%M%s")
#set $samp_str = "Sample_%s" % datetime.now().strftime("%Y%m%d%H%M%s")
#set $temp_stderr = "peptideshaker_stderr"
+ #set $bin_dir = "bin"
mkdir output_reports;
cwd=`pwd`;
- cp -r "\${PEPTIDESHAKER_JAR_PATH%/*}" bin;
- tmp_peptideshaker_jar_path=`echo "\$cwd/bin/\${PEPTIDESHAKER_JAR_PATH\#\#/*/}"`;
-
ln -s '$searchgui_input' searchgui_input.zip &&
+ jar xvf searchgui_input.zip SEARCHGUI_IdentificationParameters.par &&
######################
## PeptideShakerCLI ##
######################
- (java -cp \$tmp_peptideshaker_jar_path eu.isas.peptideshaker.cmd.PeptideShakerCLI
+ (peptide-shaker eu.isas.peptideshaker.cmd.PeptideShakerCLI
+ --exec_dir="\$cwd/${bin_dir}"
-temp_folder \$cwd/PeptideShakerCLI
-experiment '$exp_str'
-sample '$samp_str'
-replicate 1
-identification_files \$cwd/searchgui_input.zip
+ -id_params \$cwd/SEARCHGUI_IdentificationParameters.par
-out \$cwd/peptideshaker_output.cpsx
-zip \$cwd/peptideshaker_output.zip
-threads "\${GALAXY_SLOTS:-12}"
- ##Optional gene annotation parameter
- #if $species_type.species_type_selector != 'no_species_type':
- -species_type "${species_type.species_type_selector}"
- -species "${species_type.species}"
- -species_update 1
- #end if
-
##Optional processing parameters:
#if $processing_options.processing_options_selector == "yes"
-protein_fdr "${processing_options.protein_fdr}"
@@ -85,16 +79,33 @@
#if 'mzidentML' in str($outputs).split(','):
echo "Generating mzIdentML";
- (java -cp \$tmp_peptideshaker_jar_path eu.isas.peptideshaker.cmd.MzidCLI
+ (peptide-shaker eu.isas.peptideshaker.cmd.MzidCLI
+ --exec_dir="\$cwd/${bin_dir}"
-in \$cwd/peptideshaker_output.zip
-output_file output.mzid
- -contact_first_name "Proteomics"
- -contact_last_name "Galaxy"
- -contact_email "galaxyp@umn.edu"
- -contact_address "galaxyp@umn.edu"
- -organization_name "University of Minnesota"
- -organization_email "galaxyp@umn.edu"
- -organization_address "Minneapolis, MN 55455, Vereinigte Staaten"
+ #if $contact_options.contact_options_selector == "yes":
+ -contact_first_name "$contact_options.contact_first_name"
+ -contact_last_name "$contact_options.contact_last_name"
+ -contact_email "$contact_options.contact_email"
+ -contact_address "$contact_options.contact_address"
+ #if str($contact_options.contact_url).strip() != '':
+ -contact_url = "$contact_options.contact_url"
+ #end if
+ -organization_name "$contact_options.organization_name"
+ -organization_email "$contact_options.organization_email"
+ -organization_address "$contact_options.organization_address"
+ #if str($contact_options.organization_url).strip() != '':
+ -organization_url = "$contact_options.organization_url"
+ #end if
+ #else:
+ -contact_first_name "Proteomics"
+ -contact_last_name "Galaxy"
+ -contact_email "galaxyp@umn.edu"
+ -contact_address "galaxyp@umn.edu"
+ -organization_name "University of Minnesota"
+ -organization_email "galaxyp@umn.edu"
+ -organization_address "Minneapolis, MN 55455, Vereinigte Staaten"
+ #end if
2>> $temp_stderr)
&&
#end if
@@ -104,7 +115,8 @@
## and will not be passed to the command line
#if set(["0","1","2","3","4","5","6","7"]).intersection( set( str( $outputs ).split(',') ) ):
- (java -cp \$tmp_peptideshaker_jar_path eu.isas.peptideshaker.cmd.ReportCLI
+ (peptide-shaker eu.isas.peptideshaker.cmd.ReportCLI
+ --exec_dir="\$cwd/${bin_dir}"
-temp_folder \$cwd/ReportCLI
-in \$cwd/peptideshaker_output.zip
-out_reports \$cwd/output_reports
@@ -125,6 +137,7 @@
&&
#end if
+ ## # ls -l \$cwd/output_reports/* ;
#if '0' in str($outputs).split(','):
find \$cwd/output_reports -name '*Certificate_of_Analysis*' -exec bash -c 'mv "$0" "certificate.txt"' {} \;
@@ -138,6 +151,10 @@
find \$cwd/output_reports -name '*PSM_Phosphorylation_Report*' -exec bash -c 'mv "$0" "psm_phospho.txt"' {} \;
;
#end if
+ #if '8' in str($outputs).split(','):
+ find \$cwd/output_reports -name '*Extended_PSM_Report*' -exec bash -c 'mv "$0" "psmx.txt"' {} \;
+ ;
+ #end if
#if '3' in str($outputs).split(','):
find \$cwd/output_reports -name '*PSM_Report*' -exec bash -c 'mv "$0" "psm.txt"' {} \;
;
@@ -172,287 +189,6 @@
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
@@ -503,15 +239,53 @@
+
+
+
+
+
+
+
+
+ \S+.*
+
+
+ \S+.*
+
+
+ \S+@\S+
+
+
+ \S+.*
+
+
+
+
+ \S+.*
+
+
+ \S+@\S+
+
+
+ \S+.*
+
+
+
+
+
+
+
-
-
-
-
-
+
+
+
+
+
+
@@ -541,6 +315,9 @@
'3' in outputs
+
+ '8' in outputs
+
'4' in outputs
diff -r 2f8e9d0351a8 -r ffed32e397cb searchgui.xml
--- a/searchgui.xml Fri Jun 10 09:57:31 2016 -0400
+++ b/searchgui.xml Wed Jan 11 13:06:06 2017 -0500
@@ -7,7 +7,9 @@
searchgui
+
@@ -16,21 +18,25 @@
#set $exp_str = "Galaxy_Experiment_%s" % datetime.now().strftime("%Y%m%d%H%M%s")
#set $samp_str = "Sample_%s" % datetime.now().strftime("%Y%m%d%H%M%s")
#set $temp_stderr = "searchgui_stderr"
+ #set $bin_dir = "bin"
mkdir output;
mkdir output_reports;
cwd=`pwd`;
export HOME=\$cwd;
+ ## echo the search engines to run
+ echo "$engines";
+ echo "DB: ${input_database.display_name} sequences: ${input_database.metadata.sequences}";
+
##Create a searchgui.properties file for the version, which will be added to the searchgui_results if not already present
echo "searchgui.version=@SEARCHGUI_VERSION@" >> searchgui.properties;
- cp -r "\${SEARCHGUI_JAR_PATH%/*}" bin;
- tmp_searchgui_jar_path=`echo "\$cwd/bin/\${SEARCHGUI_JAR_PATH\#\#/*/}"`;
-
#for $mgf in $peak_lists:
#set $input_name = $mgf.display_name.replace(".mgf", "") + ".mgf"
ln -s -f '${mgf}' '${input_name}';
+ #set $encoded_id = $__app__.security.encode_id($mgf.id)
+ echo "Spectrums:${mgf.display_name}(API:${encoded_id}) ";
#end for
##ln -s "${input_database}" input_database.fasta;
cp "${input_database}" input_database.fasta;
@@ -40,7 +46,7 @@
###########################################
#if $create_decoy:
echo "Creating decoy database.";
- java -cp \$tmp_searchgui_jar_path eu.isas.searchgui.cmd.FastaCLI -in input_database.fasta -decoy &&
+ searchgui eu.isas.searchgui.cmd.FastaCLI --exec_dir="\$cwd/${bin_dir}" -in input_database.fasta -decoy &&
rm input_database.fasta &&
cp input_database_concatenated_target_decoy.fasta input_database.fasta &&
## ln -sf input_database_concatenated_target_decoy.fasta input_database.fasta;
@@ -49,8 +55,8 @@
#####################################################
## generate IdentificationParameters for SearchGUI ##
#####################################################
-
- (java -cp \$tmp_searchgui_jar_path eu.isas.searchgui.cmd.IdentificationParametersCLI
+ (searchgui eu.isas.searchgui.cmd.IdentificationParametersCLI
+ --exec_dir="\$cwd/${bin_dir}"
-out SEARCHGUI_IdentificationParameters.par
@GENERAL_PARAMETERS@
@@ -196,6 +202,7 @@
-andromeda_min_pep_length ${andromeda.andromeda_min_pep_length}
-andromeda_max_pep_length ${andromeda.andromeda_max_pep_length}
-andromeda_max_psms ${andromeda.andromeda_max_psms}
+ -andromeda_decoy_mode ${andromeda.andromeda_decoy_mode}
#end if
*#
@@ -245,15 +252,16 @@
#if $comet.comet_spectrum.comet_spectrum_selector == "yes"
-comet_min_peaks ${comet.comet_spectrum.comet_min_peaks}
-comet_min_peak_int ${comet.comet_spectrum.comet_min_peak_int}
- -comet_remove_prec ${comet.comet_spectrum.comet_remove_prec}
+
+ -comet_remove_prec ${comet.comet_spectrum.comet_prec.comet_remove_prec}
- #if $comet.comet_spectrum.comet_remove_prec == "1"
- -comet_remove_prec_tol ${comet.comet_spectrum.comet_remove_prec_tol}
+ #if $comet.comet_spectrum.comet_prec.comet_remove_prec == "1"
+ -comet_remove_prec_tol ${comet.comet_spectrum.comet_prec.comet_remove_prec_tol}
#end if
- #if $comet.comet_spectrum.comet_remove_prec == "2"
- -comet_remove_prec_tol ${comet.comet_spectrum.comet_remove_prec_tol}
+ #if $comet.comet_spectrum.comet_prec.comet_remove_prec == "2"
+ -comet_remove_prec_tol ${comet.comet_spectrum.comet_prec.comet_remove_prec_tol}
#end if
-comet_clear_mz_range_lower ${comet.comet_spectrum.comet_clear_mz_range_lower}
@@ -279,13 +287,44 @@
#end if
#end if
+ #if $directtag.directtag_advanced == "yes"
+ -directag_tic_cutoff ${directtag.directag_tic_cutoff}
+ -directag_max_peak_count ${directtag.directag_max_peak_count}
+ -directag_intensity_classes ${directtag.directag_intensity_classes}
+ -directag_adjust_precursor ${directtag.directag_adjust_precursor}
+ -directag_min_adjustment ${directtag.directag_min_adjustment}
+ -directag_max_adjustment ${directtag.directag_max_adjustment}
+ -directag_adjustment_step ${directtag.directag_adjustment_step}
+ -directag_charge_states ${directtag.directag_charge_states}
+ #if str($directtag.directag_output_suffix).strip() != '':
+ -directag_output_suffix ${directtag.directag_output_suffix}
+ #end if
+ -directag_ms_charge_state ${directtag.directag_ms_charge_state}
+ -directag_duplicate_spectra ${directtag.directag_duplicate_spectra}
+ -directag_deisotoping ${directtag.directag_deisotoping}
+ -directag_isotope_tolerance ${directtag.directag_isotope_tolerance}
+ -directag_complement_tolerance ${directtag.directag_complement_tolerance}
+ -directag_tag_length ${directtag.directag_tag_length}
+ -directag_max_var_mods ${directtag.directag_max_var_mods}
+ -directag_max_tag_count ${directtag.directag_max_tag_count}
+ -directag_intensity_weight ${directtag.directag_intensity_weight}
+ -directag_fidelity_weight ${directtag.directag_fidelity_weight}
+ -directag_complement_weight ${directtag.directag_complement_weight}
+ #end if
+
+ #if $novor.novor_advanced == "yes"
+ -novor_fragmentation ${novor.novor_fragmentation}
+ -novor_mass_analyzer ${novor.novor_mass_analyzer}
+ #end if
+
2> $temp_stderr)
&&
################
## Search CLI ##
################
- (java -Djava.awt.headless=true -cp \$tmp_searchgui_jar_path eu.isas.searchgui.cmd.SearchCLI
+ (searchgui -Djava.awt.headless=true eu.isas.searchgui.cmd.SearchCLI
+ --exec_dir="\$cwd/${bin_dir}"
-temp_folder `pwd`
-spectrum_files \$cwd
-output_folder \$cwd/output
@@ -355,6 +394,18 @@
-andromeda 0
#end if
+ #if 'Novor' in $engines_list:
+ -novor 1
+ #else
+ -novor 0
+ #end if
+
+ #if 'DirecTag' in $engines_list:
+ -directag 1
+ #else
+ -directag 0
+ #end if
+
## single zip file
-output_option 0
@@ -396,19 +447,26 @@
Comet and Tide shouldn't both be selected since they use a similar algoritm.
-
+
+
+
+
+
+
+
+
@@ -439,7 +497,7 @@
-
+
@@ -470,7 +528,7 @@
+ label="X!Tandem: Maximum Valid Expectation Value" type="float" value="0.01" />
+
-->
@@ -891,20 +950,22 @@
label="Comet: Minimum Number of Peaks per Spectrum" help="The minimum number of peaks per spectrum" />
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff -r 2f8e9d0351a8 -r ffed32e397cb tool_dependencies.xml
--- a/tool_dependencies.xml Fri Jun 10 09:57:31 2016 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
@@ -1,9 +0,0 @@
-
-
-
-
-
-
-
-
-