def plot_by_gene_and_domain(name, seqs, tax, id2name):
for gene in set([seq[<NUM_LIT:0>] for seq in list(seqs.values())]):<EOL><INDENT>for domain in set([seq[<NUM_LIT:1>] for seq in list(seqs.values())]):<EOL><INDENT>plot_insertions(name, seqs, gene, domain, tax, id2name)<EOL><DEDENT><DEDENT>
plot insertions for each gene and domain