############ 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("
| 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 |