Skip to content

Commit

Permalink
adjust output filename
Browse files Browse the repository at this point in the history
  • Loading branch information
jethror1 committed Sep 23, 2024
1 parent 20aa472 commit 5bcaff7
Show file tree
Hide file tree
Showing 2 changed files with 2 additions and 2 deletions.
2 changes: 1 addition & 1 deletion README.md
Original file line number Diff line number Diff line change
Expand Up @@ -11,7 +11,7 @@ This app should be executed stand-alone or as part of a DNAnexus workflow for a
The app requires a VCF file (.vcf) containing variants to be evaluated and a bed file defining the regions within which variants in the vcf should be evaluated.

## What does this app output?
The app outputs one tab delimited file `[vcf_prefix].vcf.qc`, where `vcf_prefix` is the vcf filename without extension:
The app outputs one tab delimited file `[vcf_prefix].vcf.QC`, where `vcf_prefix` is the vcf filename without extension:

* sample id
* mean het ratio (mean AAF of het variants)
Expand Down
2 changes: 1 addition & 1 deletion src/vcf_qc.py
Original file line number Diff line number Diff line change
Expand Up @@ -300,7 +300,7 @@ def main(vcf_file, bed_file):
het_hom_counts = get_het_hom_counts(vcf=tmp_vcf)
ratios = calculate_ratios(counts=het_hom_counts)

outfile = f"{re.sub(r'.vcf(.gz)?$', '', vcf_file)}.vcf.qc"
outfile = f"{re.sub(r'.vcf(.gz)?$', '', vcf_file)}.vcf.QC"

if os.path.exists("/home/dnanexus"):
write_output_file(outfile=outfile, ratios=ratios)
Expand Down

0 comments on commit 5bcaff7

Please sign in to comment.