MetaGenome -> Reads

[HadrienG]

Issue to track the progress on the Roadmap item "From a set of genomes, output perfect reads in fastq format given abundance parameters"

Test files:

• genomes.fasta
• abundance.txt

Number of genomes: 5

• Parse the file and get one SeqRecord per genome (we still assume complete genomes: one sequence per genome)
• Parse abundance.txt and calculate the necessary coverage per genome to satisfy the abundance requirements
• Generate the reads with the generator.reads() function
• Use multiprocessing to parallelise the process (1 threads per genome?)
• Write the reads to file
• write tests

[HadrienG]

Commit 3849f00 closes this, minus the multiprocessing support. Multithreading will have its own issue one the simulator will be up and running single-threaded.