This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to get list of different bases and their positions from an alignment?

I am aligning two viral genomes which I am expecting to be >90% similar. I would like a list of the changes (from one base to another) and their locations in the genome. Which program(s) can I use to achieve this? I have already tried lastz's --census option, but it only tells you if a base is aligned, not whether it is identical. Thanks.

alignment
• 2,016 views • link
updated by WouterDeCoster WouterDeCoster • written by SaltedPork SaltedPork

Hello,

is your data from a sequencing project or do you just have the reference sequences you like to compare?

fin swimmer

• 1 views • link
written by finswimmer finswimmer

Hi, just the sequences i'm comparing.

• 1 views • link
written by SaltedPork SaltedPork

Hi,

how have you done the alignment? How does the output look like?

fin swimmer

• 1 views • link
written by finswimmer finswimmer

they are about 92% similar, ~700 base difference

• 1 views • link
written by SaltedPork SaltedPork

That does not answer @finswimmer's question about format of the output.

• 1 views • link
written by GenoMax GenoMax

Hi SaltedPork,

I have reopened this question, there is not really a need to delete it.

Cheers,
Wouter

• 1 views • link
written by WouterDeCoster WouterDeCoster

1 answer

If I understand correctly you want to align to reference genomes and find the differences. If so, you could try minimap2 to align both using the -cx asm5 option. Scripts are available in that repo (paftools) to get the differences.

• 1 views • link
written by WouterDeCoster WouterDeCoster

Log in to answer this question.

More posts like this