as(tag_cluster, "CTSS") should use the dominant CTSS information.

[charles-plessy]

Ideally, the result should be on position 26453667, not 26453605.

> p[4]
TagClusters object with 1 range and 6 metadata columns:
    seqnames            ranges strand |    dominant_ctss tpm.dominant_ctss            score   nr_ctss min_density max_density
       <Rle>         <IRanges>  <Rle> |        <GRanges>             <Rle>            <Rle> <integer>   <numeric>   <numeric>
  4    chr17 26453605-26453749      + | chr17:26453667:+  288.917306039892 1090.69591575985        18   0.0110095    0.544149
  -------
  seqinfo: 26 sequences (1 circular) from danRer7 genome
> p[4] |> as('CTSS')
CTSS object with 1 position and 6 metadata columns:
      seqnames       pos strand |    dominant_ctss tpm.dominant_ctss            score   nr_ctss min_density max_density
         <Rle> <integer>  <Rle> |        <GRanges>             <Rle>            <Rle> <integer>   <numeric>   <numeric>
  [1]    chr17  26453605      + | chr17:26453667:+  288.917306039892 1090.69591575985        18   0.0110095    0.544149
  -------
  seqinfo: 26 sequences (1 circular) from danRer7 genome
  BSgenome name: unspecified