samtools ampliconstats doesn't support multi reference sequence/contig genomes (known issue) but it would be really nice if could or if the lack of support was clearly documented and handled cleanly by the code (e.g. didn't seg fault).

[sanger only example data at /lustre/scratch117/sciops/team117/npg/ces/data/ampliconstats_multiseq/ ]