Skip to content

Commit ce2b154

Browse files
committed
Issue #89 | add sw= date= to contig desc
1 parent 58df627 commit ce2b154

File tree

1 file changed

+3
-1
lines changed

1 file changed

+3
-1
lines changed

bin/shovill

Lines changed: 3 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -9,6 +9,7 @@ use File::Temp;
99
use List::Util qw(min max);
1010
use English qw(-no_match_vars);
1111
use Cwd qw(realpath getcwd);
12+
use POSIX qw(strftime);
1213

1314
# Globals
1415
my $EXE = basename($PROGRAM_NAME);
@@ -350,6 +351,7 @@ else {
350351
my $ncontigs = 0;
351352
my $nbases = 0;
352353
my $seq = read_fasta($asm);
354+
my $metadata = "sw=$EXE-$asmdir/$VERSION date=".strftime("%Y%m%d",localtime);
353355
my %len = map { ( $_ => length($seq->{$_}) ) } (keys %{$seq});
354356
for my $id (sort { $len{$b} <=> $len{$a} } keys %{$seq}) {
355357
# spades >NODE_1_length_114969_cov_29.8803_pilon
@@ -376,7 +378,7 @@ for my $id (sort { $len{$b} <=> $len{$a} } keys %{$seq}) {
376378
$cov = sprintf "%.1f", $cov;
377379
my $corr = $changes->{$id} || 0;
378380
# $id =~ s/_pilon$// unless $nocorr;
379-
my $newid = sprintf "$namefmt len=$len cov=$cov corr=$corr origname=$id", $ncontigs;
381+
my $newid = "$namefmt len=$len cov=$cov corr=$corr origname=$id $metadata";
380382
$seq->{$newid} = $seq->{$id};
381383
delete $seq->{$id};
382384
}

0 commit comments

Comments
 (0)