Hi friends!
This is more of a feature request than an issue.
Have you considered implementing fastq generation from a reference genome? I think it would be a useful feature with a lot of extensibility. It might be out of scope, but I think there is a lot someone could do with this (I can imagine bed / vcf support, etc).