wolffd@0: %bnet = mk_motif_hhmm('motif_pattern', 'acca', 'background', 't');
wolffd@0: bnet = mk_motif_hhmm('motif_pattern', 'accaggggga', 'background', []);
wolffd@0: 
wolffd@0: chars = ['a', 'c', 'g', 't'];
wolffd@0: Tmax = 100;
wolffd@0: 
wolffd@0: for seqi=1:5
wolffd@0:   evidence = cell2num(sample_dbn(bnet, 'length', Tmax));
wolffd@0:   chars(evidence(end,:))
wolffd@0: end