8) Generating five 100 AA long protein sequences and storing them in a FASTA file
This script use amino acids as random function characters, and display 100aa long sequence in a fasta file:
#!/usr/bon/perl
use strict;
use warnings;
open (FASTA, ">random_5AA_seq.fasta") or die;
my @range = qw(A G T F P R K L I H C S M Y W D E N Q V);
my $seq;
for my $i (1..5) {
print FASTA $i,". ";
for (1..100) {
$seq = $range[rand @range];
print FASTA $seq; }
print FASTA "\n";}
------------------------------------
Output FASTA file content:
1. ETLKKEINKDDYRRREHVQLKHAFWVPCGPVTVDWGPLGEDDQGLDSRHFHTLQKYTIYFEMPLINPKFATEMTDSDTMYWAVPKGGYEPNYCRSNQLIE
2. ISNWQWIQRRKTCHSASIYYSYTDWFPHDYKWSQCMCDTRLCTPNDEPMAPTGQLCYDFLAGPEVTMEHGRGAKGHGLFHWRGTTTQFYEGAQMYSNGHL
3. ATIYTKEETFVMLWIQGEYIKMLPGSCCPYFPIAYSFYVRDRHDEVWYITQGKTMNTYPILGCYVAAPSQYTMQKAVAMKRPKNTNRVDDMHMRFNLWAW
4. HEFQEVLFARDFEFKRACTSHHMDEIADQIHQWNRALMCVHEESKDCVPGNCCAQCLPGRTYYREEVPCYLYNGHQSKHHYRGRGQEYLKEWKPPMYENY
5. DACITFMPMKQHTNAENIPNHTCKCQMAHSYNYMEMKHVVQIQRHHPKLIRCMEKHRKSQLVDARNHVKMVGTQLVIEMVDNYDWRSYWGQLPAHYHVST