Commands: discover

This command uses the reads passed in quasimap and the personalised reference genome found in infer to discover new variants, using Cortex, an assembly-based variant caller.

The variants discovered get rebased, ie expressed in terms of the VCF from infer. This means you can run gramtools build using the original VCF and the one from this command, to augment the original graph with the newly found variation.