Our Sponsors
Perl script to rename the fasta file !
By Abhinav 850 days ago
#Script #1
#!/usr/bin/perl -w
use strict;
#USAGE
#perl extractPattern.pl kmerfasta > uniref100_result_broad
my %kHash;
local $/ = '>';
my $infile2 = "$ARGV[0]"; # Kmer fasta
open( FH2, '<', $infile2 ) or die $!;
while (<FH2>) {
my @allVal = split ('\n', $_);
my $namewa;
for (my $j=0; $j<scalar(@allVal); $j++) {
next if ($allVal[$j] eq ">");
if ($j == 0) {$namewa=$allVal[0]; next;}
my $num=$j;
my @lName=split(' ',$allVal[0]);
#print "$allVal[$j]\t$lName[0]\n";
#$kHash{$allVal[$j]}=$lName[0];
print ">$lName[0].$num\n$allVal[$j]\n";
}
}
