Consensus sequence
 

I compare some 454 reads with a reference genome using blastn with ungapped alignment. I have converted the blast reslts into fasta format. Now I want to convert the reads (output of blast) into consensus sequence. I tried the bioperl module Bio::SimpleAlign but that didn't work. Need some suggestions......