I am writing a C++ program to analyze BAM files. I am new to C and C++. I was wondering if anyone had some examples using the Samtools C headers in C++. So far I have been using the awesome Bamtools API, but I want to try the faidx.h header to access fasta sequences.
I have read the documentation, but I learn best by reading through code.