function wtseq = infer_wildtype (data)
% wtseq = infer_wildtype (data)
%
% Infers the wildtype sequence from the data.

nfdata = make_nfdata (data);

l = size (data.seqs, 2);

counts_A = sum (data.seqs == 'A');
counts_C = sum (data.seqs == 'C');
counts_G = sum (data.seqs == 'G');
counts_T = sum (data.seqs == 'T');

counts = [counts_A ; counts_C ; counts_G ; counts_T];

[~, nwtseq] = max (counts, [], 1);

wtseq = nseq2seq (nwtseq);

end