iqbal-lab-org / gramtools

Genome inference from a population reference graph
MIT License
92 stars 15 forks source link

Add --min-allele-freq argument to build command #53

Open sm0179 opened 7 years ago

sm0179 commented 7 years ago

add --min_freq arg to build step to be passed to the vcf_to_linear_prg script, for filtering variants in the vcf by allele frequency when constructing the linear PRG

ffranr commented 7 years ago

@sm0179 What is the use case for this feature?

sm0179 commented 7 years ago

the vcf_to_linear_prg script has an option to exclude from the PRG build process variants that occur in the population under a specified frequency (min_freq). this is because sometimes we want only a PRG with common variants. this is all already functional, we just need the build cmd to pass this argument when it calls the vcf_to_linear_prg script

ffranr commented 7 years ago

That make sense, thank you for clarifying. May I suggest: --min-allele-freq

sm0179 commented 7 years ago

definitely