docs: def dalignbed2annotationsbed(cfg):
      """
      Get annotations from the aligned BED file
      step#3
  
      :param cfg: configuration dict
      """
      datatmpd=cfg['datatmpd']
      alignmentbedp=cfg['alignmentbedp']    
      alignmentbedsortedp=alignmentbedp+'.sorted.bed'
      logging.info(basename(alignmentbedsortedp))
      if not exists(alignmentbedsortedp) or cfg['force']:
          cmd=f"{cfg['bedtools']} sort -i {alignmentbedp} > {alignmentbedsortedp}"
          runbashcmd(cmd)
      
      genomegffsortedp=cfg['genomegffp']+'.sorted.gff3.gz'
      logging.info(basename(genomegffsortedp))
      if not exists(genomegffsortedp):    
          cmd=f"{cfg['bedtools']} sort -i {cfg['genomegffp']} > {genomegffsortedp}"
          runbashcmd(cmd)
  
      annotationsbedp=f'{datatmpd}/03_annotations.bed'
      cfg['annotationsbedp']=annotationsbedp
      logging.info(basename(annotationsbedp))
      if not exists(annotationsbedp) or cfg['force']:    
          cmd=f"{cfg['bedtools']} intersect -wa -wb -loj -a {alignmentbedsortedp} -b {genomegffsortedp} > {annotationsbedp}"
          runbashcmd(cmd)
      return cfg