Splitting Abacas Single Fasta Output
1
3
Entering edit mode
9.6 years ago
Darked89 4.2k

I used ABACAS to order contigs from my draft assembly (DRAFT). The output fasta file has all 30+ chromosomes glued as one sequence with contigs separated by Ns. This to be expected, since ABACAS demands single fasta reference sequence (REF), so it has no idea about chromosomal ends. While such pseudomolecule is still usable, I want to chop it into individual chromosomes/mitochondrion.

One idea is to blast both chromosomal ends from REF against DRAFT, but if DRAFT's chromosomal ends stick out it will get more complicated than simply taking blat hits ends.

Other is to start with doctored REF sequence where individual fasta headers are replaced by i.e. 11111 Ns to make sure no gaps of this size will be present, rerun ABACAS, then chop at gaps with 11111Ns.

Any other ideas, preferably already tested?

EDIT: spell

fasta split • 2.3k views
ADD COMMENT
2
Entering edit mode
9.6 years ago
Lee Katz ★ 3.1k

CONTIGuator handles this all really well and is a wrapper around abacas. I think if you can, just re-run your comparison with CONTIGuator, and it will have one result per reference chromosome.

ADD COMMENT
1
Entering edit mode

To be fair, the current version of CONTIGuator wraps abacas just for the primer picking part, but thanks for having mentioned it :)

ADD REPLY
1
Entering edit mode

Yeah! I'm a huge fan of CONTIGuator!

ADD REPLY

Login before adding your answer.

Traffic: 1617 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