(* Chromosome analysis of unannotated sequences *)
read(chromosome:("ua15.fas"))
transform((all,dynamic_pam:(locus_breakpoint:20,locus_indel:
(10,1.5),circular:true,median:2, min_seed_length:15, 
min_rearrangement_len:45, min_loci_len:50,median:2,swap_med:1)))
build() 
swap()
select()
report("chrom",diagnosis)
report("consensustree",graphconsensus)
exit()

