'bcftools call' should output SNP/INDEL of same site to one vcf record

Open
#552 1 comment 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Assessment

Difficulty
4/5
Estimated time
3-5 days
Newbie friendliness
35/100
Issue type
Bug
Clarity
Mostly clear
Activity status
Stale
Tech stack
c

Research direction

Start by reproducing the report from a.zip with bcftools call, then run the shown bcftools norm -m+any and bcftools plugin fill-AN-AC pipeline to inspect the inconsistent FORMAT/GT values. Trace the bcftools call entry point and its handling of SNP/INDEL records at one site; done means the variants are emitted as one VCF record with consistent genotype information.

Written by the indexing model from the issue text.

Description

enhancement

'bcftools call' should output the SNP/INDEL of same site to one vcf record.

Now the SNP/INDEL of same site are output separately, and sometimes the FORMAT/GT info of two record are inconsistent.

a.zip

If we 'bcftools norm -m+any a.vcf | bcftools plugin fill-AN-AC', we can find this inconsistent in GT.

Dominant language
C
Stars
891
Forks
277
Avg merge
2d 23h
Merged PRs (30d)
2

Contributor guide

Open the contributing guide

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

More from samtools/bcftools

All issues in samtools/bcftools

Similar issues

More C issues

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.