Toward improved calling of allele-specific events in high-throughput sequencing data with a mixture of compound Dirichlet negative multinomial distributions