Trim seq

#! /usr/bin/perl -w
# *************************************************************
# This Programme does the following in sequence
# -> Parses a multi fasta
# -> Removes the sequences that has less than 250 bases in size
# -> Removes the sequence if they have >20% of N
# -> Capture the results

# Author: Rajkumar (itc@rajkumar.in)

# *************************************************************

# Libraries
	use strict;
	use Bio::Perl;
	use Bio::SeqIO;

# Variables
	my $path = "/home/raj/bio";
	my $multifasta  = "$path/tobacco/data/singlet.fa";
	
# making several Fasta files of GDNA
	
	my $in  = Bio::SeqIO->new(-file => "$multifasta", -format => 'Fasta');
	my $identifier;
	my $sequence;
	
	while ( my $seq = $in->next_seq() ) {
		$identifier = $seq->id;
  		$sequence = $seq->seq;
		
		my ($fasta) = $_;
		$fasta = ">$identifier\n$sequence\n";

		my $count_seq = length($sequence);
		my $count_N = ($sequence =~ tr/N//);
  		my $math = $count_seq * 0.2;

		if ($count_seq > 250){ next; }
		if ($count_N > $math){ next; }

		print $fasta;
	
	}



# END