When simulating thousands of genomes with multiple contigs (MAGs), issues were encountered (see #54 and #53) as processing large numbers of contigs slows down the pipeline.

We could improve performance by processing fastas in a single rule and stop splitting contigs.