silver_steve has asked for the
wisdom of the Perl Monks concerning the following question:
I have a .qual file, and I would like to find out the average score for each position across all sequences. For instance, if I have 10 sequences, each 100 bp, what is the average qual score of the values at each index across all sequences?
My strategy is to split the scores by index into their own separate arrays, and then take the average of each array.
I know there may be modules in Bio::Seq that can help with this, but I am not familiar enough with those yet to use them effectively.
Does anyone know if it is necessary to use a module, and if so, are there any existing ones that can help with this?