# 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 @@ - - - - - - - - -