Hi,
I am working with single cell amplified genomes (SAGs), which leads to a very uneven distribution of reads among the genome; some regions will have extremely high coverage and some extremely low. To make better assemblies, it is required that the high coverage are 'downsampled', i.e. discarding reads in a high coverage region until that region reaches a certain coverage. This would then ultimately lead to a more uniform coverage, easing the assembly.
Is there some way or script to do this? I use velvet as assembler.
I am working with single cell amplified genomes (SAGs), which leads to a very uneven distribution of reads among the genome; some regions will have extremely high coverage and some extremely low. To make better assemblies, it is required that the high coverage are 'downsampled', i.e. discarding reads in a high coverage region until that region reaches a certain coverage. This would then ultimately lead to a more uniform coverage, easing the assembly.
Is there some way or script to do this? I use velvet as assembler.