# Write cluster to a file fname = "clustersize{0}-No{1}.fastq".format(str(cluster.size), str(size_counter[str(cluster.size)])) output_handle = open(fname, "wb") SeqIO.write(seqs, output_handle, "fastq") elif self.reads_per_cluster_size_counter[str(cluster.size)] < self.filter_params['min_reads']: break # Sort file sortby(args.filename, reverse=True, mode='reads_per_cluster', outfile_postfix='-filtered', cutoff=args.minreads, clustersize_min=args.min, clustersize_max=args.max) # Rewrite file old_filename_parts = args.filename.name.split('.') old_filename_parts[0] += '-sortedby_' + 'reads_per_cluster' new_filename = '.'.join(old_filename_parts)