Open main menu

Biolecture.org β

Changes

no edit summary
<h3p><span style="font- Ideasize:14px"><strong>This script use&nbsp;amino acids as&nbsp;random function characters, and display 100aa long sequence in a fasta file:&nbsp;</strong></span></h3p>
<p>&nbsp#!/usr/bon/perl<br />use strict;<br />use warnings; * 20개의 amino acid를 랜덤하게 나타나게 하는 작업을 100번 반복</p>
<p><strong>&nbsp; * open (FASTA, &quot;&gt;random_5AA_seq.fasta&quot;) or die;</strong><br />&nbsp;<strong>my @range = qw(A G T F P R K L I H C S M Y W D E N Q V)</strong>;<br />&nbsp;my $seq;<br /><strong>&nbsp;for my $i (1..5)</strong> {<br />&nbsp; print FASTA $i,&quot;. &quot;;<br />&nbsp;<strong>for (1..100)</strong> {<br />&nbsp; <strong>$seq = $range[rand @range]</strong>;<br />&nbsp;print FASTA $seq; }<br />&nbsp;코드를 5번 반복print FASTA &quot;\n&quot;;}</p>
<p>&nbsp; * 반복작업은 for loop사용------------------------------------</p>
<p>&nbsp; * rand를 이용해서 chars안의 요소들을 랜덤하게 배치Output FASTA file content:</p>
<p>&nbsp;</p> <h3>- Code</h3> <hr /><p>#!/usr/bin/perl<br />&nbsp; use strict;<br />&nbsp; use warnings;<br />&nbsp;&nbsp;<br />&nbsp; my @chars=(&#39;A&#39;,&#39;C&#39;,&#39;D&#39;,&#39;E&#39;,&#39;F&#39;,&#39;G&#39;,&#39;H&#39;,&#39;I&#39;,&#39;K&#39;,&#39;L&#39;,&#39;M&#39;,&#39;N&#39;,&#39;P&#39;,&#39;Q&#39;,&#39;R&#39;,&#39;S&#39;,&#39;T&#39;,&#39;V&#39;,&#39;W&#39;,&#39;Y&#39;);<br />&nbsp; my $random_string;<br />&nbsp; my $i;<br />&nbsp; for($i = 1; $i &lt;=5 ; $i++). ETLKKEINKDDYRRREHVQLKHAFWVPCGPVTVDWGPLGEDDQGLDSRHFHTLQKYTIYFEMPLINPKFATEMTDSDTMYWAVPKGGYEPNYCRSNQLIE<br />&nbsp; {2. ISNWQWIQRRKTCHSASIYYSYTDWFPHDYKWSQCMCDTRLCTPNDEPMAPTGQLCYDFLAGPEVTMEHGRGAKGHGLFHWRGTTTQFYEGAQMYSNGHL<br />&nbsp; &nbsp; print $i,&quot;3.&quot;;ATIYTKEETFVMLWIQGEYIKMLPGSCCPYFPIAYSFYVRDRHDEVWYITQGKTMNTYPILGCYVAAPSQYTMQKAVAMKRPKNTNRVDDMHMRFNLWAW<br />&nbsp; for(14..100)HEFQEVLFARDFEFKRACTSHHMDEIADQIHQWNRALMCVHEESKDCVPGNCCAQCLPGRTYYREEVPCYLYNGHQSKHHYRGRGQEYLKEWKPPMYENY<br />&nbsp; {<br />&nbsp; &nbsp; $random_string = $chars[rand @chars];<br />&nbsp; &nbsp; print $random_string;<br />&nbsp; };5. DACITFMPMKQHTNAENIPNHTCKCQMAHSYNYMEMKHVVQIQRHHPKLIRCMEKHRKSQLVDARNHVKMVGTQLVIEMVDNYDWRSYWGQLPAHYHVST<br />&nbsp;print &quot;\n\n&quot;;<br />&nbsp;};</p> <hr /><h3>- Result</h3> <p><img alt="" src="/ckfinder/userfiles/images/random_result.png" style="height:229px; width:641px" /></p> <h3>&nbsp;</h3> <h3>- Application</h3> <p>&nbsp;* 키보드로부터 입력을 이용해 amino acid 의 갯수와 반복횟수를 입력받아 원하는 sequence의 갯수와 sequence의 길이를 출력한다.&nbsp;</p> <p>&nbsp;* 키보드로부터 입력을 받을때는 &lt;STDIN&gt;을 이용한다.</p>
<p>&nbsp;</p>
 
<h3>-Code_application</h3>
 
<hr />
<p>#!/usr/bin/perl<br />
&nbsp; use strict;<br />
&nbsp; use warnings;<br />
&nbsp;&nbsp;<br />
&nbsp; my @chars=(&#39;A&#39;,&#39;C&#39;,&#39;D&#39;,&#39;E&#39;,&#39;F&#39;,&#39;G&#39;,&#39;H&#39;,&#39;I&#39;,&#39;K&#39;,&#39;L&#39;,&#39;M&#39;,&#39;N&#39;,&#39;P&#39;,&#39;Q&#39;,&#39;R&#39;,&#39;S&#39;,&#39;T&#39;,&#39;V&#39;,&#39;W&#39;,&#39;Y&#39;);<br />
&nbsp; my $random_string;<br />
&nbsp; my $i;<br />
&nbsp; my $repeat;<br />
&nbsp; my $string_number;</p>
 
<p>print &quot;Put in the number of times to repeat : &quot;;<br />
<strong>$repeat = &lt;STDIN&gt;;</strong></p>
 
<p>print &quot;Put in the number of sequence : &quot;;<br />
<strong>$string_number = &lt;STDIN&gt;;</strong></p>
 
<p>&nbsp; for($i = 1; $i &lt;=<strong> $repeat</strong> ; $i++)<br />
&nbsp; {<br />
&nbsp; &nbsp; print $i,&quot;.&quot;;<br />
&nbsp; for(1..<strong>$string_number</strong>)<br />
&nbsp; {<br />
&nbsp; &nbsp; $random_string = $chars[rand @chars];<br />
&nbsp; &nbsp; print $random_string;<br />
&nbsp; };<br />
&nbsp;print &quot;\n\n&quot;;<br />
&nbsp;};</p>
 
<hr />
<p>&nbsp;</p>
 
<p>&nbsp;</p>
 
<p>[[Yeonsong - bioinformatics]]</p>
Anonymous user