(*search using a Timed Search*)

read("9.fas","31.ss")
read(prealigned:(aminoacids:("41_aln.aa"),tcm:"s1t2.mat"))
set(root:"t1")
search (max_time:00:12:00)
select(unique)
report("Run1a.tre",trees)
search(max_time:00:12:00)
select(unique)
report("Run1b.tre",trees)
fuse(iterations:250)
select()
swap(trees:100)
select()
report("Run1c_H86.tre",trees:(hennig,total),"Run1c_cs.tre",consensus,"Run1c.pdf",graphconsensus)
quit()