Assembling long similar contigs
I have a program which outputs ~50 assembled transcripts which are about 10k base pairs long each.
My program already filters exact duplicate sequences but many of the assembled transcripts are very similar to one another.
Is there any existing assembly program which can connect sequences that are identical over 90% of their length?
2 answers
Dedupe from the BBMap package can remove similar sequences to leave only a single copy.
dedupe.sh in=transcripts.fa out=deduped.fa minidentity=0.9 maxedits=20
• 0 views
•
link
written
by
Brian Bushnell
Log in to answer this question.
More posts like this
-
Organelle genome assembly did not circularize and I am not sure why
written by Арсений •
-
PCR-duplicates and RAD-sequencing
written by Roland •
-
About 10x chromium Gene assembly
written by r00628112 •
-
How to get list of different bases and their positions from an alignment?
written by SaltedPork
-
RNA-seq long reads vs. transcripts
written by tunl
-
Understanding de-novo transcriptome assembly quality using Transrate
written by nayaraevelin2 •
-
Mapping of De novo assembled transcripts to draft genome
written by Bioinfonext
-
Transcriptome assembly: Low GC, short contigs, low read alignment
written by jmah •
-
Newbler, Parameters To Improve My Cdna Assembly De Novo
written by aleix.arnau1990 •
-
Assembly Error In Tandem Duplicated Region
written by cacaucenturion •
Have you tried CD-HIT? It can be used for clustering and comparing protein or nucleotide sequences.