playing with GNU #guile and #htslib.

(define with_vcf (lambda (f action) (begin
(define header (bcf_hdr_read f))
(define b (bcf_init1))
(while (>= (bcf_read fp header b) 0) (action header b))
(bcf_hdr_destroy header)
(bcf_destroy b)
(hts_close f)
)))

#bioinformatics #guile #scheme #vcf #lisp

htslib 1.23 released https://github.com/samtools/htslib/releases/tag/1.23

#bioinformatics #htslib

#htslib #Bioinformatics #GenomicsIO
@yokofakun any idea what is the fastest method to get the nth bcf record using htslib or bcftools.h without explicit loops? (trying out something similar to your project https://github.com/lindenb/rbcf but with ALTREP)
I guess it is possible to guess which block to (lazy) parse if one know the blocks offsets and number of records per block

via Stephen Turner on X : "Everyone's giving NotebookLM their papers, theses, marketing materials, etc. to generate these "Podcasts".

I gave it the full text of kseq.h from seqtk. "

https://notebooklm.google.com/notebook/39e3670a-dfa5-4808-aa11-06ab8c0cfb9d/audio?original_referer=https:%2F%2Ft.co%23&pli=1

#IA #podcast #htslib #bioinformatics

new program in #htslib "annot-regs" https://raw.githubusercontent.com/pd3/utils/master/annot-regs/doc.annot-regs.pdf

PR accepted ! option "--drop-genotypes" in #bcftools concat https://github.com/samtools/bcftools/pull/1911 🥳🥳

#bioinformatics #genomics #htslib