BuKyung Create a flat text file database of protein sequences with hash function in Perl
From Biolecture.org
Revision as of 05:56, 17 June 2016 by imported>Baik BuKyung
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