SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   General (http://seqanswers.com/forums/forumdisplay.php?f=16)
-   -   How to align SNP coordinates against a circular reference genome? (http://seqanswers.com/forums/showthread.php?t=8939)

Lspoor 01-18-2011 07:33 AM

How to align SNP coordinates against a circular reference genome?
 
I have a list (tab-delimited file) of approximately 18 000 SNP coordinates for a bacterial reference genome GenBank file (When compared to one other query FASTA genome sequence using "dnadiff" function of Mummer3).

I am looking to produce a circular representation of this reference genome, with the SNP coordinates mapped against it, to look for hotspots that might indicate recombination events may have occurred. I basically want to look at the SNP distribution throughout the genome.

For smaller numbers of SNPS (1000 to 2000), I was using an open source software called BRIG but am unable to get this to work with the much larger SNP numbers?

Does anyone have any suggestions for another (vaguely user-friendly way!) I could do this, as am a relative Newbie to Linux but also have access to Windows. I have produced an OK-ish diamgram using GenomeDiagram with Biopython but it is not easy to maniupulate the diagram to get a better look at the SNP distribution )eg change the size or width of the features representing the SNPS).

Any suggestions would be very much appreciated! :confused:

houkto 02-17-2012 10:44 AM

Hi there,

I am trying to achieve the same thing. I believe almost everyone is going for http://circos.ca/ you can install it in windows or linux and it very well documented and they have a forum in which you can troubleshoot your problems. if you find something else please let me know.

Cheers


All times are GMT -8. The time now is 11:10 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.