docs: def get_merged_cds(genome):
      """
      Returns all CDS merged, used:
      - for TargQC general reports CDS coverage statistics for WGS
      - for Seq2C CNV calling when no capture BED available
      """
      bed = get_all_features(genome)
      debug('Filtering BEDTool for high confidence CDS and stop codons')
      return bed\
          .filter(lambda r: r.fields[BedCols.FEATURE] in ['CDS', 'stop_codon'])\
          .filter(high_confidence_filter)\
          .merge()