如何根据同一注释将重叠部分合并

wd_my

A文件4列，序列ID，起始位置，终止位置，注释。目的是将同一注释重叠部分合并得到B文件

Chr01    530   613   Gypsy
Chr01    552   646   Gypsy
Chr01    1605  1765  Gypsy
Chr01    1702  1929   Gypsy
Chr01    2276  2387    BEL
Chr02    2347  2405    BEL
Chr02    4551  4631    BEL
Chr06    4567  4639    DIRS
Chr06    4837  5063    DIRS
Chr06    4862  5038    DIRS
B文件
Chr01    530   646   Gypsy
Chr01    1605  1929  Gypsy
Chr01    2276  2387   BEL
Chr02    2347  2405   BEL
Chr02    4551  4631   BEL
Chr06    4567  4639   DIRS
Chr06    4837  5063   DIRS
大神们，求解答

1cpuer

回复 1# wd_my

awk 'function p(){if(b&&c&&d&&e){print b" "c" "d" "e};b=$1;c=$2;d=$3;e=$4};!a[$1]++{p();next}{if($2>=c&&$2<=d||$3>=c&&$3<=d){c=$2>=c?c2;d=$3<=d?d3}else{p()}}END{p()}' A
Chr01 552 646 Gypsy
Chr01 1605 1929 Gypsy
Chr01 2276 2387 BEL
Chr02 2347 2405 BEL
Chr02 4551 4631 BEL
Chr06 4567 4639 DIRS
Chr06 4837 5063 DIRS

awk 'function p(){if(b""c""d""e){print b" "c" "d" "e};b=$1;c=$2;d=$3;e=$4};!a[$1]++{p();next}{if($2>=c&&$2<=d||$3>=c&&$3<=d){c=$2>=c?c2;d=$3<=d?d3};p()}' A
Chr01 552 646 Gypsy
Chr01 1605 1929 Gypsy
Chr01 1702 1929 Gypsy
Chr01 2276 2387 BEL
Chr02 2347 2405 BEL
Chr02 4551 4631 BEL
Chr06 4567 4639 DIRS
Chr06 4837 5063 DIRS

wh7211

回复 1# wd_my

1. cat A
   
2. Chr01    530   613   Gypsy
   
3. Chr01    552   646   Gypsy
   
4. Chr01    1605  1765  Gypsy
   
5. Chr01    1702  1929   Gypsy
   
6. Chr01    2276  2387    BEL
   
7. Chr02    2347  2405    BEL
   
8. Chr02    4551  4631    BEL
   
9. Chr06    4567  4639    DIRS
   
10. Chr06    4837  5063    DIRS
    
11. Chr06    4862  5038    DIRS
    
12. Chr06    4864  5099    BEL
    
13. 
    
14. awk 'function p(){if(b&&c&&d&&e){print b" "c" "d" "e>"B"};b=$1;c=$2;d=$3;e=$4};!a[$1]++{p();next}{if($2>=c&&$2<=d&&$4==e||$3>=c&&$3<=d&&$4==e){c=$2>=c?c:$2;d=$3<=d?d:$3}else{p()}}END{p()}' A
    
15. 
    
16. cat B
    
17. Chr01 530 646 Gypsy
    
18. Chr01 1605 1929 Gypsy
    
19. Chr01 2276 2387 BEL
    
20. Chr02 2347 2405 BEL
    
21. Chr02 4551 4631 BEL
    
22. Chr06 4567 4639 DIRS
    
23. Chr06 4837 5063 DIRS
    
24. Chr06 4864 5099 BEL

复制代码