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