Comments
More comments
Merge procedure modified:
Pre-process individual files to
a. Move fields from INFO to FORMAT (per sample)
b. Drop un-necessary fields
Debug
1. Changed to -O3 flag in Makefile by default
2. Use processed_header in vcfmerge
Fixed memory leak in move_info_to_format
Added code for merge configuration file directives
bcf_hdr_merge now has force_samples argument
1. Debugged merge_config
2. Added sample config.sample for merge
vcfmerge works with gVCFs now. Needed a modification to synced_reader in
htslib.
Main idea: before merge_buffer fully kicks in, check whether the
intervals need to be split and aligned to a minimum end position. The
minimum is computed by the function find_end_position and alignment to
the end is performed by align_to_end. This process produces new
intervals. Through a tricky manipulation of reader->nbuffer and
reader->last_valid_record_idx, intervals are merged successfully.

FIXME:
1. The terrible hack that is is_split_record. Needs to go away as soon
as I figure out how to obtain base from reference given position
2. shake_buffer around line 1708, unsure what to do