| 1 |
gpertea |
23 |
#!/usr/bin/perl |
| 2 |
|
|
|
| 3 |
|
|
use strict; |
| 4 |
|
|
my ($file,$addcoord)=@ARGV; |
| 5 |
|
|
|
| 6 |
|
|
open(F,$file) || die ("Error opening file $file\n"); |
| 7 |
|
|
my $p=''; |
| 8 |
|
|
if ($file=~m/\.(\d+)$/) { |
| 9 |
|
|
$p=$1; |
| 10 |
|
|
# -- only needed for correct naming of gff3 entitied |
| 11 |
|
|
# produced by gmap ran on separate parts |
| 12 |
|
|
|
| 13 |
|
|
} |
| 14 |
|
|
while(<F>) { |
| 15 |
|
|
chomp; |
| 16 |
|
|
if (substr($_,0,1) eq "#") { print $_."\n";} |
| 17 |
|
|
else { |
| 18 |
|
|
if ($p) { |
| 19 |
|
|
s/\.cds(\d+)/\.cds$p\_$1/g; |
| 20 |
|
|
s/\.path(\d+)/\.p$p\_$1/g; |
| 21 |
|
|
} |
| 22 |
|
|
my @a=split; |
| 23 |
|
|
$a[3]+=$addcoord; |
| 24 |
|
|
$a[4]+=$addcoord; |
| 25 |
|
|
print join("\t", @a)."\n"; |
| 26 |
|
|
} |
| 27 |
|
|
} |
| 28 |
|
|
close(F); |
