- 论坛徽章:
- 307
|
回复 22# jiangphoenix
- #!/usr/bin/perl
- use strict;
- use warnings;
- my @aData = [split (' ', scalar <DATA>)];
- while (<DATA>){
- my ($id, $start, $end) = split;
- if ($aData[-1][0] eq $id){
- $start - $aData[-1][-1] < 3 ? ($aData[-1][-1] = $end) : push (@aData, [$id, $start, $end]);
- next;
- }
- push (@aData, [$id, $start, $end]);
- }
- print "@$_\n" for @aData;
- __DATA__
- chr1 11129265 11129312
- chr1 11132717 11132739
- chr1 11132741 11132745
- chr1 11132747 11132764
- chr1 11132766 11132818
- chr1 11132850 11132904
- chr1 11132906 11132939
- chr1 11139743 11139968
- chr1 11150021 11150225
- chr1 11172370 11172752
- chr1 11174848 11175076
- chr1 11175999 11176441
- chr1 11177029 11177309
- chr1 11177311 11177397
- chr1 11177453 11177683
- chr1 11181703 11182355
- chr1 11187182 11187371
- chr1 11193359 11193591
- chr1 11195078 11195149
复制代码
|
|