Perl & Genomics

#!/usr/bin/perl
 
# Constant: we write here the FILENAME of the multifasta 
# file containing the contigs
$multifasta = '/home/geno/tools/contigs.fna';
 
# We print a MESSAGE for the user
print STDERR "
   Pick Primer Tool v.01
 
   Parameters:
   contig1 direction1 contig2 direction2
 
   direction1,2 = U or C
 
";
 
# We populate variables with USER PARAMETERS
($contig1, $dir1, $contig2, $dir2) = @ARGV;
 
# And print them for the user
print STDERR "
contig1   $contig1
dir1      $dir1
contig2   $contig2
dir2      $dir2
";
open(I, "$multifasta") || die " ERROR: can't open $multifasta :(\n";
 
# Here we read the multifasta file
 
while (<I>) {
	chomp($_);
 
	if ($_=~/>(.*)/) {
 
	   if ($header eq $contig1) {
	   	 	$seq1 = $seq;
	   }
	   if ($header eq $contig2) {
	   		$seq2 = $seq;
	   }
 
	   $header=$1;
	   $seq = '';
	} else {
		$seq=$seq.$_;
	}
}
 
if ($header eq $contig1) {
 	$seq1 = $seq;
}
if ($header eq $contig2) {
	$seq2 = $seq;
}
 
if (length($seq1) == 0) {
	die "$contig1 was not found\n";
}
if (length($seq2) == 0) {
	die "$contig2 was not found\n";
}

#!/usr/bin/perl # Constant: we write here the FILENAME of the multifasta # file containing the contigs $multifasta = '/home/geno/tools/contigs.fna'; # We print a MESSAGE for the user print STDERR " Pick Primer Tool v.01 Parameters: contig1 direction1 contig2 direction2 direction1,2 = U or C "; # We populate variables with USER PARAMETERS ($contig1, $dir1, $contig2, $dir2) = @ARGV; # And print them for the user print STDERR " contig1 $contig1 dir1 $dir1 contig2 $contig2 dir2 $dir2 "; open(I, "$multifasta") || die " ERROR: can't open $multifasta :(\n"; # Here we read the multifasta file while (<I>) { chomp($_); if ($_=~/>(.*)/) { if ($header eq $contig1) { $seq1 = $seq; } if ($header eq $contig2) { $seq2 = $seq; } $header=$1; $seq = ''; } else { $seq=$seq.$_; } } if ($header eq $contig1) { $seq1 = $seq; } if ($header eq $contig2) { $seq2 = $seq; } if (length($seq1) == 0) { die "$contig1 was not found\n"; } if (length($seq2) == 0) { die "$contig2 was not found\n"; }