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

Issue with bcftools annotate #2297

Open
HVHunt opened this issue Sep 30, 2024 · 1 comment
Open

Issue with bcftools annotate #2297

HVHunt opened this issue Sep 30, 2024 · 1 comment

Comments

@HVHunt
Copy link

HVHunt commented Sep 30, 2024

I am trying to use bcftools annotate to rename all positions in a vcf (basically, they all need to be smaller numbers to enable the vcf to be .tbi-indexed).

I have a tab-delimited file key.txt.gz formatted like this:
CM022211.1 2329174 1
CM022211.1 2329175 2
CM022211.1 2329176 3
CM022211.1 2329177 4

And I have tried the command bcftools annotate -a key.txt.gz -c CHROM,POS,~POS in.vcf.gz after reading this thread #1545.
However, the positions remain unchanged in the output file.
Please could you advise what I’m doing wrong as after a lot of searching and trying various syntaxes, I cannot figure this out?

@pd3
Copy link
Member

pd3 commented Oct 12, 2024

Which version of bcftools are you running? The notation ~POS was deprecated in favor of -POS in the latest version 1.21.

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

No branches or pull requests

2 participants