Difference between revisions of "BuKyung Randomly generate five 100 AA long protein sequences and store them in a FASTA file"

From Biolecture.org
imported>Baik BuKyung
imported>Baik BuKyung
Line 48: Line 48:
 
<p>&nbsp;</p>
 
<p>&nbsp;</p>
  
<p>&gt;0<br />
+
<p><em>&gt;0<br />
 
AGAGCATAATATTGGGGCGGACTTGAAATCCAGTTAATTGGGGGTTCTGGAGCATTGGTCTGGTCGCCTGCTACGCTTATCACTACTGGATGAAACTTCC<br />
 
AGAGCATAATATTGGGGCGGACTTGAAATCCAGTTAATTGGGGGTTCTGGAGCATTGGTCTGGTCGCCTGCTACGCTTATCACTACTGGATGAAACTTCC<br />
 
&gt;1<br />
 
&gt;1<br />
Line 57: Line 57:
 
AGAACGGATAGGATCGTTCCCGTGATACCCGTAGCCCAGATCAAAAAAATACCACACAATATAAATGTTACAAAGCGGACGCCTAATAACCACCTAGGAG<br />
 
AGAACGGATAGGATCGTTCCCGTGATACCCGTAGCCCAGATCAAAAAAATACCACACAATATAAATGTTACAAAGCGGACGCCTAATAACCACCTAGGAG<br />
 
&gt;4<br />
 
&gt;4<br />
TCTTTCCGGCTCGACTTATCATCGGTGGATCGAATAATCAATAGGGCTTCGAACTATAGAGGGTTTCAAGGTCTGACTCTAGAGGACGTCTGGGGAGGGC<br />
+
TCTTTCCGGCTCGACTTATCATCGGTGGATCGAATAATCAATAGGGCTTCGAACTATAGAGGGTTTCAAGGTCTGACTCTAGAGGACGTCTGGGGAGGGC</em><br />
 
&nbsp;</p>
 
&nbsp;</p>
  
 
<p>&nbsp;</p>
 
<p>&nbsp;</p>
 
</div>
 
</div>

Revision as of 18:10, 12 June 2016

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