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

Non-existant samples accepted by view #121

Open
jeromekelleher opened this issue Jan 14, 2025 · 1 comment
Open

Non-existant samples accepted by view #121

jeromekelleher opened this issue Jan 14, 2025 · 1 comment
Labels
bug Something isn't working

Comments

@jeromekelleher
Copy link
Contributor

 python3 -m vcztools view performance/data/sim_10k.vcz/ -s NO_SAMPLE | head
##fileformat=VCFv4.2
##FILTER=<ID=PASS,Description="All filters passed">
##source=tskit 0.6.0
##contig=<ID=1,length=50818468>
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##vcztools_viewCommand=view performance/data/sim_10k.vcz/ -s NO_SAMPLE; Date=2025-01-14 13:49:27.031019
#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  NO_SAMPLE
1       245     0       A       G       .       PASS    AC=0;AN=2       GT      0|0
1       427     1       A       C       .       PASS    AC=0;AN=2       GT      0|0
1       479     2       A       C       .       PASS    AC=0;AN=2       GT      0|0

correct output:

$ bcftools view performance/data/sim_10k.vcf.gz -s NO_SAMPLE
Error: subset called for sample that does not exist in header: "NO_SAMPLE". Use "--force-samples" to ignore this error.
@jeromekelleher jeromekelleher added the bug Something isn't working label Jan 14, 2025
@jeromekelleher
Copy link
Contributor Author

Looking through the implementation and trying to replicate the above silent failure bug, it's not at all clear to me what's going on. I think the implementation needs to be decoupled from the CLI encoding (i.e. list of strings) and tested in isolation.

I think the best thing for us to ship a working version would be to disable the sample subsetting functionality at the front-end and return to it once that's done.

jeromekelleher added a commit to jeromekelleher/vcztools that referenced this issue Jan 14, 2025
It has a number of problems and better not to try and fix them in a
hurry pre shipping. See sgkit-dev#121 and sgkit-dev#122
jeromekelleher added a commit that referenced this issue Jan 14, 2025
It has a number of problems and better not to try and fix them in a
hurry pre shipping. See #121 and #122
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

1 participant