Write a BioPerl-based script ("pick-by-id.pl") to extract sequences by id. It will emulate the "blastdbcmd -db -entry" command. Use "GBB.pep" and "BBA15" as the two arguments. Usage: ./pick-by-id.pl GBB.pep BBA15. Your script should handle exceptions gracefully, e.g., by printing a warning message if the sequence is not found.
So, this perl one liner code can extract fasta sequences by providing a value for the id:
perl -ne 'if(/^>(\S+)/){$c=grep{/^$1$/}qw(id)}print if $c' fasta.file
The assignment is asking to create a bioperl script that can do the same, extract fasta sequences , by providing the id of the sequence as the second argument and the filename as the first argument on the command line.
Attachment:- Assignment.docx