https://github.com/gringer/bioinfscripts
Tip revision: cc40b6b0edd27642d3d40d46ade1b9dbf9686c20 authored by David Eccles (gringer) on 07 September 2019, 11:17:37 UTC
Updated README.md with GitLab message
Updated README.md with GitLab message
Tip revision: cc40b6b
blastx_boundaries.pl
#!/usr/bin/perl
use warnings;
use strict;
my %starts = ();
my %ends = ();
while(<>){
chomp;
my @F = split(/\s+/);
my $qs=$F[6];
my $qe=$F[7];
if($qe < $qs){
($qs, $qe) = ($qe, $qs);
}
my $sig = $F[0]."-".$F[1];
if(!$starts{$sig} || ($starts{$sig} > $qs)){
$starts{$sig} = $qs
}
if(!$ends{$sig} || ($ends{$sig} < $qe)){
$ends{$sig} = $qe
}
}
foreach my $key (keys(%starts)){
if($key =~ /(.*)-(.*)/){
my $target = $1;
my $query = $2;
printf("%s %s:%d-%d\n", $query, $target, $starts{$key}, $ends{$key});
}
}