Shulou(Shulou.com)06/01 Report--

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.

Welcome to subscribe "Shulou Technology Information " to get latest news, interesting things and hot topics in the IT industry, and controls the hottest and latest Internet news, technology news and IT industry trends.