comparison extract_tables.py @ 13:91a3a617d94b draft

Uploaded
author ieguinoa
date Fri, 18 Sep 2020 14:29:47 +0000
parents b85d239b1d58
children 436da36c6114
comparison
equal deleted inserted replaced
12:6694e2543e4b 13:91a3a617d94b
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'],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'].lower(),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.gz' if os.path.splitext(file_entry)[-1] == '.gz' else 'fastq.bz2' 44 file_format = 'fastq'
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()