How to resolve the issue of wgbstools mix_pat not being able to parse the hg38 genome?

[MyXu0520]

Hello, loyfer,

Thanks for developing this useful tool.

Recently, while using the wgbstools , I encountered a minor issue and would like to ask for advice on how to resolve it. When I tested the input /home/xmy/anaconda3/envs/cfmetrace/bin/wgbstools mix_pat /biocc/MyXu/cfDNA_deconv_benchmark-main/my_code_output/3/mix/GSM5652176_u/cfDNA_deconv_benchmark-main/my_code_output/3/mix/GSM5652176_Adipocytes-Z000000T7.pat.gz --rates 0.061744 0.0 -p /biocc/MyXu/cfDNA_decon-genome hg38 -c 22, I encountered the error FileNotFoundError: [Errno 2] No such file or directory: '/home/xmy/anaconda3/envs/cfmetrace/references/default'. However, I had already built the reference files using wgbstools init_genome hg38, and the files are located at /home/xmy/anaconda3/envs/cfmetrace/references/hg38, containing the hg38.fa.gz file. I also tried downloading the hg38.fa.gz file from UCSC myself, but the same error was reported. How should I resolve this issue?

Many thanks,

Best,
Manyi Xu

[yonniejon]

does wgbstools set_default_ref hg38 work?