i need to generate some artifical reads.
randomsreads.sh from BBMap seems to be the right program for me.
I don't need random reads for my whole reference genome but for several thousand regions. So I wonder if there is better/faster way instead of first creating a separate ref file for every region? Maybe a new parameter for randomreads.sh for providing a bed file with the region (@Brian Bushnell :) )