Beispiel #1
0
        
    # 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)