Rewrite the "extract.pl" using BioPerl, including the use of Bio::SeqIO to read the genome FASTA file ("mystery_seq1.fas") and the use of Bio::Seq for obtaining coding sequences and translating sequences. Your output of protein sequences should not contain stop codons except as the last codon.
Attached is the mystery_seq1.coord, mystery_seq1.coord and the extract.pl files. The extract.pl code uses the mystery_seq1.coord file to reference the open reading frame sequences in the mystery_fas1.coord fasta file. The command should be as follows (in my directory as an example):./extract.pl ../../bio425/data/mystery_seq1.fas mystery_seq1.coord
Attachment:- ASSIGNMENT.rar