Difference between revisions of "BuKyung Randomly generate five 100 AA long protein sequences and store them in a FASTA file"
imported>Baik BuKyung |
imported>Baik BuKyung |
||
Line 5: | Line 5: | ||
<hr /> | <hr /> | ||
+ | <div> | ||
<div> | <div> | ||
<div> | <div> | ||
Line 11: | Line 12: | ||
use warnings;<br /> | use warnings;<br /> | ||
open FH, ">", "outer.fasta" or die "$!\n";<br /> | open FH, ">", "outer.fasta" or die "$!\n";<br /> | ||
− | my @seq=("A", " | + | my @seq=("K", "N", "I", "M", "R", "S", "Y", "F", "C", "W", "E", "D", "V", "G", "A", "Q", "H", "L", "R", "P");</div> |
<div> </div> | <div> </div> | ||
Line 18: | Line 19: | ||
print FH ">", $i,"\n";<br /> | print FH ">", $i,"\n";<br /> | ||
for(my $j=0; $j<100;$j++){<br /> | for(my $j=0; $j<100;$j++){<br /> | ||
− | print FH $seq[int(rand( | + | print FH $seq[int(rand(20))];<br /> |
}<br /> | }<br /> | ||
print FH "\n"<br /> | print FH "\n"<br /> | ||
Line 26: | Line 27: | ||
<div>close FH;</div> | <div>close FH;</div> | ||
+ | </div> | ||
</div> | </div> | ||
</div> | </div> | ||
Line 33: | Line 35: | ||
<div> | <div> | ||
<hr /> | <hr /> | ||
− | <p><img alt="" src="/ckfinder/userfiles/images/%EC%BA%A1%EC%B2%9814.PNG" style="height: | + | <p><img alt="" src="/ckfinder/userfiles/images/%EC%BA%A1%EC%B2%9814(2).PNG" style="height:270px; width:800px" /></p> |
</div> | </div> | ||
Line 49: | Line 51: | ||
<p><em>>0<br /> | <p><em>>0<br /> | ||
− | + | GQHRCEMGHIARAMCGDRRFIMHMGIAAICHWCRRIFYVFKPGRGQFWEAVMAKHDWYGVWCGHCNCHMPKMHCPPKCCPIVHEYCHQFNPMIYWYSMQE<br /> | |
>1<br /> | >1<br /> | ||
− | + | FRDPIDEPPIYHWIWCFYHYDACSQHHKPICPVMVDMDQKHHVGRLSQRVLPNKKWIGAREIFHQDACKELHLPNEELFKSSKARMWRAERAGHANKKFS<br /> | |
>2<br /> | >2<br /> | ||
− | + | RHGGRGERWWMIAGFSVWHVVHQKIYIIASWKPRPRGSSWFSDRNQKVYNGGQYSLKRFRMEMQKLKRWGDCIPYGRDANGGQGKRNIKCMESWLHHAIK<br /> | |
>3<br /> | >3<br /> | ||
− | + | ALCRDRLLGLDLCVNIPMENMPHEQWLDRDYYRQHRKQYVEGFVNIIYRIMYHGREWHRREFVDARCADRHVYLYQKAKHRENDDGMGWFHQDVPCNLKL<br /> | |
>4<br /> | >4<br /> | ||
− | + | PDICDCCGWWELFSFGSIRHRRVCMVESGEPACFAKAIRGVINPDLNRLFAGAGLYYASYVDLCFLLYMDYNNGGVHYKYVYISRYVNMKAIKIAIEIFG</em><br /> | |
</p> | </p> | ||
− | |||
− | |||
</div> | </div> |
Latest revision as of 21:11, 16 June 2016
Back to Baik BuKyung
Source code:
use strict;
use warnings;
open FH, ">", "outer.fasta" or die "$!\n";
print FH ">", $i,"\n";
for(my $j=0; $j<100;$j++){
print FH $seq[int(rand(20))];
}
print FH "\n"
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