I am interested on the frag.dat function on seqRFLP, but this function only work for single enzyme per run (enzymes = X).
frag.dat(fil = fas, enznames = X, enzdat