m0_51268845的博客

ll /media/xiaonan/14T4/xiaohong/stats_TE/TE_stat/*sum.stat|awk -F "Zm-" '{print $2}'|awk -F "TE" '{print $1}'|tr -d '.'|sed 's/.*/"&"/' | paste -sd, -|cat

prev="";next;

ll /home/zhangsan/anaconda3/envs/*/bin/mafft

Alt + Tab Tab可以选择软件。

grep -v '^$' jop > jop1

retab。

grep -f <(less /media/spl/14T4/huj/natural_mutant_analysis/top_100.txt|awk -F "," '{print $1","}'|sed '1d'|less) /media/spl/14T4/huj/natural_mutant_analysis/hhj_gene_using_arabidopsis.annot.csv

sed 's/^[ \t]*//' |sed '1i\gene'

/home/zhangshan/anaconda3/envs/*/lib/*/site-packages/requests

less hwc.snp_mis9_maf5.vcf|grep -v "##"|grep -E '#|28824938' |cat

grep 28824938 abc.snp_mis9_maf5.vcf|tr "\t" "
"|awk -F ":" '{print $1}'|sed -n '10,$p'|sort|uniq -c

在当前目录会出现.Rhistory文件，这里面的代码会包含很多之前python的脚本，需要选取最后n行。

if info.get(key) is None:

','')','')

check.names=FALSE

grep NM_179830.1 /media/songxh/14T4/annotation_blast/data/arabidopsis/blast/target_result|cut -f 2|sed 's#model#TU#g'|awk '{print "\""$1"\""}'|tr '
' ','|cat

f = open(sys.argv[1],'r')f1 = open(sys.argv[2],'w')f2 = open(sys.argv[3],'r')

储存所有exon位置信息。# 储存所有exon位置信息。# 计算多个外显子长度。# 计算多个外显子长度。# 信息保存在字典里。# 打开测试 gtf 文件。# 取并集(删除重复元素)# 读取注释类型文件。

if($1!

sort -k 1,1 -k 5,5g target_result|less -SN