Code snippet for including Erik Garrison's tabixpp located at https://github.com/ekg/tabixpp

readTabix.cpp

#include "tabix.hpp"

void anyfunc(){

  //tabixed vcf
  string inVCF = "input.vcf.gz";
  
  // region 
  string region = "20:1-4000000";
  
  // create tabix object
  Tabix tbx(inVCF);
  
  // set region to extract from vcf
  if (tbx.setRegion(region) != true)
    throw std::runtime_error("Malformed region string: " + region);
  
  // get the header
  string header;
  tbx.getHeader(header);
  
  // read each line in region
  while (tbx.getNextLine(line)) {
    // parse vcf record as string
  }
  
}