带行头与删除重复

带行头一句话搞定！！！！awk -F"\t" ‘{if(/^#/){print $0}else{if(($3-$2)==10000){print $0}}}’ merged_loops.bedpe > merged_loops_10k.bedpe
合并
tail postprocessed_pixels_3000.bedpe
tail postprocessed_pixels_5000.bedpe
cat postprocessed_pixels_3000.bedpe postprocessed_pixels_5000.bedpe >prost35 ###这种是没有去头的！！

1、删除行头空格键和tab键
cat seddlt | sed ‘s/¹*//g’
sed删除空行,#开头的行,以及第一行
删除空行
sed “/^$/d”
删除#开头的行
sed “/^#/d”
删除第一行:
sed -i “1d” a.txt

模式用//括起，扩展正则要加 -r选项。

sed -n '/模式/p’ 输入文件
显示包含数字2的行：
[root@liupeng lp]# sed -n ‘/2/p’ sed.txt
2
12
sed -n ‘/#/p’ prost35
(base) [zhangyong@cluster dm]$ sed -n ‘/#/p’ prost35
##chr1 x1 x2 chr2 y1 y2 name score strand1 strand2 color observed expectedBL expectedDonut expectedH expectedV fdrBL fdrDonut fdrH fdrV numCollapsed centroid1 centroid2 radius

juicer_tools version 1.9.9

##chr1 x1 x2 chr2 y1 y2 name score strand1 strand2 color observed expectedBL expectedDonut expectedH expectedV fdrBL fdrDonut fdrH fdrV numCollapsed centroid1 centroid2 radius

juicer_tools version 1.9.9

删除空行和注释：

[root@liupeng lp]# sed -r '/^KaTeX parse error: Expected group after '^' at position 2: |^̲#/d' sed.txt -…指空行；^# 以#开头的行，即注释行

sed ‘5,$s/aa/AA/g’ test.txt
表示对第5行直到文件末尾的所有行进行搜索替换操作，同样s命令的地址限定也支持使用正则表达式限定符合条件的行，然后在这些行中进行字符串的搜索替换操作，例如：
sed ‘/²/s/aa/AA/g’ test.txt

打印从第5行开始到第一个包含以test开始的行之间的所有行：
sed -n ‘5,/^test/p’ file

(base) [zhangyong@cluster dm]$ grep -n “#” prost35
1:##chr1 x1 x2 chr2 y1 y2 name score strand1 strand2 color observed expectedBL expectedDonut expectedH expectedV fdrBL fdrDonut fdrH fdrV numCollapsed centroid1 centroid2 radius
2:# juicer_tools version 1.9.9
11799:##chr1 x1 x2 chr2 y1 y2 name score strand1 strand2 color observed expectedBL expectedDonut expectedH expectedV fdrBL fdrDonut fdrH fdrV numCollapsed centroid1 centroid2 radius
11800:# juicer_tools version 1.9.9
(base) [zhangyong@cluster dm]$

[zhangyong@cluster dm]$ wc -l prost35
26835 prost35
sed -i “11799,11800d” prost35 ###好像直接操作了源文件，重定向没用！！！确实是这样的！！！把那两行删除了
(base) [zhangyong@cluster dm]$ wc -l prost35
26833 prost35
(base) [zhangyong@cluster dm]$ grep -n “#” prost35
1:##chr1 x1 x2 chr2 y1 y2 name score strand1 strand2 color observed expectedBL expectedDonut expectedH expectedV fdrBL fdrDonut fdrH fdrV numCollapsed centroid1 centroid2 radius
2:# juicer_tools version 1.9.9
(base) [zhangyong@cluster dm]$

---

1. \t ↩︎

2. 0-9 ↩︎