利用TBtools工具进行本地blast大致可分为三大步骤：1、得到目的序列所在的scaffold及其起始和终止位置；2、根据左右延长的起始终止位置得到所有外显子的序列；3、将所有外显子序列与query作两两比对，按照query位置将其拼接，从而完成CDS的提取。
对于多外显子基因，此过程较为繁琐，以下命令行则避免了重复性操作，直接将多个物种的各外显子按顺序首尾连接。
（readme：运行目录下包括各个物种的文件，文件内容为该物种该基因的每个外显子序列，其格式为.fasta）
ls > all;for i in `cat all`;do sed -i -e '/>/d' $i;sed -i 'N;s/\n/ /g' $i;sed -i -e 's/\s//g' $i;sed -i -e "1i>$i" $i;done;for j in *.txt;do cat $j >>output.fas;echo >> output.fas;done
然后利用多序列比对，寻找GT-AG的方法找到剪切位点，使其简便而准确的获取目的CDS。
2021-02-28 1.appearance: 'none' 去掉浏览器默认样式
照例先附参考文件： 1. RepeatMasker的安装与使用 2021-02-28 https://www.ji...
周日 2021-02-28 22:00 - 06:30 阴雨 07h31m 记录闪现的灵感(inspiration...
今天的知识点 (2021-02-28) —— 第684天 (我也要出题[http://www.h-camel.co...
【2021-02-28】 持续更新，请关注微信公号 AFoolishPear（微信号：muen_11） 读书笔记0...
【2021-02-28日精进 第284天】表现：9分 今天没有就学习的内容去输出，why？因为华伟需要澄清...
时间：2021-02-28 打卡人：婷婷 打卡天数：89Day 名言警句：再烂的专科也有人挑灯夜读 再好的本科也有...
2021-02-28精进第40天 | 没有记录就没有发生 2021年的事业梦想个人目标描述：达成MDRT业绩目标。...
作者：洞见HeiFa 来源：洞见（ID：2021-02-28） 中国历史上被誉为“圣人”的人物屈指可数，王阳明就是...
2021-02-28，周日，阴雨 人与人的链接、相互吸引其实都是能量吸引。就像磁力，大磁石总是会把小磁石吸到周边。...
本文标题：2021-02-28
本文链接：https://www.haomeiwen.com/subject/gndgfltx.html
网友评论