############ outer loop################### foreach my $fileSeq(@sequences) # loop through all the sequences { $counter++; $seqobj = ''; print("
$counter
"); $seqobj = Bio::PrimarySeq->new ( -seq => $fileSeq, -alphabet => 'dna'); # create an object of the sequence print("
Sequence object created....all going well so far...
"); ############ inner loop ################### foreach $element(@enzymes) # look through all the enzymes for each sequence { $newObject = ''; $newObject = new Bio::Tools::RestrictionEnzyme(-NAME =>$element); # create an object of the enzyme print("
Enzyme object created for enzyme : $element
"); $cutSeq = $newObject->site(); # show the cut site for the enzyme $anotate = $newObject->annotate_seq($seqobj); # annotate the sequence at cut site print("
Sequence virtually labelled....
"); @cut = $newObject->cut_seq($seqobj); # virtualy cut the sequnce using the enzyme print("
Sequence virtually cut....
"); print(""); # produce output print(""); print(""); print(""); print("
Sequence number $counter
"); print("The cleavage site is $cutSeq for enzyme $element denoted by \" ^ \"
The annotated sequence is given below along with the digested sequence
"); print("Annotated:$anotate
"); print("Digested:@cut


"); } ############ END inner loop ################### print("Loop finished... returning to the start
"); }############ END outer loop ###################