- or download this
{
open (IN, $file) || die "Can't open $file\n";
my $tag; my $seq;
...
$protein{$tag} = uc($seq); # add this line to get the last sequence
close IN;
}
- or download this
>chr1 RecName: Full=Putative transcription factor 001R;
MAFSAEDVLKEYDRRRRMEALLLSLYYPNDRKLLDYKEWSPPRVQVECPKAPVEWNNPPS
EKGLIVGHFSGIKYKGEKAQASEVDVNKMCCWVSKFKDAMRRYQGIQTCKIPGKVLSDLD
...
KCINRASDPVYQKVKTLHAYPDQCWYVPCAADVGELKMGTQRDTPTNCPTQVCQIVFNML
DDGSVTMDDVKNTINCDFSKYVPPPPPPKPTPPTPPTPPTPPTPPTPPTPPTPRPVHNRK
VMFFVAGAVLVAILISTVRW
- or download this
>chr1 RecName: Full=Putative transcription factor 001R; 9.425903320
+3125
>chr2 RecName: Full=Uncharacterized protein 002L; 8.2193603515625
- or download this
chr1 9.793558
chr2 7.791410
- or download this
#!/usr/bin/perl
use strict;
use warnings;
...
my $iep = $calc->iep;
print join("\t", $seq->id, $iep), "\n";
}