Philippe's Expérience Note

Name:
CutPeptide.pl

Description:
將ensemble的 peptide sequences (fasta files),個別輸出成單筆sequence

Source Code:

#!/usr/bin/perl
##############
# Edit by Philippe
# 05.14.08
##############
# input FASTA files
open(IN, "/home/philippe/CutPeptide/Homo_sapiens.NCBI36.49.pep.all.fas");
while (line = )
{
        chomp(line);
        # regular expression, match >
        if (line = />(.*)/)
        {
                ENSP = 1;
                print "ENSP\n";
                # output FASTA files
               open (OUT, "+>/home/philippe/CutPeptide/pep_all_fasta/ENSP.fas");
                print OUT ">ENSP\n";
        }
        else
        {
                print OUT "line\n";
        }
}
close IN;
close OUT;