[Bioperl-l] Count or weight matrix in bioperl?
Torsten Seemann
torsten.seemann at infotech.monash.edu.au
Thu Feb 9 02:54:39 UTC 2006
> Say I have an array of nucleotide sequences of of length N. I want to calculate the count matrix (weight matrix). That is for each position 1..N, I want to know how many As, Cs ,Ts and Gs there are. Is the code to do this already written in bioperl to build this matrix if I pass it those strings?
> Please excuse my lack of knowledge as I am a new comer to bioinformatics.
Use the Bio::Tools::SeqStats module. The PDoc documentation even has an
example similar to what you want to do:
http://doc.bioperl.org/releases/bioperl-1.5.0-RC1/Bio/Tools/SeqStats.html
--Torsten Seemann
More information about the Bioperl-l
mailing list