diff --git a/workflow/scripts/pool_and_psuedoreplicate.py b/workflow/scripts/pool_and_psuedoreplicate.py index a0a6a988a0199c7a49257f4b00b32b1dfcd0cce0..cdf4d5b40a060094d97f322cbce4f4ee8da9626b 100644 --- a/workflow/scripts/pool_and_psuedoreplicate.py +++ b/workflow/scripts/pool_and_psuedoreplicate.py @@ -230,17 +230,19 @@ def generate_design(paired, cutoff_ratio, design_df, cwd, no_reps, no_unique_con design_new_df.at[3, 'tag_align'] = design_new_df.at[0, 'tag_align'] # Make 2 self psuedoreplicates - pool_pseudoreplicates_dict = {} - for rep, tag_file in zip(design_df['replicate'], design_df['tag_align']): + self_pseudoreplicates_dict = {} + for rep, tag_file in zip(design_df['replicate'], design_df['tag_align']): replicate_prefix = experiment_id + '_' + str(rep) - pool_pseudoreplicates_dict[rep] = \ - self_psuedoreplication(tag_file, replicate_prefix, paired) + self_pseudoreplicates_dict = self_psuedoreplication(tag_file, replicate_prefix, paired) + # Update design to include new self pseudo replicates - for rep, pseudorep_file in pool_pseudoreplicates_dict.items(): + pool_pseudoreplicates_dict = {} + for rep, pseudorep_file in self_pseudoreplicates_dict.items(): path_to_file = cwd + '/' + pseudorep_file replicate = rep + 1 design_new_df.loc[replicate, 'tag_align'] = path_to_file + pool_pseudoreplicates_dict[replicate] = path_to_file # Drop index column design_new_df.drop(labels='index', axis=1, inplace=True)