Bash useful operation
awkawk '$4>1000' a.txt ##过滤第四列大于1000的,输出来这些行awk '$1>1000 && $2=="mRNA" {print $1,$2,$3}' a.txt #过滤第一列大于1000,第二列等于mRNA,输出只带第一列,第二列,第三列的所有行awk '$2~/^scaffold/' a.txt #输出第二列开头匹配scaffold的行
grep!warning: grep -c can only count the number of times a string
...