SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
comparing large sets of sequences gsgs Bioinformatics 11 01-10-2013 04:28 AM
Comparing trees joxcargator73 Bioinformatics 1 12-12-2012 01:25 PM
comparing alignments sooper Bioinformatics 1 10-11-2012 04:46 AM
Comparing Transcriptomes poregide Bioinformatics 1 08-25-2010 11:23 PM
comparing sequences and mapping annotations? gaster Bioinformatics 0 07-06-2009 01:33 PM

Reply
 
Thread Tools
Old 02-15-2013, 08:48 AM   #1
ramirob
Member
 
Location: Vermont

Join Date: Apr 2012
Posts: 14
Default Comparing sequences between themselves

Hello,

We ran an experiment where we sequenced WT and mutant isolates from the same particular organism. We have aligned these samples to reference (using BWA) and then found SNPs (using GATK). However, aside from identifying SNPs by comparing to reference, we would like to compare the aligned samples between themselves to look for important differences (eg. WT to Mutant samples). Does anyone have any suggestions on how to approach that? It can certainly be done indirectly by comparing the SNP calls, but I was wondering if there is any existing software that does something like that (e.g. find SNPs between aligned samples as opposed to find SNPs when compared to reference).

Thanks in advance,
Ramiro
ramirob is offline   Reply With Quote
Old 02-15-2013, 08:52 AM   #2
Richard Finney
Senior Member
 
Location: bethesda

Join Date: Feb 2009
Posts: 700
Default

Try blseq
You can do it online.
http://blast.ncbi.nlm.nih.gov/Blast....SPEC=blast2seq
Richard Finney is offline   Reply With Quote
Old 02-15-2013, 08:58 AM   #3
ramirob
Member
 
Location: Vermont

Join Date: Apr 2012
Posts: 14
Default

Thnaks Richard!

We are talking about huge genome alignments. The samples have already been aligned to reference (this was part of an NGS experiment). My understanding is that bl2seq is more for aligning smaller sequences that have not been aligned before.
ramirob is offline   Reply With Quote
Old 02-15-2013, 09:20 AM   #4
Richard Finney
Senior Member
 
Location: bethesda

Join Date: Feb 2009
Posts: 700
Default

How huge?

By the way, there is a command line version for those big batch jobs. It might also handle large jobs that web version can't (maybe?).
Richard Finney is offline   Reply With Quote
Old 02-16-2013, 01:05 AM   #5
Liam_Gallagher
Member
 
Location: Europe

Join Date: Oct 2011
Posts: 18
Default

Hmm....I would try to use the WT (assembled de-novo) as reference for aligning the mutant and then use a variant caller, to identify variants between two genomes (using the WT as reference).

Also you can take a look at Jspecies. It is a java-based tool that measure the probability if two genomes belonging to the same species or not....and I think you can use your multifasta sequences as input and test your samples for differences. Basically it computes some blast between the samples and for each read you should see the percent of similarity. Then you can investigate the fragments with low similarity....
Liam_Gallagher 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 07:32 PM.


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