Entering edit mode
6.7 years ago
alesssia
▴
580
Hello everyone,
I need to simulate a few metagenomics samples and I was thinking to use BBMap's simulator, randomreads.sh, in "metagenome" mode.
The problem is that I need to simulate the effect of PCR amplification, and I thought of including identical duplicates but with different quality scores.
Does anyone have any suggestion on how to do this? Or knows if BBMap randomreads.sh already simulates PCR duplicates (and if yes, how do I modify their proportion)?
Thanks a lot!
It probably does not simulate PCR dups. You could go the route you mention above.
Thanks, I think I will just generate the metagenomics samples, then randomly select varying percentage of reads, duplicate them, and finally shuffle their quality score. Will this work in your opinion?
You could check if there are PCR dups after generating the metagenome (using A: Introducing Clumpify: Create 30% Smaller, Faster Gzipped Fastq Files ). If not then proceed with your plan as stated above.
Good point. Since I need to keep the percentage of duplicates equal to fixed values I should first check if there are any and add/remove if needed! Thanks!