Perl subroutine to read genome/reads fasta file !

Rahul Nayak — Tue, 03 Jul 2018 04:56:15 -0500

sub readSeqFromFasta{
	my $file = $_[0];
	my (%map,$id,$seq);
	open(IN,$file) or die("Cannot open file for reading $file:$!\n");
	while(){
		chomp;
		if($_ =~ /^>/){
			if(defined $id && defined $seq){
				$map{$id} = $seq;
				$seq="";
			}
			($id) = split(' ',$_);
			$id=~s/>//;
		}
		else {
			$seq=$seq.$_;
		}
	}
	$map{$id} = $seq;
	close(IN);
	return \%map;

}

BOL: Perl subroutine to read genome/reads fasta file !

Perl subroutine to read genome/reads fasta file !