第一步:
就是取开头结尾并存为一个文件
![](https://img.haomeiwen.com/i18429961/68aa32131ba2555d.png)
去重也加进文件, 没重新截图
![](https://img.haomeiwen.com/i18429961/69ed979a3d2e8a9e.png)
查看结果发现,都不是chr,因此转一下:
![](https://img.haomeiwen.com/i18429961/29969db425c65168.png)
脚本:rewrite_chr_in_ncbi_intron_anno.py
需要先做corresponding_relationship_v3.txt对应文件
输出为:rewrite_chr_in_30_min_direction_uniq_v1.bed
![](https://img.haomeiwen.com/i18429961/4f80b2fa8702f691.png)
![](https://img.haomeiwen.com/i18429961/18c7688c6c1e1d60.png)
第二步:往intron上比,去重,再添加+-
先是比对,脚本和输入输出为:
![](https://img.haomeiwen.com/i18429961/960fb2540c36629c.png)
之后的都是按脚本,/mnt/T30/wus/brantch_point_human/GRCH37/new/SRR1049830_1.7/bp
结果更少,
网友评论