Mercurial > repos > mcharles > rapsosnp
view rapsodyn/fastaGroomerForMakeBlastdb.pl @ 4:9074a5104cdd draft
Uploaded
author | mcharles |
---|---|
date | Wed, 17 Sep 2014 04:20:08 -0400 |
parents | 442a7c88b886 |
children | 3f7b0788a1c4 |
line wrap: on
line source
#!/usr/bin/perl my $inputfasta = $ARGV[0]; open(IB, $inputfasta) or die ("Can't open $inputfasta \n"); while (my $ligne = <IB>){ if ($ligne=~/\[.*?\=.*?\]/){ $ligne =~ s/[\[\]]//g; print $ligne; } elsif ($ligne =~/^\>/){ print $ligne; $ligne = <IB>; $ligne =~ s/N/a/g; print $ligne; } else { print $ligne; } } close (IB);