code: def plotall(xargs):
      
      p = OptionParser(plotall.__doc__)
      add_allmaps_plot_options(p)
      opts, args, iopts = p.set_image_options(xargs, figsize="10x6")
  
      if len(args) != 1:
          sys.exit(not p.print_help())
  
      inputbed, = args
      pf = inputbed.rsplit(".", 1)[0]
      agpfile = pf + ".chr.agp"
      agp = AGP(agpfile)
      objects = [ob for ob, lines in agp.iter_object()]
      for seqid in natsorted(objects):
          plot(xargs + [seqid])
docs: %prog plotall input.bed
  
      Plot the matchings between the reconstructed pseudomolecules and the maps.
      This command will plot each reconstructed object (non-singleton).