Optimal minimum read length PE HiSeq reads for denovo assembly

Question: Optimal minimum read length PE HiSeq reads for denovo assembly

5.8 years ago by

Netherlands

Lesley Sitter • 550 wrote:

Hi everyone,

I have a pair ended HiSeq dataset obtained from a sort of RADseq experiment which i want to use to do a denovo assembly with. Now i have trimmed the reads using TrimGalore, but this left me with a bunch of reads of size 20 (the default minimum read length in TrimGalore).

I don't know if the reads are too short to actually use downstream. These small reads lower the optimal kmer estimations and will probably reduce quality of the assembly. I was wondering if there is some sort of method to decide the optimal minimum read length.

Thanks in advance,

With kind regards,
Lesley

read length pe reads denovo assembly hiseq • 1.9k views

ADD COMMENT • link •

modified 5.8 years ago by dylan.storey • 60 • written 5.8 years ago by Lesley Sitter • 550

5.8 years ago by

dylan.storey • 60

United States

dylan.storey • 60 wrote:

If most of your reads are that short , your assembly is going to be bad no question about it. Seeing as it is RADSeq data though I don't think you're overly interested in getting long contigs. I would suggest running your assembler with multiple parameters to test for the configuration that works best. Something like this should get you started(I'm assuming you're working on a Linux/Unix distribution with a bash like shell):

for k in {10..100}; do ./My_assembler --kmer_setting $k --output_directory_setting $k/; done;

Changing the integers in the curly brackets ({low...high}) will allow you to adjust the range.

Then you simply need to figure out how you want to assess your outputs for quality and pick the best one.

ADD COMMENT • link written 5.8 years ago by dylan.storey • 60

Thanks!

I hoped that there was some sort of tool for it, but this works just as well!

Muchos Gracias :D

ADD REPLY • link written 5.8 years ago by Lesley Sitter • 550

5.8 years ago by

Brian Bushnell ♦ 17k

Walnut Creek, USA

Brian Bushnell ♦ 17k wrote:

Reads shorter than than the kmer length used in assembly will be ignored, at least for contig-building. They could be used for scaffolding. I would throw away reads shorter than about 35bp if you are going to do scaffolding, to reduce spurious joins in low-complexity or repeat areas.

ADD COMMENT • link written 5.8 years ago by Brian Bushnell ♦ 17k

Please log in to add an answer.

Similar posts • Search »