Question: Splitting Abacas Single Fasta Output
gravatar for Darked89
8.3 years ago by
Barcelona, Spain
Darked894.2k wrote:

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.1k views
ADD COMMENTlink modified 8.0 years ago by Lee Katz3.0k • written 8.3 years ago by Darked894.2k
gravatar for Lee Katz
8.2 years ago by
Lee Katz3.0k
Atlanta, GA
Lee Katz3.0k wrote:

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 COMMENTlink written 8.2 years ago by Lee Katz3.0k

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

ADD REPLYlink written 8.1 years ago by mgalactus730

Yeah! I'm a huge fan of CONTIGuator!

ADD REPLYlink written 8.1 years ago by Lee Katz3.0k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2056 users visited in the last hour