trim_galore双端数据过滤nohup trim_galore -q 20 --phred33 --stringency 3 --leng...

将基因fa或蛋白fa序列文件由多行转换为一行awk 'NR==1{a[NR]=$0}NR>1{if(a[NR-1] ~ ">"){print a...

查看自己的文件 cat -A mmgs.txt 去除M,注意：M的输入方式是 Ctrl + v ，然后Ctrl + M sed -i 's/^M...

查找TPS相关基因 grep -wi --color terpene.*synthase combined_function_annotatio...

自己测序的基因组，需要查找对应基因号 用awk命令查找对应基因号及下一行 awk '$0~/PTG008298L/{print $0; getl...

香榧基因组index构建 激活conda环境 conda activate rna 构建hisat2索引后台运行 nohup hisat2-bu...