Hi all,
I have a weird question! I have some contiges and want to convert them to all possible illumina reads (suppose with length 100). I am wondering if there is any tool to do this?
Actually my problem is vise versa of the assembly problem. In assembly problem we have the reads and want to find the contigs. But here I have the contiges and want to find all possible small reads that can produce this contig.
Thank you,
While there are many tools out there, you can use randomreads.sh from BBMap suite with your contig. See the in-line program help for usage information.
Additionally, for the specific problem where you want all possible reads, try:
shred.sh in=contigs.fa out=reads.fq length=100 minlength=100 overlap=99
Log in to answer this question.
