I would like to do a pairwise/multiple sequence alignment for a gene from two/three species and then record in a separate file the exact positions of perfect identity between the sequences. Is there a tool that does this already? Or is there any option in blast/clustal that could help in recovering this information ?
Will be very thankful for your time and responses.