Mercurial > repos > ieguinoa > ena_upload
comparison extract_tables.py @ 14:436da36c6114 draft
Uploaded
author | ieguinoa |
---|---|
date | Tue, 22 Sep 2020 20:36:41 +0000 |
parents | 91a3a617d94b |
children | d85e35a36715 |
comparison
equal
deleted
inserted
replaced
13:91a3a617d94b | 14:436da36c6114 |
---|---|
32 sample_alias = 'sample_'+str(sample_index) | 32 sample_alias = 'sample_'+str(sample_index) |
33 samples_table.write('\t'.join([sample_alias,action,'ena_accession',sample['title'],sample['tax_name'], sample['tax_id'],sample['description'],'ENA_submission_date'])+ '\n') | 33 samples_table.write('\t'.join([sample_alias,action,'ena_accession',sample['title'],sample['tax_name'], sample['tax_id'],sample['description'],'ENA_submission_date'])+ '\n') |
34 for exp_index,exp in enumerate(sample['experiments']): | 34 for exp_index,exp in enumerate(sample['experiments']): |
35 exp_alias = 'experiment_'+str(exp_index)+'_'+str(sample_index) | 35 exp_alias = 'experiment_'+str(exp_index)+'_'+str(sample_index) |
36 lib_alias = 'library_'+str(exp_index)+'_'+str(sample_index) | 36 lib_alias = 'library_'+str(exp_index)+'_'+str(sample_index) |
37 experiments_table.write('\t'.join([exp_alias,action,'accession_ena',exp['title'],study_alias,sample_alias,exp['experiment_design'],lib_alias,exp['library_strategy'],exp['library_source'],exp['library_selection'],exp['library_layout'].lower(),exp['insert_size'],exp['library_construction_protocol'],exp['platform'],exp['instrument_model'],'submission_date_ENA']) + '\n') | 37 experiments_table.write('\t'.join([exp_alias,action,'accession_ena',exp['title'],study_alias,sample_alias,exp['experiment_design'],lib_alias,exp['library_strategy'],exp['library_source'],exp['library_selection'],exp['library_layout'],exp['insert_size'],exp['library_construction_protocol'],exp['platform'],exp['instrument_model'],'submission_date_ENA']) + '\n') |
38 run_index = 0 | 38 run_index = 0 |
39 # exp['runs'] is a list of lists | 39 # exp['runs'] is a list of lists |
40 for run in exp['runs']: | 40 for run in exp['runs']: |
41 run_index += 1 | 41 run_index += 1 |
42 run_alias = '_'.join(['run',str(exp_index),str(sample_index),str(run_index)]) | 42 run_alias = '_'.join(['run',str(exp_index),str(sample_index),str(run_index)]) |
43 for file_entry in run: | 43 for file_entry in run: |
44 file_format = 'fastq' | 44 file_format = 'fastq.gz' if os.path.splitext(file_entry)[-1] == '.gz' else 'fastq.bz2' |
45 runs_table.write('\t'.join([run_alias,action,'ena_run_accession',exp_alias,file_entry,file_format,'file_checksum','submission_date_ENA']) + '\n') | 45 runs_table.write('\t'.join([run_alias,action,'ena_run_accession',exp_alias,file_entry,file_format,'file_checksum','submission_date_ENA']) + '\n') |
46 | 46 |
47 studies_table.close() | 47 studies_table.close() |
48 samples_table.close() | 48 samples_table.close() |
49 experiments_table.close() | 49 experiments_table.close() |