view rapsodyn/fastaGroomerForMakeBlastdb.pl @ 15:56d328bce3a7 draft default tip

Uploaded
author mcharles
date Thu, 29 Jan 2015 08:54:06 -0500
parents 3f7b0788a1c4
children
line wrap: on
line source

#!/usr/bin/perl
#V1.0.0
my $inputfasta = $ARGV[0];

open(IB, $inputfasta) or die ("Can't open $inputfasta \n");

while (my $ligne = <IB>){
	if ($ligne=~/\[.*?\=.*?\]/){
		$ligne =~ s/[\[\]]//g;
		print $ligne;
	}
	elsif ($ligne =~/^\>/){
		print $ligne;
		$ligne = <IB>;
		$ligne =~ s/N/a/g;
		print $ligne;
	}
	else {
		print $ligne;
	}
}

close (IB);