#!usr/bin/perl -w

### Script generates random DNA sequence given 
# GC composition and length of sequence ###

use strict;
use library;

print "\nEnter GC percentage: ";
my $gc = <STDIN>;
chomp $gc;

print "\nEnter length of sequence: ";
my $lenf = <STDIN>;
chomp $lenf;

my $seq = '';

for (my $i = 0; $i < $lenf; ++$i)
	{my $newbase = rndmbp ($gc);
	$seq .= $newbase;
	}

print "\n$seq\n";

print "\nPlease enter filename for outfile.\n\n";
my $outputfile = <STDIN>;
chomp $outputfile;
unless 	(open(OUTFILE, ">$outputfile"))
	{print "Cannot open file \"OUTFILE\" to write to";
	exit;}
print OUTFILE ">Random Sequence ($gc% GC).\n";
print OUTFILE "$seq\n";
close (OUTFILE);


exit;

#####################################