code: def plot_by_gene_and_domain(name, seqs, tax, id2name):
      
      for gene in set([seq[0] for seq in list(seqs.values())]):
          for domain in set([seq[1] for seq in list(seqs.values())]):
              plot_insertions(name, seqs, gene, domain, tax, id2name)
docs: plot insertions for each gene and domain