downsampling high coverage regions in SAGs

novi · 06-30-2011, 02:28 AM

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.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
When viewing a BAM, can you specify regions of a certain coverage?	oiiio	Bioinformatics	3	01-10-2012 07:59 AM
Agilent SureSelect - coverage of high GC regions	NGSfan	Genomic Resequencing	28	11-07-2011 02:13 PM
"coverage" of introns, intergenic regions for RNASEQ	PFS	Bioinformatics	2	09-07-2011 02:51 PM
low 454 coverage combined with high solexa coverage	strob	Bioinformatics	7	10-07-2010 11:14 AM
PubMed: High-throughput sequencing of microdissected chromosomal regions.	Newsbot!	Literature Watch	0	11-06-2009 03:00 AM

06-30-2011, 02:28 AM	#1
novi Junior Member Location: Uppsala Join Date: Jun 2011 Posts: 2	downsampling high coverage regions in SAGs 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.