This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to determine bacteria relatedness from two variant calls

When the NGS reads of an unknown bacteria are mapped to two reference genomes for variant calls, how do we determine the relatedness the unknown bacteria to each reference genome? From SNPs, assuming mappings fully cover two reference genomes, is there a metric that determines which reference is closer to the unknown bacteria?

next-gen genome snp

1 answer

You could try SNPRelate (https://bioconductor.org/packages/release/bioc/html/SNPRelate.html). You can create a PCA plot of the samples based on their VCF files.

Thank you very much for your suggestions!

Log in to answer this question.