Beefy Boxes and Bandwidth Generously Provided by pair Networks
more useful options

Re: PROTEIN FILE help me pleaseee

by pvaldes (Chaplain)
on Jan 30, 2013 at 15:13 UTC ( #1016081=note: print w/replies, xml ) Need Help??

in reply to PROTEIN FILE help me pleaseee

Deliberately incomplete, after all this is (your) homework, but you have some of your problems solved here... only because Oxytricha is a cute green thing

#!/usr/bin/perl -w use strict; my @protein = (); my $phenylcounter = 0; open (my $PROTEFILE, '<', $ARGV[0]) or die $!; # If the file is not existing it gives an error message # yup, I avoid the <STDIN> idea, solved by you. while ($line = <$PROTEFILE>) { # We read the lines of the file next if $line =~ m/^>/; # If the line starts with ">" + is not considered. $line =~ s/\s//g; # all white spaces removed from the line # Using the translate command, the program counts the number of F in +the sequence, assigning it to a variable. # left for you... tr/F//... $phenylcounter++} } # the while loop terminates close $PROTEFILE; # and the program continues closing the file. print "The aminoacid sequence", $ARGV[0], " contains ", $phenylcounter +, " Phenylalanine aminoacids"; print "\n"; # followed by new line

Updated: ($phenilcounter != $phenylcounter), fixed now

Log In?

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://1016081]
and all is quiet...

How do I use this? | Other CB clients
Other Users?
Others examining the Monastery: (9)
As of 2018-05-22 14:17 GMT
Find Nodes?
    Voting Booth?