my @tables = (); my $counter = 0; foreach my $fileSeq (@sequences) { # loop through all the sequences my $seqobj = Bio::PrimarySeq->new ( -seq => $fileSeq, -alphabet => 'dna'); # create an object of the sequence foreach my $element (@enzymes) { # look through all the enzymes for each sequence my $newEnzyme = new Bio::Tools::RestrictionEnzyme(-NAME => $element); push @tables,{ COUNTER => $counter++, CUTSEQ => $newEnzyme->site(), ELEMENT => $element, ANNOTATE => $newEnzyme->annotate_seq($seqobj), CUT => $newEnzyme->cut_seq($seqobj), ##### ???? ###### }; } } .. and in .tmpl file
Sequence number
The cleavage site is for enzyme denoted by " ^ "
The annotated sequence is given below along with the digested sequence
Annotated:
Digested: