BuKyung Randomly generate five 100 AA long protein sequences and store them in a FASTA file

Back to Baik BuKyung

---

Source code:

---

#!/usr/bin/perl

 use strict;
 use warnings;
 open FH, ">", "outer.fasta" or die "$!\n";

my @seq=("K", "N", "I", "M", "R", "S", "Y", "F", "C", "W", "E", "D", "V", "G", "A", "Q", "H", "L", "R", "P");

 

for(my $i=0; $i<5;$i++){

 print FH ">", $i,"\n";
 for(my $j=0; $j<100;$j++){
  print FH $seq[int(rand(20))];
 }
 print FH "\n"

}

 

close FH;

 

---

---

Result

After the 3.pl is executed, the outer.fasta file is generated.

The content is

 

>0
GQHRCEMGHIARAMCGDRRFIMHMGIAAICHWCRRIFYVFKPGRGQFWEAVMAKHDWYGVWCGHCNCHMPKMHCPPKCCPIVHEYCHQFNPMIYWYSMQE
>1
FRDPIDEPPIYHWIWCFYHYDACSQHHKPICPVMVDMDQKHHVGRLSQRVLPNKKWIGAREIFHQDACKELHLPNEELFKSSKARMWRAERAGHANKKFS
>2
RHGGRGERWWMIAGFSVWHVVHQKIYIIASWKPRPRGSSWFSDRNQKVYNGGQYSLKRFRMEMQKLKRWGDCIPYGRDANGGQGKRNIKCMESWLHHAIK
>3
ALCRDRLLGLDLCVNIPMENMPHEQWLDRDYYRQHRKQYVEGFVNIIYRIMYHGREWHRREFVDARCADRHVYLYQKAKHRENDDGMGWFHQDVPCNLKL
>4
PDICDCCGWWELFSFGSIRHRRVCMVESGEPACFAKAIRGVINPDLNRLFAGAGLYYASYVDLCFLLYMDYNNGGVHYKYVYISRYVNMKAIKIAIEIFG