Difference between revisions of "BuKyung Extract a sequence MKKTGIKG from ASMKATAHQMKKTGIKGMSTYALLRL and print it out in a file"

From Biolecture.org
imported>Baik BuKyung
imported>Baik BuKyung
 
Line 62: Line 62:
  
 
<p><span style="font-size:16px">is in the outer.fasta file.</span></p>
 
<p><span style="font-size:16px">is in the outer.fasta file.</span></p>
 
<p>&nbsp;</p>
 
 
<p>&nbsp;</p>
 
 
<p>&nbsp;</p>
 
 
<p>&nbsp;</p>
 
 
<p>&nbsp;</p>
 
 
<p>&nbsp;</p>
 
 
<p><span style="font-size:16px">The content of outer.fasta file&nbsp;is</span></p>
 
 
<p>&nbsp;</p>
 
 
<p>&nbsp;</p>
 
 
<p>&nbsp;</p>
 
 
<p>&nbsp;</p>
 
 
<div>
 
<div><em>&gt;0</em>
 
 
<p><em>ACCACTACTAAGCGCATGAACGACTGTTAGGTTTCCGATGGCTGCTTGCGTTCCGTGTTCCAGCTGACTGGGCTGAACTATTTGTAATGTTGGTTGCACT<br />
 
&nbsp;GC content is:0.49<br />
 
&gt; 1<br />
 
CAGGTACACGGACTGTTTGGTTTGCCCAATTAATTGGCGGGTCGTAAACCGGTTTTTCGTTGGGCGCGGAGTTGTCGTAAACGGTCGGTATTAACTACCT<br />
 
&nbsp;GC content is:0.5<br />
 
&gt; 2<br />
 
ATATTCTGTTCGAAGGCGAGGCCTTAATAAACGGGCTCACACTATACGTTTCTAGCGTGCCAGTACGCGTATGCCCTGAGCAGCATCTTGAATAGTCCTT<br />
 
&nbsp;GC content is:0.48<br />
 
&gt; 3<br />
 
CACGTCTTGAGGCATGCTCACATAACTTGGGATTGATACAATCGGGGGACGGTAGCGGGGCTAGTGGGCATCGTCGGCGGTCTACGAGCAAAAGTATCAG<br />
 
&nbsp;GC content is:0.55<br />
 
&gt; 4<br />
 
CAGGACGTGAACCGAAAGCTGCACACCTATACTATCGTAGTATACCACCGTTCCGTAAATCCATCGCTGATCCTGCCATGAAGGGCTAAGTACGCATGAG</em></p>
 
<em> &nbsp;GC content is:0.5</em></div>
 
 
<div>&nbsp;</div>
 
 
<div>&nbsp;
 
<p><em>Average sequence length is&nbsp;:100<br />
 
&nbsp;GC contents:0.504</em></p>
 
<em>&nbsp;AT contents:0.496</em></div>
 
 
</div>
 
</div>
 
<div>&nbsp;</div>
 
 
<div>&nbsp;</div>
 
 
<div>&nbsp;</div>
 
</div>
 
 
<p><!-- end content --></p>
 
 
<div class="visualClear">&nbsp;</div>
 
 
<div class="portlet" id="p-logo">&nbsp;</div>
 
 
<p><!-- end of the left (by default at least) column --></p>
 
 
<div class="visualClear">&nbsp;</div>
 

Latest revision as of 21:42, 16 June 2016

Back to Baik BuKyung


Source code:


#!/usr/bin/perl

 use strict;
 use warnings;
 open FH, ">", "outer.fasta" or die "$!\n";
 
 my $Original_sequence="ASMKATAHQMKKTGIKGMSTYALLRL";
 my $Extracted_sequence="MKKTGIKG";
 if($Original_sequence =~ /$Extracted_sequence/){
  print FH ("Original sequence is :",$Original_sequence,"\nExtracted sequence is :",$Extracted_sequence);
}
 else{
  print FH ($Extracted_sequence, " is not in the ", $Original_sequence);

 }

 


 

Result

 

 

After the 7.pl is executed, the outer.fasta file is created.

 

In outer.fasta file,

 

Original sequence is :ASMKATAHQMKKTGIKGMSTYALLRL
Extracted sequence is :MKKTGIKG

 

If MKKTGIKG is not in ASMKATAHQMKKTGIKGMSTYALLRL,

 

MKKTGIKG is not in the ASMKATAHQMKKTGIKGMSTYALLRL

 

is in the outer.fasta file.