第一步:
就是取开头结尾并存为一个文件
1去重也加进文件, 没重新截图
1查看结果发现,都不是chr,因此转一下:
1脚本:rewrite_chr_in_ncbi_intron_anno.py
需要先做corresponding_relationship_v3.txt对应文件
输出为:rewrite_chr_in_30_min_direction_uniq_v1.bed
1 1第二步:往intron上比,去重,再添加+-
先是比对,脚本和输入输出为:
1之后的都是按脚本,/mnt/T30/wus/brantch_point_human/GRCH37/new/SRR1049830_1.7/bp
结果更少,
网友评论