Advertisement

计算反转录转座子插入时间三 MEGA批量化处理

阅读量:

在获取成对的LTR序列之后,可以根据两条LTR序列之间的差异性,依据K-2模型计算每个核苷酸位点的平均替换次数,该计算过程可通过MEGA软件完成。由于涉及的文件数量较多,需要进行多次操作,因此采用MEGA-cc工具进行批量处理。具体的安装方式可参考

完成安装后,启动MEGA程序,点击界面右下角的PROTOTYPE选项,进入模拟模式,选择序列类型(此处选用Nucleotide non-coding),在ALIGN选项中选择MUSCLE算法,并采用默认参数确认后生成比对配置文件*.mao。

在命令行环境下调用MEGA-CC工具:
M6CC.exe -a *.mao -d inputfile -o outputfile
鉴于文件数量庞大,将所有文件路径及名称记录于mytest.txt文档中,并再次调用MEGA-CC工具:
M6CC.exe -a *.mao -d mytest.txt -o outputfile
处理完成后,在MEGA中切换至模拟模式,选择DISTANCE作为序列类型,并根据以下参数生成用于计算遗传距离的配置文件:

在这里插入图片描述

基于前一步比对所得结果作为

全部评论 (0)

还没有任何评论哟~