Skip to content

Commit

Permalink
Merge pull request #59 from gymrek-lab/test-simgenotype
Browse files Browse the repository at this point in the history
minor cleanups to simgenotype output messages
  • Loading branch information
mlamkin7 authored Jun 30, 2022
2 parents 062afdc + f021d3a commit 9e63e06
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 1 deletion.
3 changes: 2 additions & 1 deletion haptools/__main__.py
Original file line number Diff line number Diff line change
Expand Up @@ -70,9 +70,10 @@ def simgenotype(invcf, sample_info, model, mapdir, out, popsize, seed, chroms):
Simulate admixed genomes under a pre-defined model.
Example:
haptools simgenotype \
--model ./tests/data/outvcf_gen.dat \
--mapdir ./tests/map/ \
--mapdir ./tests/data/map/ \
--chroms 1,2 \
--invcf ./tests/data/outvcf_test.vcf \
--sample_info ./tests/data/outvcf_info.tab \
Expand Down
2 changes: 2 additions & 0 deletions haptools/sim_genotype.py
Original file line number Diff line number Diff line change
Expand Up @@ -34,6 +34,8 @@ def output_vcf(breakpoints, model_file, vcf_file, sampleinfo_file, out):
output prefix
"""

print(f"Outputting VCF file {vcf_file}")

# details to know
# vcf file: how to handle samples and which sample is which haplotype block randomly choose out of current population types
# need to go line by line of the vcf when creating the new vcf
Expand Down

0 comments on commit 9e63e06

Please sign in to comment.