Mercurial > repos > mcharles > rapsosnp
annotate rapsodyn/fastaGroomerForMakeBlastdb.pl @ 2:761fecc07fa9 draft
Uploaded
| author | mcharles |
|---|---|
| date | Thu, 11 Sep 2014 03:10:47 -0400 |
| parents | 442a7c88b886 |
| children | 3f7b0788a1c4 |
| rev | line source |
|---|---|
| 0 | 1 #!/usr/bin/perl |
| 2 my $inputfasta = $ARGV[0]; | |
| 3 | |
| 4 open(IB, $inputfasta) or die ("Can't open $inputfasta \n"); | |
| 5 | |
| 6 while (my $ligne = <IB>){ | |
| 7 if ($ligne=~/\[.*?\=.*?\]/){ | |
| 8 $ligne =~ s/[\[\]]//g; | |
| 9 print $ligne; | |
| 10 } | |
| 11 elsif ($ligne =~/^\>/){ | |
| 12 print $ligne; | |
| 13 $ligne = <IB>; | |
| 14 $ligne =~ s/N/a/g; | |
| 15 print $ligne; | |
| 16 } | |
| 17 else { | |
| 18 print $ligne; | |
| 19 } | |
| 20 } | |
| 21 | |
| 22 close (IB); |