This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Get the divergence between two species with a ref genome and reads

Hi !

I was wondering if it was possible to get the mean divergence (in percentage of nucleotide) between two species using a genome reference and reads. I explain my self :

I have a genome of reference for the species 1 and i have illumina reads for the species 2. I have mapped those reads to the reference genome and produced a VCF file. I would like now to get the mean divergence between those two species (for example 1% of bases are not the same etc...) and is it biologically relevant given that the reference genome is more likely to not be completely covered ?

Is it possible ?

Thanks for your help,

Maxime Policarpo

assembly vcf mapping

2 answers

Divergence is not a uniform value across the genome: coding regions are restricted by stronger selection force against fast and big changes over time, while non-coding regions are less constrained. Given you could map reads from one species to reference genome from another species, your comparison seems to be between two relatively related species. I'd suggest you to start with coding regions, following the standard practices in what is called phylogenomics, comparing gene sequences across species at a whole genome levels, and calculate divergences at the nucleotide level. Non-coding regions are a whole different ball game, you would need much more effort to find differences in non-coding regions, especially for structural variations.

Thanks for your answer ! I will try to do that :D

Thanks for your answer ! I will try to do that :D

Log in to answer this question.