hap.py
hap.py copied to clipboard
Get something error in running the example.
Hi, I install the hap.py through anaconda. When I run the test run: $ hap.py example/happy/PG_NA12878_chr21.vcf.gz example/happy/NA12878_chr21.vcf.gz -r example/chr21.fa -f example/happy/PG_Conf_chr21.bed.gz -o test
I get something error. May someone help me? Thank you a million.
2021-08-08 20:27:09,450 WARNING No reference file found at default locations. You can set the environment variable 'HGREF' or 'HG19' to point to a suitable Fasta file.
2021-08-08 20:27:09,450 WARNING No reference file found at default locations. You can set the environment variable 'HGREF' or 'HG19' to point to a suitable Fasta file.
[W] overlapping records at chr21:10993857 for sample 0
[W] Symbolic / SV ALT alleles at chr21:15847469
[W] Variants that overlap on the reference allele: 144
[W] Variants that have symbolic ALT alleles: 14
[I] Total VCF records: 65402
[I] Non-reference VCF records: 65402
[W] overlapping records at chr21:24024261 for sample 0
[W] Variants that overlap on the reference allele: 5
[I] Total VCF records: 101524
[I] Non-reference VCF records: 101524
2021-08-08 20:27:40,250 ERROR [stderr] regex_error
2021-08-08 20:27:40,257 ERROR Command 'quantify /tmp/hap.py.result.mkd_6q.vcf.gz -o test.roc.tsv -r example/chr21.fa --threads 16 --output-vtc 0 --output-rocs 1 --type xcmp --qq IQQ --qq-header QUAL --roc-delta 0.500000 --clean-info 1 --fix-chr-regions 0 -v test.vcf.gz -R 'CONF_VARS:/tmp/tmp9lwgPm.bed' -R 'CONF:example/happy/PG_Conf_chr21.bed.gz' --roc-regions ''' returned non-zero exit status 1
2021-08-08 20:27:40,258 ERROR Traceback (most recent call last):
2021-08-08 20:27:40,258 ERROR File "/home/user/.conda/envs/happy/bin/hap.py", line 511, in