# HG changeset patch # User iuc # Date 1619784565 0 # Node ID 57251c760cab375068ac82928349c2ce28d2462c # Parent 382518f24d6d33be8739e0c59568c8865ebd51d3 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/ena_upload commit ffea061c1ad6e7291abfe220230dbdbe8d19a2bd" diff -r 382518f24d6d -r 57251c760cab ena_upload.xml --- a/ena_upload.xml Sat Nov 28 09:45:44 2020 +0000 +++ b/ena_upload.xml Fri Apr 30 12:09:25 2021 +0000 @@ -1,24 +1,31 @@ - + - 0.2.4 + 0.2.7 samples_macros.xml ena-upload-cli xlrd + + + + '$output'; #end if #if $action_options.input_format_conditional.input_format != "user_generated_tables": - cp $studies_table_path $studies_table_out; - cp $samples_table_path $samples_table_out; - cp $experiments_table_path $experiments_table_out; - cp $runs_table_path $runs_table_out; - #if $action_options.input_format_conditional.dry_run == "true": - #set $dry_run_option = "True" - #end if + cp $studies_table_path $studies_table_out && + cp $samples_table_path $samples_table_out && + cp $experiments_table_path $experiments_table_out && + cp $runs_table_path $runs_table_out && +#else: + ln -s '$action_options.input_format_conditional.experiments_users_table' $experiments_table_path && + ln -s '$action_options.input_format_conditional.studies_users_table' $studies_table_path && + ln -s '$action_options.input_format_conditional.runs_users_table' $runs_table_path && + ln -s '$action_options.input_format_conditional.samples_users_table' $samples_table_path && #end if -#if $dry_run_option == "False" and $action_options.test_submit == "False": +#if $action_options.test_submit_parameters.dry_run == "false" and $action_options.test_submit == "False": webin_id=`grep 'username' $credentials`; if [ "\$webin_id" = "" ]; then ## No credentials in user defined preferences ## Fallback to global defined credentials (if exist) #import os #if os.path.isfile(os.environ.get('GALAXY_ENA_SECRETS', '')): - credentials_path=\${ENA_SECRETS}; + credentials_path=\${GALAXY_ENA_SECRETS}; webin_id=`grep 'username' \$GALAXY_ENA_SECRETS`; if [ "\$webin_id" = "" ]; then echo "No global credentials defined. Check your GALAXY_ENA_SECRETS file or set your credentials via: User -> Preferences -> Manage Information"; @@ -91,15 +93,41 @@ #end for #end for #else: - #for $file in $action_options.input_format_conditional.data: - #set $safename_reads_file = re.sub('[^\w\-_\.]', '_', $file.element_identifier) - ln -s '$file' $safename_reads_file && - $files_to_upload.append(str($safename_reads_file)) - #end for + #if $action_options.input_format_conditional.run_input_format_conditional.run_input_format == 'paired_list': + #for $pair in $action_options.input_format_conditional.run_input_format_conditional.paired_end_collection: + #set $safename_reads_file = re.sub('[^\w\-_\.]', '_', $pair.name) + #if $pair.forward.is_of_type('fastq.gz', 'fastqsanger.gz'): + #set $safename_fwd_reads_file = $safename_reads_file + '_1.fastq.gz' + #elif $pair.forward.is_of_type('fastqsanger.bz2', 'fastq.bz2'): + #set $safename_fwd_reads_file = $safename_reads_file + '_1.fastq.bz2' + #else: + #set $safename_fwd_reads_file = $safename_reads_file + '_1.fastq' + #end if + #if $pair.reverse.is_of_type('fastq.gz', 'fastqsanger.gz'): + #set $safename_rev_reads_file = $safename_reads_file + '_2.fastq.gz' + #elif $pair.reverse.is_of_type('fastqsanger.bz2', 'fastq.bz2'): + #set $safename_rev_reads_file = $safename_reads_file + '_2.fastq.bz2' + #else: + #set $safename_rev_reads_file = $safename_reads_file + '_2.fastq' + #end if + + ln -s '$pair.forward' $safename_fwd_reads_file && + $files_to_upload.append(str($safename_fwd_reads_file)) + ln -s '$pair.reverse' $safename_rev_reads_file && + $files_to_upload.append(str($safename_rev_reads_file)) + #end for + #end if + #if $action_options.input_format_conditional.run_input_format_conditional.run_input_format == 'multiple_selection_list': + #for $file in $action_options.input_format_conditional.run_input_format_conditional.data: + #set $safename_reads_file = re.sub('[^\w\-_\.]', '_', $file.element_identifier) + ln -s '$file' $safename_reads_file && + $files_to_upload.append(str($safename_reads_file)) + #end for + #end if #end if -#if $dry_run_option == "False": +#if $action_options.test_submit_parameters.dry_run == "false": ena-upload-cli --tool 'ena-upload-cli v@VERSION@ @ Galaxy' --action '$action_options.action' @@ -109,19 +137,15 @@ #for $dataset in $files_to_upload: '$dataset' #end for +--experiment '$experiments_table_path' +--study '$studies_table_path' +--run '$runs_table_path' +--sample '$samples_table_path' #if $action_options.input_format_conditional.input_format == "user_generated_tables": - --experiment '$action_options.input_format_conditional.experiments_users_table' - --study '$action_options.input_format_conditional.studies_users_table' - --run '$action_options.input_format_conditional.runs_users_table' - --sample '$action_options.input_format_conditional.samples_users_table' #if "$action_options.input_format_conditional.viral_submission" == "true": --vir #end if #else: - --experiment '$experiments_table_path' - --study '$studies_table_path' - --run '$runs_table_path' - --sample '$samples_table_path' #if $action_options.input_format_conditional.input_format == "build_tables": #if $action_options.input_format_conditional.conditional_viral_metadata.viral_sample == "true": --vir @@ -133,12 +157,19 @@ #end if #end if -#if $action_options.submit_dev == "true": +#if $action_options.test_submit_parameters.submit_dev == "true": -d #end if - > '$output' + >> '$output'; + echo -e 'center_name\t$action_options.center' >> '$output'; + echo -e 'action_option\t$action_options.action' >> '$output'; + #if $action_options.input_format_conditional.input_format != "user_generated_tables": + cp $studies_table_path_updated $studies_table_out 2>/dev/null; + cp $samples_table_path_updated $samples_table_out 2>/dev/null; + cp $experiments_table_path_updated $experiments_table_out 2>/dev/null; + cp $runs_table_path_updated $runs_table_out 2>/dev/null; + #end if #else: - echo "" exit 0; #end if @@ -170,9 +201,9 @@ #set $safename_reads_file = re.sub('[^\w\-_\.]', '_', $file.element_identifier) $run_files.append(str($safename_reads_file)) #end for - $runs.append($run_files) + $runs.append((str($run.run_base_name),$run_files)) #end for - + $experiments.append({'title':str($experiment.experiment_title),'experiment_design':str($experiment.experiment_design),'library_strategy':str($experiment.library_strategy),'library_source':str($experiment.library_source),'library_selection':str($experiment.library_selection),'library_layout':str($experiment.library_layout),'insert_size':str($experiment.insert_size),'library_construction_protocol':str($experiment.library_construction_protocol),'platform':str($experiment.platform),'instrument_model':str($experiment.instrument_model),'runs':$runs}) #end for #if $action_options.input_format_conditional.conditional_viral_metadata.viral_sample == "true": @@ -196,17 +227,18 @@ - + + - + action_options['input_format_conditional']['input_format'] == "build_tables" or action_options['input_format_conditional']['input_format'] == "excel_tables" @@ -225,13 +257,18 @@ - +
+ + +
- - + + + +
@@ -271,13 +308,18 @@ - +
+ + +
- - + + + +
@@ -311,14 +353,50 @@
+ + + + +
+ + +
+ + + + + + + + + + + + + + + + + + +
+ + + + + + +
- +
+ + +
- @@ -383,12 +461,14 @@ - +
+ + +
- - + @@ -428,13 +508,15 @@ - +
+ + +
- - + @@ -457,6 +539,7 @@ + @@ -472,18 +555,20 @@ - +
- +
+ + +
- @@ -518,6 +603,7 @@ + @@ -534,7 +620,7 @@ - +
diff -r 382518f24d6d -r 57251c760cab extract_tables.py --- a/extract_tables.py Sat Nov 28 09:45:44 2020 +0000 +++ b/extract_tables.py Fri Apr 30 12:09:25 2021 +0000 @@ -3,6 +3,11 @@ import pathlib from datetime import datetime +""" +Parse the configfile generated by the Galaxy tool. +This file is JSON-formatted and should be converted to a set of tabular files. +""" + FILE_FORMAT = 'fastq' parser = argparse.ArgumentParser() @@ -78,11 +83,15 @@ 'submission_date_ENA']) + '\n') run_index = 0 # exp['runs'] is a list of lists - for run in exp['runs']: + for (base_run, run_files) in exp['runs']: run_index += 1 - run_alias = '.'.join(['run_' + str(run_index), str(exp_index), str(sample_index)]) \ - + '_' + timestamp - for file_entry in run: + if base_run != '': + run_alias = base_run + else: + # no alias provided, generated a unique one + run_alias = '_'.join(['run_' + str(run_index), str(exp_index), + str(sample_index)]) + '_' + timestamp + for file_entry in run_files: runs_table.write('\t'.join([run_alias, action, 'ena_run_accession', exp_alias, file_entry, FILE_FORMAT, 'file_checksum', 'submission_date_ENA']) + '\n') diff -r 382518f24d6d -r 57251c760cab process_xlsx.py --- a/process_xlsx.py Sat Nov 28 09:45:44 2020 +0000 +++ b/process_xlsx.py Fri Apr 30 12:09:25 2021 +0000 @@ -3,7 +3,7 @@ import sys import xlrd - +import yaml FILE_FORMAT = 'fastq' @@ -36,15 +36,45 @@ sheet_col_index = sheet_columns[expected_columns[col]] row_dict[expected_columns[col]] = xl_sheet.cell(row_id, sheet_col_index).value # should check for duplicate alias/ids? - data_dict[xl_sheet.cell(row_id, index_col).value] = row_dict + if xl_sheet.cell(row_id, index_col).value in data_dict.keys(): + tmp = data_dict[xl_sheet.cell(row_id, index_col).value] + data_dict[xl_sheet.cell(row_id, index_col).value] = [tmp] + data_dict[xl_sheet.cell(row_id, index_col).value].append(row_dict) + else: + data_dict[xl_sheet.cell(row_id, index_col).value] = row_dict return data_dict +def paste_xls2yaml(xlsx_path): + print('YAML -------------') + xls = xlrd.open_workbook(xlsx_path) + content_dict = {} + for sheet_name in xls.sheet_names(): + if sheet_name == 'controlled_vocabulary': + continue + xls_sheet = xls.sheet_by_name(sheet_name) + sheet_contents_dict = {} + colnames = [] + for col in range(xls_sheet.ncols): + colnames.append(xls_sheet.cell(0, col).value) + # skip first 2 rows (column names and suggestions) + for row_id in range(2, xls_sheet.nrows): + row_dict = {} + for col_id in range(0, xls_sheet.ncols): + row_dict[colnames[col_id]] = xls_sheet.cell(row_id, col_id).value + # should check for duplicate alias/ids? + sheet_contents_dict[row_id] = row_dict + content_dict[sheet_name] = sheet_contents_dict + yaml.dump(content_dict, sys.stdout) + print('YAML -------------') + + parser = argparse.ArgumentParser() parser.add_argument('--form', dest='xlsx_path', required=True) parser.add_argument('--out_dir', dest='out_path', required=True) parser.add_argument('--action', dest='action', required=True) parser.add_argument('--vir', dest='viral_submission', required=False, action='store_true') +parser.add_argument('--verbose', dest='verbose', required=False, action='store_true') args = parser.parse_args() xl_workbook = xlrd.open_workbook(args.xlsx_path) @@ -77,9 +107,11 @@ xl_sheet = xl_workbook.sheet_by_name('ENA_experiment') if xl_sheet.nrows < 3: raise ValueError('No experiments found in experiments sheet') -exp_columns = ['alias', 'title', 'study_alias', 'sample_alias', 'design_description', 'library_name', - 'library_strategy', 'library_source', 'library_selection', 'library_layout', - 'insert_size', 'library_construction_protocol', 'platform', 'instrument_model'] +exp_columns = ['alias', 'title', 'study_alias', 'sample_alias', 'design_description', + 'library_name', 'library_strategy', 'library_source', 'library_selection', + 'library_layout', 'insert_size', 'library_construction_protocol', + 'platform', 'instrument_model'] + experiments_dict = extract_data(xl_sheet, exp_columns) # PARSE RUNS SHEET @@ -123,6 +155,8 @@ # ADD A TIMESTAMP TO THE ALIAS? SEEMS LIKE ENA REQUIRES ALL ENTRIES FOR A WEBIN TO HAVE UNIQUE IDS? # dt_oobj = datetime.now(tz=None) # timestamp = dt_oobj.strftime("%Y%m%d_%H:%M:%S") +runs_included = [] +exp_included = [] for study_alias, study in studies_dict.items(): # study_alias = study_alias + '_' + timestamp studies_table.write('\t'.join([study_alias, action, 'ENA_accession', study['title'], @@ -162,12 +196,37 @@ exp['library_construction_protocol'], exp['platform'], exp['instrument_model'], 'submission_date_ENA']) + '\n') + exp_included.append(exp_alias) for run_alias, run in runs_dict.items(): - if run['experiment_alias'] == exp_alias: - runs_table.write('\t'.join([run_alias, action, 'ena_run_accession', exp_alias, - run['file_name'], FILE_FORMAT, 'file_checksum', - 'submission_date_ENA']) + '\n') + # check that the experiments library_layout is set to paired + # when multiple entries are associated with the same run alias + if not isinstance(run, list): + runs_list = [run] + else: + runs_list = run + for run_entry in runs_list: + if run_entry['experiment_alias'] == exp_alias: + runs_table.write('\t'.join([run_alias, action, 'ena_run_accession', + exp_alias, run_entry['file_name'], + FILE_FORMAT, 'file_checksum', + 'submission_date_ENA']) + '\n') + runs_included.append(run_alias) + +# check if any experiment or run was not associated with any sample +for run in runs_dict.keys(): + if run not in runs_included: + print(f'The run {run} is listed in the runs section but not associated with any \ + used experiment') + +for exp in experiments_dict.keys(): + if exp not in exp_included: + print(f'The experiment {exp} is listed in the experiments section but not associated \ + with any used sample') + studies_table.close() samples_table.close() experiments_table.close() runs_table.close() + +if args.verbose: + paste_xls2yaml(args.xlsx_path) diff -r 382518f24d6d -r 57251c760cab samples_macros.xml --- a/samples_macros.xml Sat Nov 28 09:45:44 2020 +0000 +++ b/samples_macros.xml Fri Apr 30 12:09:25 2021 +0000 @@ -1,5 +1,24 @@ - + +
+ + +
+
+ + + + + + + + + + + + + + @@ -9,22 +28,20 @@ - - - + + - + - - + @@ -51,8 +68,11 @@ - - (^[0-9]{4}(-[0-9]{2}(-[0-9]{2}(T[0-9]{2}:[0-9]{2}(:[0-9]{2})?Z?([+-][0-9]{1,2})?)?)?)?(/[0-9]{4}(-[0-9]{2}(-[0-9]{2}(T[0-9]{2}:[0-9]{2}(:[0-9]{2})?Z?([+-][0-9]{1,2})?)?)?)?)?$)|(^not colected$)|(^not provided$)|(^restricted access$) + + + + + (^[0-9]{4}(-[0-9]{2}(-[0-9]{2}(T[0-9]{2}:[0-9]{2}(:[0-9]{2})?Z?([+-][0-9]{1,2})?)?)?)?(/[0-9]{4}(-[0-9]{2}(-[0-9]{2}(T[0-9]{2}:[0-9]{2}(:[0-9]{2})?Z?([+-][0-9]{1,2})?)?)?)?)?$)|(^not collected$)|(^not provided$)|(^restricted access$) @@ -118,6 +138,7 @@ +
@@ -166,7 +187,7 @@ - + @@ -182,6 +203,7 @@ +
diff -r 382518f24d6d -r 57251c760cab test-data/2.fastqsanger.gz Binary file test-data/2.fastqsanger.gz has changed