Helo everyone, After blastx Unigenes.fa from de novo RNA-Seq assemblies to a reference proteome, I obtained the result that multicontigs belong to one gene. So I just try to uniform my data that to remove shorter congtigs and keep the longest contig representing one gene. example is as below:
From origal data
ENSDARP00000143232.1 GGCTCCTCTTTTTCAACTGGACATCCTTAAAACTGTATGAAAGGGGCGGAGCCTTTTGCTACTTGCATACTTAAGCTCCTTCACATTCCTCTAGCCCTTTACGAA ENSDARP00000143232.1 GGCTCCTCTTTTTCAACTGGACATCCTTAAAACTGTATGAAAGGGGCGGAGCCTTTTGCTACTTGCATACTTAAGCTCCTTCAC ENSDARP00000143232.1 GGCTCCTCTTTTTCAACTGGACATCCTTAAAACTGTATGAAAGGGGCGGAGCCTTTTGC
To what I want
ENSDARP00000143232.1 GGCTCCTCTTTTTCAACTGGACATCCTTAAAACTGTATGAAAGGGGCGGAGCCTTTTGCTACTTGCATACTTAAGCTCCTTCACATTCCTCTAGCCCTTTACGAA
Could you give me some suggestions or some scrips to help me, thanks!
0 answers
No answers yet.
Log in to answer this question.
-
Repeat Modeler Database Generation
written by SomeOne •
-
Repeat Modeler Database Generation
written by SomeOne •
-
RNA de novo assembly - blasts - KEGG - GO
written by dimitrischat
-
Trinity assembly way too many transcripts!
written by karsa692 •
-
Transcriptome blasting doubts
written by pablo61991 •
-
extracting longest unigene from de novo RNA-Seq assemblies
written by bioming •
-
Trinity, longest contigs
written by loly.pearl86 •
-
How efficiently to use information of two de novo assembled transcriptome for further analysis?
written by seta
-
How efficiently to use information of two de novo assembled transcriptome for further analysis?
written by seta
-
Optimizing De Novo Transcriptome assemblies for non model organisms
written by giorgiocasaburi
Have u tried
cdhit. Use query length to filter from blast file. How did you obtain unigenes?