Add functionality to calculate expected coverage of kmers. Fix vcf in ctx_out name