BuKyung Extract a sequence MKKTGIKG from ASMKATAHQMKKTGIKGMSTYALLRL and print it out in a file

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.

 

 

 

 

 

 

The content of outer.fasta file is

 

 

 

 

>0

ACCACTACTAAGCGCATGAACGACTGTTAGGTTTCCGATGGCTGCTTGCGTTCCGTGTTCCAGCTGACTGGGCTGAACTATTTGTAATGTTGGTTGCACT
 GC content is:0.49
> 1
CAGGTACACGGACTGTTTGGTTTGCCCAATTAATTGGCGGGTCGTAAACCGGTTTTTCGTTGGGCGCGGAGTTGTCGTAAACGGTCGGTATTAACTACCT
 GC content is:0.5
> 2
ATATTCTGTTCGAAGGCGAGGCCTTAATAAACGGGCTCACACTATACGTTTCTAGCGTGCCAGTACGCGTATGCCCTGAGCAGCATCTTGAATAGTCCTT
 GC content is:0.48
> 3
CACGTCTTGAGGCATGCTCACATAACTTGGGATTGATACAATCGGGGGACGGTAGCGGGGCTAGTGGGCATCGTCGGCGGTCTACGAGCAAAAGTATCAG
 GC content is:0.55
> 4
CAGGACGTGAACCGAAAGCTGCACACCTATACTATCGTAGTATACCACCGTTCCGTAAATCCATCGCTGATCCTGCCATGAAGGGCTAAGTACGCATGAG

 GC content is:0.5

 

 

Average sequence length is :100
 GC contents:0.504

 AT contents:0.496