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
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
Comment