Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Script to generate template VCF file for assembly #166

Closed
timothymillar opened this issue Jul 11, 2023 · 2 comments
Closed

Script to generate template VCF file for assembly #166

timothymillar opened this issue Jul 11, 2023 · 2 comments
Assignees
Labels
enhancement New feature or request
Milestone

Comments

@timothymillar
Copy link
Collaborator

Simple script to generate a template VCF for MCHap assemble.

There are several motivations for this:

  • Simplicity, we only need position and alleles of putative segregating sites
  • Speed, less compute and using same input files
  • Ability to mask reference allele to lower parameter space Allow masking of SNV reference alleles #165
@timothymillar timothymillar added the enhancement New feature or request label Jul 11, 2023
@timothymillar timothymillar self-assigned this Jul 11, 2023
@timothymillar timothymillar added this to the v0.9.0 milestone Jul 11, 2023
timothymillar added a commit to timothymillar/MCHap that referenced this issue Jul 12, 2023
timothymillar added a commit to timothymillar/MCHap that referenced this issue Jul 12, 2023
timothymillar added a commit to timothymillar/MCHap that referenced this issue Jul 12, 2023
timothymillar added a commit to timothymillar/MCHap that referenced this issue Jul 12, 2023
timothymillar added a commit to timothymillar/MCHap that referenced this issue Jul 12, 2023
timothymillar added a commit to timothymillar/MCHap that referenced this issue Jul 13, 2023
timothymillar added a commit to timothymillar/MCHap that referenced this issue Jul 13, 2023
@timothymillar
Copy link
Collaborator Author

timothymillar commented Oct 20, 2023

Idea for arguments:

  • --maf: min population allele frequency (default = 0.0)
  • --mad: min population allele depth (default = 0)
  • --ind-maf: min individual allele frequency (default = 0.1)
  • --ind-mad: min individual allele depth (default = 3)
  • --ind-min: min number of individuals passing --ind-maf and --ind-mad (default = 1)

All filters are combined (intersection). This gives good flexibility and makes few assumptions about population size/structure for defaults.

@timothymillar
Copy link
Collaborator Author

Fixed in #169

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

1 participant