This is a test version of Biostars. For the public version, visit https://www.biostars.org.
merge PE reads for Transcript De-Novo Assembly ??

Dear all,

Do you recommend merging pair-end Illumina reads for Transcript De-Novo assembly? Thank you for any advice and sharing your experiences..

rna-seq transcript assembly de-novo

Yes, I got better results after joining the reads, but some bioinformaticians would be probably disagree. So , better is to go for de novo assembly by both ways, then check the average contig length, total number of contigs and length of assembly and then you can argue that which one worked better for you

hth

Thank you for comment - probably it would be the best solution...

2 answers

It all depends on how many reads overlap If the majority of the reads overlap then there is little that should be gained from treating them as paired end. In fact it should be counterproductive to do so as the system has to deal with more and redundant data.

Logic dictates that providing more information to the system ought to make it perform better. In this case the extra information is that the reads are overlapping and an external tool solved that problem.

Now in reality and practice, the way algorithms are built, tuned and released, depending on the tool and version it just might be that you end up with unexpected performance when choosing one option vs the other. Hence as Manvendra Singh suggests I think it is best to be cautious and evaluate both methods.

My PE overlap 90%, I merged them and got really good results with high k-mers. I only ued merged PEs in the assembly, avoiding many chimeric scaffolds.

Thank you guys for the comment and sharing your experiences.. What tools is the best to evaluate if my reads are overlap?

This really depends on the assembler (as well as the merging program!). I have found that merging reads improves Ray assemblies, and makes Soap assemblies dramatically worse.

Log in to answer this question.