SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
alignment questions...pooled-seq, multiple references, pseudogenes, plastid genomes jullee Bioinformatics 3 04-16-2014 12:18 AM
vcf compare LilianK Bioinformatics 0 12-24-2013 04:49 AM
How compare together 3 partial genomes? mm.perrineau De novo discovery 3 12-11-2012 02:53 PM
Compare Genomes infrared1983 Bioinformatics 6 10-12-2011 01:29 PM
best way to compare two (or more) genoms NicoBxl Bioinformatics 1 03-24-2011 07:16 AM

Reply
 
Thread Tools
Old 01-07-2016, 09:48 AM   #1
ootunaoo
Junior Member
 
Location: US

Join Date: Jan 2016
Posts: 1
Default Questions regarding how to compare 2 genomes

Hi Seqanswers,

I would like to ask several questions regarding how to compare 2 genomes in order to find differences: Assume I have 2 dataset of sequencing data from 2 plants of a same species (e.g. arabidopsis) - 1 plant has normal phenotype, the other has disease phenotype. Theoretically, the disease phenotype is known to be controlled by a single gene, and these 2 plants should have similar genome accept the region that responsible for the phenotype. I would want to somehow compare the 2 genomes to find out differences between 2 plants (in order to find the disease gene).

I'm a newbie in Bioinformatics (also newbie of Seqanswers), I do not know where to start. Would you mind providing me some guides to help me find an approach for my problem? - Projects or publications that have similar object; Documents, internet or books, that I should read; Maybe a suggested pipeline would be great...

Thank you for reading.

P/S: And I am sorry if my english is horrible.
ootunaoo is offline   Reply With Quote
Old 01-07-2016, 09:58 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,989
Default

The following tool may or may not work in your case but it may be worth taking a look: http://genetics.bwh.harvard.edu/snptrack/ There is a paper linked at the site as well.
GenoMax is offline   Reply With Quote
Old 01-10-2016, 12:23 AM   #3
gsgs
Senior Member
 
Location: germany

Join Date: Oct 2009
Posts: 140
Default

mathematically speaking,

suppose you have the two mappings a:{1,..,n}-->{A,C,G,T} and b:{1,..,m}-->{A,C,G,T}
representing the two genomes.

pick L (e.g. L=16) and compute
f:{1,..,n-L}-->{0,1}
with
f(x)=1, iff exists y such that a(x+i)=b(y+i) , i=0..L-1

this can quickly be computed by marking all values of b in a 4^L table.
[you may add the inverse complement of b() here]

then plot moving averages of f, the number of values in the averages being
approximately the length of the expected gene.

this gives an overview of the matching-quality by genome-region

you should see a "valley" in a nonmatching region


[is there a name for this function ?]
gsgs 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 08:12 AM.


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