Maybe add a check here that all contig ids were found since there should be a one to one mapping between the fasta file contig ids and the cluster output contig ids? Just to make sure all contigs are in a cluster output file. So make sure len(seqs) == len(contig_ids) ??