SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
scaffold and contig chengeng Bioinformatics 5 07-04-2013 06:15 AM
Help:scaffold 2 contigs! rogerholmes.novogene Bioinformatics 11 07-28-2011 10:05 PM
SRMA Problem SAMRecord contig does not match the current reference sequence contig gavin.oliver Bioinformatics 5 07-05-2011 05:28 AM
export scaffold or link information anyone1985 Bioinformatics 2 12-19-2010 11:39 PM
Scaffold using Phusion assembler e-summer-3 Bioinformatics 0 12-05-2010 08:59 PM

Reply
 
Thread Tools
Old 08-18-2010, 11:50 PM   #1
fahmida
Member
 
Location: Australia

Join Date: Aug 2010
Posts: 54
Default How To: Contig to Scaffold?

Hi,

I would be grateful to have suggestions regarding the followings:

1) I have a set of contigs from the de novo assembly of ~2GB 454Flx data using the 'newbler' assembler (version 2.0.01). This version of 'newbler' do not create any scaffold from the contigs. I just wonder what tool I may use to make scaffolds using the contigs. I have the following output files from 'newbler': 454contigs.fna, 454contigs.ace and 454contigs.qual

2) For the same species, I have contigs from the assembler 'SoapDeNovo' using Illumina paired-end reads. I just have the 'Soap_contigs.fna' contig file. I have blasted 'newbler' contigs against 'SoapDeNovo' contigs. How I may use contigs from two different platforms to make scaffolds.

Both the assemblies are done outside as we don't have in-house HPC facility. For the task described above I just have my own desktop machine (Mac OSX 10.6, 2.26Ghz Core2Duo, 4GB RAM).

Thanks in advance for any help.

Regards.
fahmida is offline   Reply With Quote
Old 08-23-2010, 01:10 AM   #2
odysseus
Junior Member
 
Location: The Republic of Korea

Join Date: Aug 2009
Posts: 8
Default

Hi~

1) Are contigs created from paired-end reads? If they are, you should assemble them using paired-end assembly function of Newbler. After paired-end assembly, you can get a set of scaffold.

2) In my think, if you have raw reads, you can hybrid assemble them using hybrid assembler such as MIRA3. If you solve this difficulty, please inform me about the solution.

Cheers! ^^
odysseus is offline   Reply With Quote
Old 08-23-2010, 01:49 AM   #3
francesco.vezzi
Member
 
Location: Udine (Italy)

Join Date: Jan 2009
Posts: 50
Default

MIRA3 is probably what you need...
Remember that SOAPdenovo gives as output a file named *.scafSeq that contains the scaffolds (from my experience the contig file contains a lot of very short and repetitive contigs often shorter than the read length!!!)

Another solution is to give the output of the two assemblers as input for a "standard" assembler like PCAP or ARACHNE...

Anyway I think that you resources are really too low...
francesco.vezzi is offline   Reply With Quote
Old 08-22-2013, 04:41 AM   #4
hugh_hang
Member
 
Location: Hangzhou, China

Join Date: Jan 2013
Posts: 28
Default

Quote:
Originally Posted by francesco.vezzi View Post
MIRA3 is probably what you need...
Remember that SOAPdenovo gives as output a file named *.scafSeq that contains the scaffolds (from my experience the contig file contains a lot of very short and repetitive contigs often shorter than the read length!!!)

Another solution is to give the output of the two assemblers as input for a "standard" assembler like PCAP or ARACHNE...

Anyway I think that you resources are really too low...
if my data is 454 sequencing data from cDNA, is it necessary to generate scaffolds? OR can the scaffolds make any sense?
hugh_hang is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:28 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO