Re: [Samtools-help] segfault: bcftools stats

2016-05-27 Thread John Marshall
On 27 May 2016, at 06:14, Peter Johansson wrote: > Seems the problematic entry had a negative QUAL (-2147480064), which looks > very suspicious. Need to go upstream and look for the source of that. QUAL is Phred-scaled, so by definition valid values are non-negative. So the upstream source mig

Re: [Samtools-help] segfault: bcftools stats

2016-05-26 Thread Peter Johansson
Hi again, Seems the problematic entry had a negative QUAL (-2147480064), which looks very suspicious. Need to go upstream and look for the source of that. I have no idea why negative QUAL causes a segfault and whether it would make sense with a more gracious exit. I have attached the problema

[Samtools-help] segfault: bcftools stats

2016-05-26 Thread Peter Johansson
Hi, I get a segfault when I run 'bcftools stats'. Not sure if it's a bug or incorrect input data. This is with bcftools 1.3.1 with htslib 1.3.1 The command is: bcftools stats --fasta-ref /usr/local/share/genomes/human_g1k_v37.fasta KI006.anno.vcf.gz I get the backtrace below: Program recei