From 27e2d8a6fa339141effde502af17225a6b716b31 Mon Sep 17 00:00:00 2001
From: daisieh_local <daisieh@gmail.com>
Date: Tue, 8 Apr 2014 14:24:12 -0700
Subject: [PATCH] Fixes Issue #132
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Makes a temp file that has the reference sequence named “reference” so
that percentcoverage doesn’t choke.
---
 Postprocessing/PercentCoverage.pl | 12 +++++++++++-
 1 file changed, 11 insertions(+), 1 deletion(-)

diff --git a/Postprocessing/PercentCoverage.pl b/Postprocessing/PercentCoverage.pl
index e001e9b..a450a80 100755
--- a/Postprocessing/PercentCoverage.pl
+++ b/Postprocessing/PercentCoverage.pl
@@ -5,6 +5,8 @@
 use File::Basename qw(basename);
 use lib "$FindBin::Bin/../lib";
 use Postprocessing;
+use Parsing qw(parsefasta);
+use File::Temp qw(tempfile);
 
 my $ref_file = shift @ARGV;
 my $contigs_file = shift @ARGV;
@@ -24,7 +26,15 @@
 	die "Couldn't find contigs file $contigs_file";
 }
 
-my $contigs = percentcoverage ($ref_file, $contigs_file, $out_name, $aligner);
+my ($reffasta, $reffastaarray) = parsefasta ($ref_file);
+my ($fh, $filename) = tempfile();
+print $fh ">reference\n";
+foreach my $r (@$reffastaarray) {
+	print $fh "$reffasta->{$r}\n";
+}
+close $fh;
+
+my $contigs = percentcoverage ($filename, $contigs_file, $out_name, $aligner);
 
 if (defined $contigs) {
 	my $refseq = delete $contigs->{reference};