Generate consensus from MultiFasta

[Antony03]

Hi,

I have a multifasta like:

>sequence1
ATCGATCATCG
>sequence2
ATCGATNNTCG
>sequence3
ATCGATCATCG

All sequences are the same, except there are some N. I want to make a consensus with the N if there are.

Do you know I can do this?
Thanks!

[Antony03]

A little up!

Please

[gringer]

here's my rough 5-minute attempt:

#!/usr/bin/perl
use warnings;
use strict;
my @sequence = ();
while(<>){
  if(!/^>/){
    chomp;
    my @bases = split(//,$_);
    for(my $i = 0; $i < @bases; $i++){
      if($bases[$i] eq "N"){
        $sequence[$i] = "N";
      } elsif($sequence[$i] ne "N") {
        $sequence[$i] = $bases[$i];
      }
    }
  }
}

print(">consensus\n".join("",@sequence)."\n");

[may not compile, will probably contain bugs, will probably be very slow]