Mercurial > repos > mcharles > rapsosnp
comparison rapsodyn/fastaGroomerForMakeBlastdb.pl @ 0:442a7c88b886 draft
Uploaded
author | mcharles |
---|---|
date | Wed, 10 Sep 2014 09:18:15 -0400 |
parents | |
children | 3f7b0788a1c4 |
comparison
equal
deleted
inserted
replaced
-1:000000000000 | 0:442a7c88b886 |
---|---|
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); |