SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Comparing and merging genome assemblies megh Bioinformatics 5 07-23-2014 03:58 AM
Finding *new* regions of DNA in genome assemblies green tree De novo discovery 5 02-20-2012 02:19 PM
How to align SNP coordinates against a circular reference genome? Lspoor General 1 02-17-2012 10:44 AM
PubMed: [Letter to the editor] Whole-genome amplification of sodium bisulfite-convert Newsbot! Literature Watch 0 04-14-2011 07:00 AM
PubMed: Integrating genome assemblies with MAIA. Newsbot! Literature Watch 0 09-09-2010 02:00 AM

Reply
 
Thread Tools
Old 05-10-2011, 12:30 PM   #1
TuA
Junior Member
 
Location: CH

Join Date: May 2011
Posts: 9
Post How to convert genome coordinates from two assemblies

Hi!

Iím working with two closely related species ( < 1 MY) and got a reference genome for both of them. I re-sequenced several individuals per species and mapped the reads to the reference genome of the particular species. Now, I would like to look at SNPs between and within the species. For this, I obviously need to have the same genome coordinates for the SNP positions. Does anyone have a good advice how the convert the coordinates? My main concerns are how to deal with indels and how I could still be able to use additional available data of one of the species such as genome annotation.

Your help is greatly appreciated!
TuA is offline   Reply With Quote
Old 05-10-2011, 11:16 PM   #2
simonandrews
Simon Andrews
 
Location: Babraham Inst, Cambridge, UK

Join Date: May 2009
Posts: 871
Default

You should look at the liftOver tool from UCSC which can be used to covert between different coordinate systems. You can create a set of 'chain' files which allow you to convert between coordinates in your different assemblies. The matching is based on either blat (if the species are close enough for a DNA level comparison) or blastz for more distant relationships.
simonandrews is offline   Reply With Quote
Old 05-16-2011, 03:57 AM   #3
TuA
Junior Member
 
Location: CH

Join Date: May 2011
Posts: 9
Default

Thanks Simon for your helpful post. It seems to be pretty time intense to create these liftOver chain files. Do you have any experience on the accuracy, i.e. how well the coordinates can be transformed (blat performance)?

Many thanks
TuA is offline   Reply With Quote
Old 05-16-2011, 05:03 AM   #4
simonandrews
Simon Andrews
 
Location: Babraham Inst, Cambridge, UK

Join Date: May 2009
Posts: 871
Default

Quote:
Originally Posted by TuA View Post
Thanks Simon for your helpful post. It seems to be pretty time intense to create these liftOver chain files. Do you have any experience on the accuracy, i.e. how well the coordinates can be transformed (blat performance)?
The liftOver chains might take a long time to calculate, but once you have them the speed with which coordinates can be transformed is blazingly quick. I guess the efficiency of the conversion will depend on the degree of identity between your genomes. We mostly use the system for converting between different assemblies of the same genome, and it's very accurate and quick there. If your two genomes are reasonably high identity then you should have no problem.
simonandrews is offline   Reply With Quote
Old 05-23-2011, 01:17 AM   #5
TuA
Junior Member
 
Location: CH

Join Date: May 2011
Posts: 9
Default

Thanks for your help! I'll give it a try...
TuA is offline   Reply With Quote
Old 11-20-2013, 07:43 PM   #6
liguow
Member
 
Location: Houston

Join Date: Apr 2009
Posts: 12
Default CrossMap

CrossMap is a program for convenient conversion of genome coordinates between assemblies. It supports most commonly used file formats including SAM/BAM, Wiggle/BigWig, BED, GFF/GTF, VCF.

http://crossmap.sourceforge.net/
liguow is offline   Reply With Quote
Reply

Tags
assembly, coordinates

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


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