Is there a way to have bcftools stats output Non-Reference Discordance by non-reference allele frequency?

[njbernstein]

It would be great if bcftools stats could generate this output as well. I dont see it in the docs nor when I use bcftools stats. Tool works great though thanks!

Here is my invocation:

bcftools stats --af-tag "RAF" \
                --af-bins 0.00000,0.00100,0.00200,0.00500,0.01000,0.05000,0.10000,0.20000,0.50000,1 \
                -s ${sample_name} ${sample}.imputed.chr${chrom}.bcf ${sample_name}.chr${chrom}.1000genomes.bcf \
                > bcftools_stats/${sample_name}.chr${chrom}.bcfstats.txt