get_homologues how remove a genome from future runs

Question

[cloutier.dani12]

I've been using get_homologues for E. coli pangenome analysis of approximately 20 genomes comprising different ecotypes. As input for the get homologues.pl script I have been using a directory containing my genomes in genbank format. As the manual states: "The advantage of this option is that new files can be added to the input folder in the future and previous calculations will be conserved". As my genomes of interest became available, I was able to add them to the directory without having to re-run the previous comparisons. The problem that I'm having is when I try to remove a genome from the directory. When the analysis is re-run, an error message indicates that the removed genome cannot be found and the entire analysis has to be recomputed. Is there are way to by pass this? Thanks!

Answer 1

[bcontreras]

Hi Dani, option -I allows you to pass a list of genomes/proteomes to be analyzed which must be a subset of those in the input folder. You should this option for your purpose, instead of removing the original input file, since that breaks the internal renaming of sequences as integers. Hope this helps. Please use the issue manager at github (https://github.com/eead-csic-compbio/get_homologues/issues) for these questions, cheers, Bruno