Compare Genomes

infrared1983 · 06-28-2011, 11:50 AM

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

jmw86069 · 07-03-2011, 09:36 PM

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.

lh3 · 07-04-2011, 06:34 PM

For mammalian genomes, use bwa-sw. Two human assemblies in several CPU days.

HerefordGuy · 08-05-2011, 07:42 PM

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.

narain · 10-12-2011, 12:29 PM

But is it good to use the contigs files or to use the scaffolds files for comparison?

HerefordGuy · 10-12-2011, 01:03 PM

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.

narain · 10-12-2011, 01:29 PM

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!

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

06-28-2011, 11:50 AM	#1
infrared1983 Junior Member Location: DALLAS Join Date: Apr 2011 Posts: 1	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

08-05-2011, 07:42 PM	#4
HerefordGuy Junior Member Location: Missouri Join Date: Apr 2010 Posts: 6	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.

10-12-2011, 01:03 PM	#6
HerefordGuy Junior Member Location: Missouri Join Date: Apr 2010 Posts: 6	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.

07-03-2011, 09:36 PM	#2
jmw86069 Member Location: RTP, NC, USA Join Date: Jun 2009 Posts: 28	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.

07-04-2011, 06:34 PM	#3
lh3 Senior Member Location: Boston Join Date: Feb 2008 Posts: 693	For mammalian genomes, use bwa-sw. Two human assemblies in several CPU days.

10-12-2011, 12:29 PM	#5
narain Member Location: Washington DC Join Date: Aug 2011 Posts: 78	But is it good to use the contigs files or to use the scaffolds files for comparison?

10-12-2011, 01:29 PM	#7
narain Member Location: Washington DC Join Date: Aug 2011 Posts: 78	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!