Hello fellow bioinformaticians!

I was wondering if anyone had any clue as to whether or not a tool exists for generating dummy/fake reads from an already assembled genome? So far my searches have yielded very little...

Essentially, the ideal software would take a single-contig assembled genome and generate simulated raw reads based on some input criteria such as desired average coverage, level of base quality, etc... .

People are not always nice enough to release the raw reads from their sequencing projects, hence simulated raw reads is the second best thing.

I realize this might be a very difficult task, but assuming no tool exists for this, do you think it would be possible to make a decent one?

Thank you for the input!
[BBMap][1]'s 'randomreads' does exactly what you need.


  [1]: https://sourceforge.net/projects/bbmap/