F0Z7V0F01A03EB_210
F0Z7V0F01A03EB_180
F0Z7V0F01A03EB_136
F0Z7V0F01A03EB_362
F0Z7V0F01A03GP_298
F0Z7V0F01A03GP_205
etc...........
####
>GJKKTUG01DYDGC
GGGTATTCCTTCTCCACCTTGCAGCTAACATCAGTGTTTCGTCTACTCAAGCACGCCAAC
ACGCCCTAGAGCGCCCTGTCCAGGGGATGGCAACCAACTCTGACCCTGCAAGTGCAGCAG
ACATGAGGAATACAAACTACAATCTTTTACTTGATGATGCAATGCCGGACAAACTCTAGA
>F0Z7V0F01EDB3V
AAGGCGAGNGGTATCACGCAGTAAGTTACGGTTTTCGGGTAACGCGTCNGNGGNACTAAC
CCACGGNGGGTAACCCGTCNCTACCGGTATAGGACTAAGGTTACCGGAACGTCGTGGGGT
ACCCCCCGGACGGGGACCGTCCCCTCATANAGTCAACNGTNTGAGATGGACTAACTCAAA
CCTAGTTTCAAGTACTATTTAACTTACTTACGTTACCCGTAATTTCGGCGTTTAGAGGCG
etc....................
##
##
#!/usr/bin/perl
use strict;
use warnings;
use Bio::SeqIO;
use Data::Dumper;
my %seq_id;
my $fasta_id;
open HEADER , "){
chomp $_;
#print $_."\n";
$fasta_id = $_;
$fasta_id =~ s/_.*//g ;
#print $fasta_id."\n";
%seq_id = ("$fasta_id" => "$_");
# print Dumper (\%seq_id);
my $infile=$ARGV[0] || die ("Please give me an input fasta file\n");
my $inseq = new Bio::SeqIO(-format => 'fasta',
-file => $infile);
while (my $seq_obj = $inseq->next_seq ) {
my $id = $seq_obj->id ;
chomp $id;
# print $id."\n";
my $seq = $seq_obj->seq ;
if (exists ($seq_id{$id})) {
print ">";
print $seq_id{$fasta_id};
print "\n".$seq."\n";
}
}
}