1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
#!/usr/local/env perl
# Copyright (c) 2010-201%, Daniel S. Standage and CONTRIBUTORS
#
# The AEGeAn Toolkit is distributed under the ISC License. See
# the 'LICENSE' file in the AEGeAn source code distribution or
# online at https://github.com/standage/AEGeAn/blob/master/LICENSE.
#
# lcl2glb.pl: convert glb2lcl.pl output from locus based coordinates to sequence
# (chromosome or scaffold) based coordinates.
#
# Usage: perl lcl2glb.pl < locus-features-local.gff3 > locus-features-global.gff3
use strict;
while(<STDIN>)
{
if(m/^#/)
{
print unless(m/^##sequence-region/);
next;
}
next if(m/\tlocus\t/);
my @fields = split(/\t/);
my ($ilocuspos) = $fields[8] =~ m/iLocus_pos=([^;]+)/;
my ($ilseqid, $ilstart, $ilend) = $ilocuspos =~ m/(.+)_(\d+)-(\d+)/;
my $offset = $ilstart - 1;
$fields[0] = $ilseqid;
$fields[3] += $offset;
$fields[4] += $offset;
$fields[8] =~ s/iLocus_pos=([^;]+);//;
print join("\t", @fields);
}
|
