Variable-order reference-free variant discovery with the Burrows-Wheeler Transform