code: def explore(args):
      
      logger.info("reading sequeces")
      data = load_data(args.json)
      logger.info("get sequences from json")
      #get_sequences_from_cluster()
      c1, c2 = args.names.split(",")
      seqs, names = get_sequences_from_cluster(c1, c2, data[0])
      loci = get_precursors_from_cluster(c1, c2, data[0])
      logger.info("map all sequences to all loci")
      print("%s" % (loci))
      map_to_precursors(seqs, names, loci, os.path.join(args.out, "map.tsv"), args)
      #map_sequences_w_bowtie(sequences, precursors)
      logger.info("plot sequences on loci")
      #get_matrix_position()
      #plot_sequences()
      logger.info("Done")
docs: Create mapping of sequences of two clusters