Removing a symbolic allele from VCF header
Nobody has claimed this yet.
Assessment
- Difficulty
- 5/5
- Estimated time
- Over a week
- Newbie friendliness
- 30/100
- Issue type
- Feature
- Clarity
- Needs clarification
- Activity status
- Stale
- Tech stack
- c
- Domain
- bioinformatics
Research direction
Start with bcf_hdr_remove and the referenced vcf.c sections around lines 755 and 761, then review how BCF_HL_GEN and BCF_HL_STR represent ##ALT entries. Determine whether a direct symbolic-allele removal operation is needed; done means the intended API behavior is agreed and covered for VCFv4.1 headers.
Written by the indexing model from the issue text.
Description
I'm trying to remove the VCF header line for a symbolic alternate allele. Since ##ALT header lines belong to a specific BCF_HL_* type, it seems like my only options are removing all ##ALT= lines, or remove using the ID. Does this look right? Am I missing a more direct way?
bcf_hdr_append(hdr, "##ALT=<ID=DEL,Description=\"Deletion\">");
Remove all ##ALT= lines:
bcf_hdr_remove(hdr, BCF_HL_GEN, "ALT");
Remove by matching to the ID.
bcf_hdr_remove(hdr, BCF_HL_STR, "DEL");
Currently working in VCFv4.1
- Dominant language
- C
- Stars
- 950
- Forks
- 475
- Avg merge
- 3d 13h
- Merged PRs (30d)
- 11
Contributor guide
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
More from samtools/htslib
-
Difficulty 3/5 1-2 days Newbie friendliness 74/100
-
Difficulty 3/5 1-2 days Newbie friendliness 76/100
-
tabix performance Open
Similar issues
-
[adam] AdamNet network read doesn't cap to MAX_ADAM_PACKET_LEN, overflows client receive buffers Open
Difficulty 2/5 1-3 hours Newbie friendliness 88/100
FujiNetWIFI/fujinet-firmware#1649 · 2 comments ·
-
Difficulty 2/5 1-3 hours Newbie friendliness 68/100
HarbourMasters/Shipwright#7229 ·
-
Difficulty 2/5 1-3 hours Newbie friendliness 72/100
riscv-software-src/riscv-isa-sim#2435 · 1 comment ·
-
bug Self Built Image SNAPSHOT Supported Device target/ramips
Difficulty 2/5 1-3 hours Newbie friendliness 88/100
-
Difficulty 2/5 1-3 hours Newbie friendliness 76/100