code: def comparebed(args):
      
      p = OptionParser(comparebed.__doc__)
      opts, args = p.parse_args(args)
  
      if len(args) != 2:
          sys.exit(not p.print_help())
  
      abed, bbed = args
      abed = Bed(abed)
      bbed = Bed(bbed)
      query_links(abed, bbed)
      query_links(bbed, abed)
docs: %prog comparebed AP.chr.bed infer.bed
  
      Compare the scaffold links indicated in two bed files.