How to calculate the Molecular weight of protein in batch by perl

In this article, the editor introduces in detail "perl how to batch calculate protein molecular weight", the content is detailed, the steps are clear, and the details are handled properly. I hope this article "perl how to batch calculate protein molecular weight" can help you solve your doubts.

The perl script is as follows:

Die "perl $0" unless (@ ARGV==2); use Bio::SeqIO;use Bio::Seq;use Bio::Tools::SeqStats;use Bio::Tools::pICalculator;use Data::Dumper;# read sequence my $in = Bio::SeqIO- > new (- file = > "$ARGV [0]",-format = > 'Fasta'); open OUT, "> $ARGV [1]" or die "$!"; print OUT "# ID\ tlength\ tMV (Da)\ tpI\ n" My $calc = Bio::Tools::pICalculator- > new (- places = > 2in- Kset = > 'EMBOSS'); # read sequence while (my $seq = $in- > next_seq ()) {my ($id, $sequence, $desc) = ($seq- > id, $seq- > seq, $seq- > desc); my $weight = Bio::Tools::SeqStats-> get_mol_wt ($seq); $calc- > seq ($seq); my $iep = $calc- > iep Print OUT sprintf ("% s\ t% s\ n", $seq- > id, $seq- > length, "$weight- > [0]", $iep);} $in- > close (); close (OUT) After reading this, the article "how to batch calculate the molecular weight of protein in perl" has been introduced. If you want to master the knowledge of this article, you still need to practice and use it to understand it. If you want to know more about related articles, welcome to follow the industry information channel.