View source

From NMR Wiki

for One2three - perl script. One letter to three letter protein sequences

You do not have permission to do that, for the following reason:

The action you have requested is limited to users in the group user.

You can view and copy the source of this page:

Save this file as one2three.

This script can be useful to prepare input for [[QUEEN]] software.

Usage:
 one2three file > output_file

==Script==
<source lang="perl">
#!/usr/bin/perl

my %AA = (ALA=>'A',TYR=>'Y',MET=>'M',LEU=>'L',CYS=>'C',GLY=>'G',
         ARG=>'R',ASN=>'N',ASP=>'D',GLN=>'Q',GLU=>'E',HIS=>'H',TRP=>'W',
         LYS=>'K',PHE=>'F',PRO=>'P',SER=>'S',THR=>'T',ILE=>'I',VAL=>'V');

my %aa = reverse %AA;

if (scalar @ARGV == 0)
{
        print "\n\tconvert one-letter aminoacid sequence file to three-letter with one residue per line\n";
        use File::Basename;
        print "\t",basename($0), " <one-letter seq file>\n";
        print "\tempty spaces and multiline files are allowed\n\n";
        exit 1;
}

my $file = $ARGV[0];
open F, "<$file" or die $!;
my @lines = <F>;
chomp foreach @lines;

my $let = join('',@lines);
$let =~ s/\s//g;
my @aa = split /|/,$let;

print ">A\n"; #use dummy chain name
foreach my $a (@aa)
{
        print $aa{$a}, "\n";
}
</source>

[[Category:Software]]
[[Category:Scripts]]
[[Category:Perl]]

Return to One2three - perl script. One letter to three letter protein sequences.

Retrieved from "http://nmrwiki.org/wiki/index.php?title=One2three_-_perl_script._One_letter_to_three_letter_protein_sequences"

View source

From NMR Wiki

Views

Personal tools

navigation

community

Contents

Toolbox