Hi,
I am trying to run VariantDatasetCombiner like this:
path_to_input_list = ‘input_files.txt’ # a file with one GVCF path per line
gvcfs =
with hl.hadoop_open(path_to_input_list, ‘r’) as f:
for line in f:
print(line)
gvcfs.append(line.strip())
combiner = hl.vds.new_combiner(
output_path=‘dataset.vds’,
temp_path=‘my-temp-bucket’,
gvcf_paths=gvcfs,
use_genome_default_intervals=True,
reference_genome=‘GRCh38’,
)
combiner.run()
but an error ocures:
Current key: { locus: { contig: chrM, position: 302 }, alleles: [3; A, AC, ACC] }
Previous key: { locus: { contig: chrM, position: 302 }, alleles: [3; A, AC, C] }
This error can occur after a split_multi if the dataset
contains both multiallelic variants and duplicated loci.
Do you have some recomendation how what to change in my code and how I can run VariantDatasetCombiner?
hl.experimental.run_combiner was working for me with the following code:
output_file = ‘output.mt’ # output destination
temp_bucket = ‘my-temp-bucket’ # bucket for storing intermediate files
hl.experimental.run_combiner(inputs, out_file=output_file, tmp_path=temp_bucket, reference_genome=‘GRCh38’, use_genome_default_intervals=True, overwrite=True)
, but I would like to access metods that VariantDataset has.