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

From Biolecture.org
Revision as of 18:10, 12 June 2016 by imported>Baik BuKyung

Back to Baik BuKyung


Source code:


#!/usr/bin/perl

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

my @seq=("A", "T", "G", "C");
 
for(my $i=0; $i<5;$i++){

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

}
 
close FH;
 


Result

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

The content is

 

>0
AGAGCATAATATTGGGGCGGACTTGAAATCCAGTTAATTGGGGGTTCTGGAGCATTGGTCTGGTCGCCTGCTACGCTTATCACTACTGGATGAAACTTCC
>1
ACGGCCAAGGTGTGTGTCATAGCTGTTCCGTCGTCTATCCTGAACGAGCCTATAAATACATGAACCTTGGCAGGTTACAAACGCCGGAAGCGTTCTGGGC
>2
CACGGCCGCCCACATGGGACAACCTGGCATCGGCGGTAAAGGCTCGTACCATAGTTGAAGCGAGACGGCAGTATAAGCACCGTCCGAATCGTCACCCTTT
>3
AGAACGGATAGGATCGTTCCCGTGATACCCGTAGCCCAGATCAAAAAAATACCACACAATATAAATGTTACAAAGCGGACGCCTAATAACCACCTAGGAG
>4
TCTTTCCGGCTCGACTTATCATCGGTGGATCGAATAATCAATAGGGCTTCGAACTATAGAGGGTTTCAAGGTCTGACTCTAGAGGACGTCTGGGGAGGGC