Range of quality of base calls at each position in my alignment of 454 reads

[trasver]

Hi all

I have an alignment of 454 reads to a template sequence. I also have the quality file.

What I'd like to do is identify the range of Phred quality scores at each nucleotide position in my alignment using the information contained in the quality file.

Does anyone know of any software out there that might do this?

(bearing in mind that there will be gaps inserted in the reads as well as some stripped nucleotides.)

Thanks for any suggestions

[flxlex]

The quality file will only give you one quality value. I am not sure if there are programs that will give you the range for each position.