Beefy Boxes and Bandwidth Generously Provided by pair Networks
Keep It Simple, Stupid

Re: Counting amino acids

by space_monk (Chaplain)
on Jul 19, 2013 at 05:53 UTC ( #1045259=note: print w/replies, xml ) Need Help??

in reply to Counting amino acids

As the previous poster has said, you have not given us much to go on, but if you are just starting out and have nothing so far, here's a plan of attack for the problem

One way to do it is to write your own routine which would read FASTA format and then simply greps or counts the number of amino acids. This is not the recommended way however!

A better way is to look on CPAN and see if someone has done some or all of the work for you.Bio::DB:Fasta and Bio::SeqIO, FastaParse and Bio::Phylo::IO are all starting points for reading FASTA files and parsing FASTA data.

If you spot any bugs in my solutions, it's because I've deliberately left them in as an exercise for the reader! :-)

Log In?

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

How do I use this? | Other CB clients
Other Users?
Others perusing the Monastery: (3)
As of 2018-02-24 18:34 GMT
Find Nodes?
    Voting Booth?
    When it is dark outside I am happiest to see ...

    Results (310 votes). Check out past polls.