reorder contigs given by nucmer
2
1
Entering edit mode
7.7 years ago
Phil S. ▴ 690

Hi there,

I used nucmer to place (de novo assembled) contigs onto a close relative and as far as I can tell from the dotplot it turned out really really well. Now I want to reorder the contigs into the ordering found by nucmer. Is there any tool or script out there which is able to do that?

Thanks,
phil

assembly nucmer guided genome • 4.1k views
ADD COMMENT
3
Entering edit mode
7.7 years ago
SES 8.5k

There is a tool called ABACUS that is designed for that exact purpose, and in my experience it works really well. This program actually uses nucmer to align the contigs to the reference, so this sounds like a good fit. ABACUS is also a part of the larger PAGIT pipeline for assembly improvement.

ADD COMMENT
0
Entering edit mode

Hi,

thanks for the fast reply. Unfortunatelly it seems like the abacas does not get quite the results using nucmer alone. I dont know why but using nucmer alone I end up getting something like this:

Using abacas I get just around half of the bases placed and the rest is filled up with Ns. I can't see any switch for the amount of Ns inserted or when they are inserted etc...

edit: I just saw that MAUVE also has a 'move contig' option which is able to do that. Will keep you posted if that worked better. Anyways if someone has a different solution I'm more than happy to hear about it!

Cheers

ADD REPLY
0
Entering edit mode

Have you tried modifying the options to ABACUS? I would check your settings for running mummer alone, and then give those settings to ABACUS. You can change the percent identity, coverage, etc. of matches, so you should be able to get close to the same results. There will likely be differences because of the post-processing of the matches.

ADD REPLY
0
Entering edit mode

I have a look into that right now. Thanks!

ADD REPLY
0
Entering edit mode

@SES But ABACAS does not stitch the tips which are overlapping. It just reorders and places N's in between contigs. In that case, there will be same region twice in a scaffold. Isn't it?

ADD REPLY
0
Entering edit mode

Yep it seems like ABACAS does not stitch. I modified ABACAS that it nearly produced the same results. Though I am figuring out how to stitch manually....(meanwhile testing the 'move contigs' function MAUVE provides..)

Thanks for all your thoughts / ideas

ADD REPLY
1
Entering edit mode
6.0 years ago
Andreas ★ 2.5k

Hi,

I'm a bit late to the party but:

Mummer's show-tiling does the trick. It can also generate a concatenated pseudo-sequence (gaps filled with Ns) if you use option -p

Andreas

ADD COMMENT

Login before adding your answer.

Traffic: 1476 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6