SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to compare two WIG files?? tfcheng Bioinformatics 1 08-28-2013 01:23 AM
Compare two bam files. fabrice Bioinformatics 4 10-03-2011 03:20 PM
compare samples jgSoton Bioinformatics 2 10-03-2011 12:38 AM
best way to compare two (or more) genoms NicoBxl Bioinformatics 1 03-24-2011 06:16 AM
Compare SNPs bair Bioinformatics 5 01-27-2010 02:38 AM

Reply
 
Thread Tools
Old 06-28-2011, 10:50 AM   #1
infrared1983
Junior Member
 
Location: DALLAS

Join Date: Apr 2011
Posts: 1
Default Compare Genomes

Hi everyone, I did Genome assembly by SOAPdenovo and obtained contigs and scaffolds, now I wanna compare these scaffolds to some reference genomes downloaded from Genbank. I am wondering what tools I can use to do this job, I used BLAST but it is not in good visualization, any other advice on such tools and programs?
I am new to this, any suggestion is appreciated.
Thanks
infrared1983 is offline   Reply With Quote
Old 07-03-2011, 08:36 PM   #2
jmw86069
Member
 
Location: RTP, NC, USA

Join Date: Jun 2009
Posts: 28
Default

Use NUCmer from the MUMmer package (available from Sourceforge.) For large reference genomes, you have to split the reference into smaller pieces for it to fit in memory, then align everything to those pieces. It is relatively fast and efficient, and seems to be a commonly used tool as the basis for other assembly comparison pipelines. There are some recently available wrapper scripts which automate some of the steps, but it's not tough to run through the steps yourself. The online docs for MUMmer as a whole are very good.
jmw86069 is offline   Reply With Quote
Old 07-04-2011, 05:34 PM   #3
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

For mammalian genomes, use bwa-sw. Two human assemblies in several CPU days.
lh3 is offline   Reply With Quote
Old 08-05-2011, 06:42 PM   #4
HerefordGuy
Junior Member
 
Location: Missouri

Join Date: Apr 2010
Posts: 6
Default De novo to reference comparison

I have had good luck with both MUMmer and MAUVE.
MAUVE is very straight forward and MUMmer I simply follow the recipes in the manuel.
HerefordGuy is offline   Reply With Quote
Old 10-12-2011, 11:29 AM   #5
narain
Member
 
Location: Washington DC

Join Date: Aug 2011
Posts: 78
Default

But is it good to use the contigs files or to use the scaffolds files for comparison?
narain is offline   Reply With Quote
Old 10-12-2011, 12:03 PM   #6
HerefordGuy
Junior Member
 
Location: Missouri

Join Date: Apr 2010
Posts: 6
Default Compare Genomes

I suppose it depends what your question is.
But, typically the goal of genome assembly is to create one scaffold per chromosome (assemble the entire genome/chromosome). With this in mind, I think it would make more sense to compare scaffolds to scaffolds.
HerefordGuy is offline   Reply With Quote
Old 10-12-2011, 12:29 PM   #7
narain
Member
 
Location: Washington DC

Join Date: Aug 2011
Posts: 78
Default

However, scaffolds are nothing but a combination of contigs filled with "estimates" of gaps in between. Does that not mean then that it is more sensible to compare contigs than to compare scaffolds as it has estimate of gaps in between which may be highly erroneous!
narain is offline   Reply With Quote
Reply

Tags
soap denovo

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 06:29 PM.


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