Mercurial > repos > mcharles > rapsosnp
diff rapsodyn/fastaGroomerForMakeBlastdb.pl @ 0:442a7c88b886 draft
Uploaded
| author | mcharles |
|---|---|
| date | Wed, 10 Sep 2014 09:18:15 -0400 |
| parents | |
| children | 3f7b0788a1c4 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/rapsodyn/fastaGroomerForMakeBlastdb.pl Wed Sep 10 09:18:15 2014 -0400 @@ -0,0 +1,22 @@ +#!/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);