This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Comparing contigs files and recover similar contigs

Hi, I need help with a problem. I got two files containing contigs from two different assemblers. I would like to know if it is possible to compare the contigs of go out to keep only those that are similar to 95%. I was thinking of using a BLAST, aligning the contigs of the first file on the scond but I can not extract those with 95% identity. I thank you in advance for your help and wish you a great day

genome sequence
• 1,549 views • link
updated by h.mon h.mon • written by ben.blais • ben.blais

Depending on the size of the contigs you may be able to use CD-HIT.

• 0 views • link
written by GenoMax GenoMax

Hi, thanks you, my contigs measured between 500 and 200 000 bases. Thanks you for your help

• 0 views • link
written by ben.blais • ben.blais

Sounds like you need to use an assembly reconciliation tool. This is a recent enough review that you may want to look through.

• 0 views • link
written by GenoMax GenoMax

You need global alignment for this. Try needle from the EMBOSS suite.

• 0 views • link
written by Joe Joe

What is the actual end goal here? A 'hybrid' assembly of best contigs or something?

• 0 views • link
written by Joe Joe

1 answer

You have two options:

1) You may try an assembly reconciliation or assembly merging. There are several tools for that, an internet search should help you find some. I haven't done an assembly reconciliation yet, but a recent review doesn't seem encouraging about its usefulness:

A comparative evaluation of genome assembly reconciliation tools

2) The funannotate pipeline (a fungal annotation pipeline) includes wrapper scripts to remove redundant contigs from an assembly using minimap2. It has parameters to set percent overlap and identity. You could concatenate the two assemblies and run this script - It could take a long time, though.

• 0 views • link
written by h.mon h.mon

Log in to answer this question.

More posts like this