linux处理基因组文件（统计行数，求交集，输出行，求和，去重复行）

1. 统计文件行数

wc -l a.bed

2. 求文件区间交集

bedtools intersect -wo -a a.bed -b b.bed > c.bed

3. 输出文件某些行

sed -n '1,2p' c.bed

4. 某一列去重

cat a.bed | awk -F ' ' '{print $2}' | sort | uniq > c.bed

5. 求某一列之和

cat a.bed | awk '{print $4}' | egrep -v "^$"| paste -sd+|bc

6. 筛选某一列的值

cat result_nmi_avg.bed | awk -F ' ' '$4>=10 && $4<=75{print}' > result_nmi_avg2.bed