BuKyung Extract a sequence MKKTGIKG from ASMKATAHQMKKTGIKGMSTYALLRL and print it out in a file
Back to Baik BuKyung
Source code:
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
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
Average sequence length is :100
GC contents:0.504