Merging Multiple De Novo Assembly Fasta Files
1
1
Entering edit mode
11.0 years ago
rob234king ▴ 610

I have a series of assemblies in fasta format from SOAP2 using different k values. I would like to merge the resulting assemblies to one file because some k values assemble parts that others don't, to create hopefully larger contigs. I thought I might be able to use oases to merge the files as this is what happens for the velvet/oases transcriptome assembly but not sure how I could do that. The commands don't appear to lend themselves to this application. Any ideas would be welcomed. I also have another three re-sequencing assemblies of different k values which could also try merging to this one re-sequencing data set.

A reference exists although as ~30 contigs, some are very large so it's quite good (there should be only four chromosomes in the genome).

I don't want to manually try to do this, I have already manually plugged some large N gaps in the first reference contig but too time costly and have not looked to stick contigs together yet either.

assembly • 6.9k views
ADD COMMENT
0
Entering edit mode

you can have a look https://github.com/ghyao/GAA

ADD REPLY
0
Entering edit mode

thanks, I've been looking around and looks like the most user friendly and maybe use it to merge one good assembly with the reference and hopefully it will at least replace the N's in the reference sequence.

ADD REPLY
0
Entering edit mode

Please let me know the performance of GAA

ADD REPLY
0
Entering edit mode

I cant understand why this program needs contigs to be named in a certain fashion. The contigs from 2 different assemblies are completely unrelated in order or naming.

ADD REPLY
0
Entering edit mode
11.0 years ago

Minimus2 ?

ADD COMMENT
0
Entering edit mode

Looks like may be useful with multiple rounds of merging two files but looks like would have to install the whole amos package. I'll give the installation a try and see if works. Thanks

ADD REPLY
0
Entering edit mode

Had to convert the configure file to linux and then ran but got below error? any ideas? error: cannot run /bin/sh config/config.sub

ADD REPLY

Login before adding your answer.

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