I have assembled a transcriptome using TRINITY with 100,000 transcripts. I have perform blastx against NR database (-outfmt 5). I want to remove redundancy in the assembled transcript before proceed to further processing. 

1) How can I cluster them together into unigenes and remove redundancy? 

2) Or how to select transcipts with longest sequences if they return same hit? 

3) Any software or program suggestion for doing this?

4) Is it necessary for removing redundancy?

Thank you very much.




You will get transcripts after assembling short reads. Depending on alternative splicing,one gene can have more than one transcript. 

Clustering or grouping puts all these similar sequences together and help to make a set of transcript for one gene. 

On the other hand assembly process can create some transcripts that are not real (sequences/transcripts with more than 95% identity to another sequence) and clustering helps identifying them. Therefore, the flow goes like:

Assembly -> Transcripts -> clustering transcripts -> unigenes -> further downstream processing.

You can use [CD-HIT][1] for this purpose.


  [1]: http://weizhongli-lab.org/cd-hit/