常用命令

1.NR注释作图 /home/xwm/basic/RNA_seq.NR.sh

2.上下调gene获取命令：

awk '$10~/down/ && $12~/yes/{print $1}' genes.diff_exp_result > 526_vs_332_down.list

awk '$10~/up/ && $12~/yes/{print $1}' genes.diff_exp_result > 526_vs_332_up.list

3.NR注释统计作图/home/xwm/basic/RNA_anno_NR.KEGG.GO_level.pdf.sh

4.运行trinity butterfly部分： /mnt/lustre/users/liubinxu/workspace/run_trinity_pbs.pl

5.KEGG pathway 分级统计 /mnt/lustre/users/liubinxu/workspace/pathway_sort.pl -pathway pathway_table.xls

6.pathway top 20 perl /home/wangyan/script/kegg_pathway_top20bars.pl -t pathway_table.xls

7.KEGG_add_exp2.pl的补充程序for i in $(find ./ -name "*.m"); do /mnt/lustre/users/renyi/bio-softs/matlab2013/bin/matlab -nojvm -nodisplay < $i; done

8.互换两列内容：cat genes.diff_exp_result  |awk -F "\t" 'BEGIN{OFS="\t"};{a=$8;$8=$10;$10=a; print $0}' > genes.diff_exp_result1