docs: def matrix(fasta_path: 'path to tictax annotated fasta input',
             scafstats_path: 'path to BBMap scaftstats file'):
      '''
      Generate taxonomic count matrix from tictax classified contigs
      '''
      records = SeqIO.parse(fasta_path, 'fasta')
      df = tictax.matrix(records, scafstats_path)
      df.to_csv(sys.stdout)