SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Bioinformatics tool/method to accurately identify plasmid/phages in bacterial genome EulnayM Bioinformatics 0 11-06-2014 09:37 PM
Tool to take mapping to gene list jarthur Bioinformatics 1 01-16-2013 03:31 PM
Tool for short nucleotide sequences Dolphin22 General 1 09-27-2012 10:01 PM
Getting a list of all index sequences Mouth_Breather Illumina/Solexa 3 07-12-2012 02:26 AM
PubMed: PositionMatcher: A Fast Custom-Annotation Tool for Short DNA Sequences. Newsbot! Literature Watch 1 02-28-2011 08:00 AM

Reply
 
Thread Tools
Old 02-08-2015, 03:31 AM   #1
Antony03
Member
 
Location: Canada, Quebec

Join Date: Apr 2012
Posts: 53
Default Is a bioinformatics tool to list mutations between two short bacterial sequences?

Hello,

I have two bacterial sequences of 10 kbp in fasta format and I am looking for a bioinformatics tool which is able to list the mutations (Indels and point mutations) between both sequences.

I performed a global pairwise alignment of my sequences using needle (Needleman-Wunsch algorithm) from the Emboss package, which is known to produces an optimal alignment.

In a second time, I used MUMmer as indicated in the manual for snps/indels dectection: http://mummer.sourceforge.net/manual/#snpdetection

However, there some differences (mainly for indels) between the alignment produced by needle (that I use as a quality control) and the result generated by MUMmer. I suspect that the differences come from the initial alignment produced by MUMmer.

Is there another tool (like MUMmer) for which I can give my own alignment?

Thank
Antony03 is offline   Reply With Quote
Old 02-08-2015, 09:31 AM   #2
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

Quote:
I performed a global pairwise alignment of my sequences using needle (Needleman-Wunsch algorithm) from the Emboss package, which is known to produces an optimal alignment.
Needleman-Wunsch does not yield optimal alignments. For optimality, under a model that assumes fewer mutation events are more likely than more mutation events, you need affine transforms to encourage insertions or deletions to be made adjacent rather than scattered around at random, which Needleman-Wunsch or Smith-Waterman will do because they are greedy algorithms.

That's just an aside, though. I don't know what tools would help you in this case.
Brian Bushnell 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 01:48 AM.


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