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.
•
written
by
jernej.kralj96 •
1 answer
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.
Log in to answer this question.
-
Help with Finding SNPs in H. pylori Assembled Genomes
written by vrao •
-
Clustal Omega Installation on M1 Mac (Ventura 13.1)
written by mboursnell •
-
Running SNIPPY on assembled genomes (calling variations from assembled genomes)
written by Kausik •
-
What reference to be used in snippy
written by Optimist
-
Is there a hybrid aligner that uses both long and short reads to align to a reference genome?
written by NPL •
-
Illumina polishing after PacBio assembly and polishing?
written by ikangkim •
-
how many genomes from BLAST hit should be considered for phylogeny tree analysis?
written by DanielC
-
Unicycler generating short contigs during hybrid assembly
written by hermoinejean •
-
Clustal Omega and bootstrapping
written by Antonio R. Franco
-
Assembling PacBio reads using a reference genome
written by celeste •
