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;

創作者介紹

Philippe's Expérience Note

philip 發表在 痞客邦 PIXNET 留言(0) 人氣()