Question: Splitting Abacas Single Fasta Output
3
gravatar for Darked89
7.2 years ago by
Darked894.2k
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 • 1.9k views
ADD COMMENTlink modified 7.0 years ago by Lee Katz2.9k • written 7.2 years ago by Darked894.2k
2
gravatar for Lee Katz
7.2 years ago by
Lee Katz2.9k
Atlanta, GA
Lee Katz2.9k 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 7.2 years ago by Lee Katz2.9k
1

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 7.0 years ago by mgalactus720
1

Yeah! I'm a huge fan of CONTIGuator!

ADD REPLYlink written 7.0 years ago by Lee Katz2.9k
Please log in to add an answer.

Help
Access

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