Divide by 0 error when passing multiple files to plot-bcfstats
0
0
Entering edit mode
7.3 years ago

I'm trying to run bcftools stats on multiple files and then have plot-vcfstats give me a single plot that shows the data for each file. If I run plot-vcfstats on one file, it's fine. If I try to run it on multiple files, it breaks.

Input:

plot-vcfstats -p Chr_all/ Chr12.vcf.gz_stats.txt ChrX.vcf.gz_stats.txt

Error message:

Parsing bcftools stats output: Chr12.vcf.gz_stats.txt
Parsing bcftools stats output: ChrX.vcf.gz_stats.txt
Illegal division by zero at /usr/local/Cluster-Apps/bcftools/1.2/bin/plot-vcfstats line 480.

This happens regardless of what combination of input files I use. All input files I have tested will run fine in plot-vcfstats individually, but will throw this divide by zero exception if I pass more than one as an argument to plot-vcfstats

I have tried merging them in the hope of avoiding this, but got the exact same error as above:

Illegal division by zero at /usr/local/Cluster-Apps/bcftools/1.2/bin/plot-vcfstats line 480.

Command for merge:

plot-vcfstats -m Chr12.vcf.gz_stats.txt ChrX.vcf.gz_stats.txt

Can anyone help me sort this out, or suggest a workaround if this is a bug?

Thanks,

-James

vcf bcftools matplotlib • 1.4k views
ADD COMMENT

Login before adding your answer.

Traffic: 1538 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6