docs: def annotate_diamond(fasta_path: 'path to fasta input',
                       diamond_path: 'path to Diamond taxonomic classification output'):
      '''
      Annotate fasta headers with taxonomy information from Diamond
      '''
      records = tictax.parse_seqs(fasta_path)
      annotated_records = tictax.annotate_diamond(records, diamond_path)
      SeqIO.write(annotated_records, sys.stdout, 'fasta')