This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Core genome analysis for viruses

Hello everyone! I've been trying to build core genomes of HIV and SIV viruses using Roary. But I keep getting 0 core genes. I've tried changing "minimum percentage identity for blastp" parameter. I used the default value first and here's the result :

Core genes  (99% <= strains <= 100%)    0
Soft core genes (95% <= strains < 99%)  0
Shell genes (15% <= strains < 95%)  0
Cloud genes (0% <= strains < 15%)   511
Total genes (0% <= strains <= 100%) 511

Then I changed it to 70 and here's what I got:

Core genes  (99% <= strains <= 100%)    0
Soft core genes (95% <= strains < 99%)  0
Shell genes (15% <= strains < 95%)  14
Cloud genes (0% <= strains < 15%)   78
Total genes (0% <= strains <= 100%) 92

I don't understand why shell and cloud genes numbers varied!! nor why the core genes number is 0! Is it because these viruses are highly mutable? Should I change the default value of other parameters? if that so , what s are the values of parameters i should use , also , Could anyone recommend me any other tool more specific for core genome viruses analysis to construct core genomes of these viruses? Many thanks

assembly sequence sequencing gene

0 answers

No answers yet.

Log in to answer this question.