This is a test version of Biostars. For the public version, visit https://www.biostars.org.
why cd-hit-est tool could not recognize identical contigs in the assembled transcriptome input file?

Hi everybody,

I have almost high quality Illumina reads that were trimmed and assembled using CLC genomic workbench software, I tried different K-mer size and got the best assembly, in terms of some basic parameters, like N50, the number of contigs and the percentage of mapped back reads (about 90%) in K-mer of 64. This assembly (44.8 MB in size) was subjected to cd-hit tool with strict 100% identity for the alignments since plants have many paralogues with have high sequence identity, but, the size of output file has not been changed. I also tried this analysis with stringency of 0.9, and the output size file (43.5MB) did not change significantly as compared with the original input file (44.8 MB). Could anybody please let me know whether these results are usual or there is something wrong? Thanks for any comments

rna-seq assembly blast
• 2,736 views • link
updated by Ram Ram • written by seta seta

Why would you expect to have identical contigs in your assembly? The point of an assembler is to collapse and extend overlapping sequences.

• 0 views • link
updated by Ram Ram • written by Istvan Albert Istvan Albert

1 answer

I have no experience in this filed and expect some redundancy just based on similar published works. So, it's normal in your view

• 0 views • link
updated by Ram Ram • written by seta seta

I don't think any assembler should ever generate exactly identical contigs. I can understand how similar contigs could be present but it would not be reasonable to expect identical contigs.

• 0 views • link
written by Istvan Albert Istvan Albert

Log in to answer this question.

More posts like this