Yeah, this script is sort of made in a haste so it probably needs proof reading. If the contig id is not in the original fasta file the script will crash. If the a contig is in the fasta file but not in the clustering it will be missing in the output files. That should be the case since we are not clustering all contigs but only the ones greater than the threshold.