The recently added -H option for samtools bedcov generates a header describing the columns, but it only includes the chromosome, start, and end fields. If additional fields are present in the BED file, they are printed in the output but not reflected in the header. As a result, the column names in the header do not correspond to all the output columns.

Tested with the latest version, Samtools 1.21.

Thanks.