2
$\begingroup$

I'm new to bioinformatics, and I'm looking for a tool that takes two FASTA files, one containing a reference, let's call it A.fasta, the other another sequence, let's call it B.fasta, and outputs just the mutations that B contains with respect to A.

I'm open as to the output format. Is there a standard for this?

Background: I'm trying to analyse Sars-CoV-2 sequences and identify mutations with respect to the reference. I've done a google search for a tool that does what I want to but couldn't find something straightaway. Is there no standard tool for this?

Improve this question
$\endgroup$
1
  • $\begingroup$ I would say parsimony, but thats just me $\endgroup$
    – M__ ♦
    Feb 25, 2022 at 20:57

1 Answer 1

Reset to default
2
$\begingroup$

If it's just a comparison of two sequences that you're looking for, then diffseq from the emboss toolkit should be okay:

https://www.bioinformatics.nl/cgi-bin/emboss/diffseq

I think this tool requires the sequences to be oriented in the same direction. If you discover a lot more variants than you expect, try reverse-complementing one of the sequences.

Improve this answer
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.