problems with snippy on galaxy

Question

[jernej.kralj96]

Hello

I have a few questions regarding snippy. I need to analyze around 20 genomes and compare them to reference genome. I used unicycler to assemble my reference genome from illumina short reads and Pacbio long reads and then I used snippy to find changes in genomes. I am getting very confusing results, for instance in some of the genomes I get exact same mutations in the same position. For instance indel on a gene for immunity protein caused same change in almost every analzyed genome. Then I used prokka to annotate some of the genomes and clustal omega to align genes of interest to reference genome and found out that they are identical 100 % to reference genome. So what could be the isssues ? Can anyone explain.

Have a nice day.

Answer 1

[Istvan Albert]

Visualize your data in IGV, It is possible to get both a variation and identity if you have a population of cells some that match identically the references and subset does not.

Align your reads to the reference, align your assembly to a reference and see what it tells you.