ga4gh / ga4gh/benchmarking-tools

CalledProcessError

Open
#33 0 comments 0 reactions 0 assignees View on GitHub
Dominant language
HTML
Stars
213
Forks
47
PR merge metrics
No merged PRs in 30d

Description

I am getting a CalledProcessError regarding the benchmark vcf, even though it has been indexed.

```
2021-05-05 18:27:33,391 WARNING No reference file found at default locations. You can set the environment variable 'HGREF' or 'HG19' to point to a suitable Fasta file.
Hap.py v0.3.14
[E::vcf_hdr_read] no sample line
**_**Failed to open or file not indexed: /u/home/n/neharajk/scratch/HG002_GRCh37_1_22_v4.2.1_benchmark_19.vcf**_**

2021-05-05 18:27:41,571 ERROR Command 'vcfcheck /u/home/n/neharajk/scratch/HG002_GRCh37_1_22_v4.2.1_benchmark_19.vcf --check-bcf-errors 1' returned non-zero exit status 1
2021-05-05 18:27:41,571 ERROR Traceback (most recent call last):
2021-05-05 18:27:41,571 ERROR File "./hap.py", line 540, in
2021-05-05 18:27:41,573 ERROR main()
2021-05-05 18:27:41,573 ERROR File "./hap.py", line 316, in main
2021-05-05 18:27:41,573 ERROR convert_gvcf_to_vcf=convert_gvcf_truth)
2021-05-05 18:27:41,573 ERROR File "/u/scratch/m/mdistler/hap.py-build/bin/pre.py", line 125, in preprocess
2021-05-05 18:27:41,576 ERROR mf = subprocess.check_output("vcfcheck %s --check-bcf-errors 1" % pipes.quote(vcf_input), shell=True)
2021-05-05 18:27:41,576 ERROR File "/u/local/apps/python/2.7.2/lib/python2.7/subprocess.py", line 544, in check_output
2021-05-05 18:27:41,599 ERROR raise CalledProcessError(retcode, cmd, output=output)
2021-05-05 18:27:41,600 ERROR CalledProcessError: Command 'vcfcheck /u/home/n/neharajk/scratch/HG002_GRCh37_1_22_v4.2.1_benchmark_19.vcf --check-bcf-errors 1' returned non-zero exit status 1
```

Any help?

Contributor guide

No contributing guide indexed for this repository

Research direction

Start with hap.py and the preprocess entry point in pre.py, especially the vcfcheck invocation shown in the traceback. Reproduce the command with the referenced benchmark VCF and determine whether the failure is caused by the file, its index, or missing VCF header content; done means identifying a reproducible cause and a documented resolution.

Written by the indexing model from the issue text.

Assessment

Tech stack
python
Domain
bioinformatics
Issue type
Bug
Difficulty
3/5
Estimated time
1-2 days
Activity status
Stale
Clarity
Needs clarification
Newbie friendliness
20/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.