 
         uint2 vcf_db_range = ctx.snp_filter.active_vcf_ranges[read_index];
 
-
-
         // traverse the VCF range and mark corresponding read BPs as inactive
         for(uint32 feature = vcf_db_range.x; feature <= vcf_db_range.y; feature++)
         {

        vcf_range.x = vcf_start - db.reference_window_start.begin();
        vcf_range.y = vcf_range.x;


        // do a linear search to find the end of the VCF range
        // (there are generally very few VCF entries for an average read length --- and often none --- so this is expected to be faster than a binary search)
        while(vcf_range.y < db.reference_window_start.size() - 1 &&

        : bqsr_lambda(ctx, batch)
    { }

public:
    CUDA_HOST_DEVICE void operator() (const uint32 read_index)
    {
