#pragma once

// -----------------------------------------------------------------------
// Phylogenetics

typedef struct bio·Node
{
    string name;
    string comment;
    double dist;
    double support;
    int    nchild;
    int    ndescendent;
    struct bio·Node *parent;
    struct bio·Node *child[2];
    struct bio·Node *sibling;
} bio·Node;

error phylo·addchild(bio·Node *parent, bio·Node* child);
bool  phylo·isleaf(bio·Node *node);

typedef struct bio·Tree
{
    bio·Node *root;
} bio·Tree;

/* newick i/o */
bio·Tree bio·readnewick(io·Peeker stream, void*, mem·Allocator heap, void*);
error    bio·writenewick(bio·Tree tree, io·Putter out, void*);

// -----------------------------------------------------------------------
// Sequences

/* i/o */
typedef struct bio·FastaReader bio·FastaReader;
typedef struct bio·FastqReader bio·FastqReader;

typedef struct bio·Seq
{
    int  len;
    byte *name;
    byte *s;
    byte *q;
} bio·Seq;

bio·FastaReader *bio·openfasta(io·Reader stream, void *s, mem·Allocator heap, void *h);
error            bio·readfasta(bio·FastaReader *rdr, bio·Seq *seq);
error            bio·closefasta(bio·FastaReader *rdr);

bio·FastqReader *bio·openfastq(io·Reader stream, void *s, mem·Allocator heap, void *h);
error            bio·readfastq(bio·FastqReader *rdr, bio·Seq *seq);
error            bio·closefastq(bio·FastqReader *rdr);

/* alignment */
enum
{
    aln·K = 20,     // kmer size (k <= 32)
    aln·L = 3,      // number of kmers / hash
    aln·N = 1000,   // number of hashes
};

error aln·sketch(byte *seq, int l, uint64 *phi[aln·N], int *locs[aln·N]);
error aln·sort(int len, int l, uint64 *vals);