my $genome_file = @ARGV[0]; open(MULTIFASTA,"<",$genome_file) or die "Can't open $genome_file!\n"; my $header; my $sequence; $/ = "\n"; my %hash; my @fasta = `awk '/^>/{f=++d".fasta"} {print > f}' $genome_file`;