Hi, I'm trying to run bcftools stats on multiple chromosomes and then have plot-vcfstats give me a single plot that shows the data for each of of the chromosomes. I've got one vcf file per chromosome. If I run plot-vcfstats on one file, it's fine. If I try to fun on multiple files, it breaks.
Input: plot-vcfstats -p Chr_all/ Chr12-1-133851895.vcf.gz_stats.txt ChrX-1-155270560.vcf.gz_stats.txt Error message: Parsing bcftools stats output: Chr12-1-133851895.vcf.gz_stats.txt Parsing bcftools stats output: ChrX-1-155270560.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. Command for merge: plot-vcfstats -m Chr12-1-133851895.vcf.gz_stats.txt ChrX-1-155270560.vcf.gz_stats.txt Can anyone help me sort this out, or suggest a workaround if this is a bug? -James
------------------------------------------------------------------------------
_______________________________________________ Samtools-help mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/samtools-help
