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