使用linux批量引物设计,使用SSRMMD便捷、迅速与准确地进行:SSR位点检测,多态性SSR筛选,与批量SSR引物设计...
1.简介
SSRMMD (Simple Sequence Repeat Molecular Marker Developer) 是用Perl语言编写的软件,可以从组装序列(FASTA格式,例如:基因组或转录组)中检测完美的SSR位点和候选的多态性SSR。该软件还包含了一个名为connectorToPrimer3的程序,此程序提供了SSRMMD到Primer3的接口,使得可以轻松地进行批量SSR引物设计。
2.软件下载与安装
下载地址:https://github.com/GouXiangJian/SSRMMD
安装方法:不管是在linux还是在windows中运行,都是解压即用。
3.软件使用:SSR位点检测
3.1在linux中运行
[重要参数]
解压软件后,将当前工作目录切进软件内部(cd SSRMMD-master),在命令行输入perl SSRMMD.pl -h,即可查看该软件的帮助文档,其中有几个参数对于SSR位点的检测比较重要,例如:
-f1 : 用于检测SSR位点的FASTA格式的文件 (必须提供该参数!)
-e : 指定检测SSR位点的方法 (默认设置: 0, 可选的值: 0 [这种方法速度更快], 1 [这种方法类似于MISA软件使用的方法])
-mo : SSR基序的阈值 (默认设置: 1=10,2=7,3=6,4=5,5=4,6=4 [其中, 等号左边是基序的长度, 等号右边是最小重复次数])
-l : SSR的侧翼序列的长度 (默认设置: 100)
-ss : 是否输出SSR的统计文件 (默认设置: 0, 可选的值: 0 [表示不输出], 1 [表示输出])
-t : 运行软件时使用的线程数 (默认设置: 1)
[运行示例]
以软件自带的例子文件example1.fa举例,如果想要:(1) 修改SSR基序的阈值设置,(2) 将SSR的侧翼序列设置为200bp,(3) 输出SSR的统计文件,(4) 使用2线程来运行程序,则命令行可以这样输入:
perl SSRMMD.pl -f1 example/input/example1.fa -mo 2=6,3=5,4=4 -l 200 -ss 1 -t 2
[运行输出]
运行完成后,会在当前目录下创建一个名为SSRMMDOUT的目录,里面包含SSR信息文件example1.fa.SSRs,以及SSR统计文件example1.fa.stat。
3.2在windows中运行
如果想在windows中用上述同样的命令运行SSRMMD,则需要事先安装Perl的解释器,下载地址为:https://www.perl.org/。好消息是,该软件包中提供了一个已经编译好的版本SSRMMD.exe,位于bin目录中。因此,可以直接使用它而不必再安装Perl的解释器了。
解压软件后,将当前工作目录切进软件内部(cd SSRMMD-master/bin),在命令行输入SSRMMD.exe -h,即可查看该软件的帮助文档。
[运行示例]
若用上述同样的参数配置来运行,则命令行可以这样输入:
SSRMMD.exe -f1 ../example/input/example1.fa -mo 2=6,3=5,4=4 -l 200 -ss 1 -t 2
[运行输出]
类似地,运行完成后,会在当前目录下创建一个名为SSRMMDOUT的目录,里面包含SSR信息文件example1.fa.SSRs,以及SSR统计文件example1.fa.stat。
4.软件使用:多态性SSR筛选
4.1在linux中运行
[重要参数]
如果想要进一步筛选多态性SSR,则必须准备两个fasta格式的组装文件。SSRMMD会在这两个文件之间,对SSR的侧翼序列率先进行保守性评估,随后进行唯一性评估。
类似地,帮助文档中有几个参数对于多态性SSR筛选比较重要,例如:
-me : 指定检测侧翼序列保守性的算法 (默认设置: NO, 可选的值: NO [仅仅使用HASH, 该方法非常快, 但是只能检索出侧翼序列绝对保守的SSR], LD [使用Levenshtein-Distance算法], NW [使用Needleman–Wunsch算法])
-d : 当参数-me设置为LD时, 设置由Levenshtein-Distance算法计算的侧翼序列保守性的阈值 (默认设置: 0.05)
-i : 当参数-me设置为NW时, 设置由Needleman-Wunsch算法计算的侧翼序列保守性的阈值 (默认设置: 0.95)
-st : 指定检测侧翼序列唯一性的计算方式 (默认设置: 0, 可选的值: 0 [节约时间], 1 [节约内存])
[运行示例]
以软件自带的例子文件example1.fa和example2.fa举例,如果想要:(1) 使用Needleman–Wunsch算法来评估侧翼序列的保守性,(2) 以节约时间的方式来评估侧翼序列的唯一性,(3) 使用2线程来运行程序,则命令行可以这样输入:
perl SSRMMD.pl -f1 example/input/example1.fa -f2 example/input/example2.fa -p 1 -me NW -st 0 -t 2
[运行输出]
运行完成后,会在当前目录下创建一个名为SSRMMDOUT的目录,里面包含SSR信息文件example1.fa.SSRs和example2.fa.SSRs,以及多态SSR记录文件example1.fa-and-example2.fa.compare。
4.2在windows中运行
[运行示例]
若用上述同样的参数配置来运行,则命令行可以这样输入:
SSRMMD.exe -f1 ../example/input/example1.fa -f2 ../example/input/example2.fa -p 1 -me NW -st 0 -t 2
[运行输出]
类似地,运行完成后,会在当前目录下创建一个名为SSRMMDOUT的目录,里面包含SSR信息文件example1.fa.SSRs和example2.fa.SSRs,以及多态SSR记录文件example1.fa-and-example2.fa.compare。
5.软件使用:批量SSR引物设计
5.1在linux中运行
目录connectorToPrimer3中的脚本connectorToPrimer3.pl可以轻松地将SSRMMD与Primer3结合起来,以实现批量SSR引物设计。在命令行输入perl connectorToPrimer3/connectorToPrimer3.pl -h,即可查看该脚本的帮助文档。
若用默认的参数配置设计引物,则命令行可以这样输入:
#Design primers for all SSRs (SSRMMD option '-p' = 0)
perl connectorToPrimer3/connectorToPrimer3.pl -i SSRMMDOUT/example1.fa.SSRs -o all.txt
#Design primers for candidate polymorphic SSRs (SSRMMD option '-p' = 1)
perl connectorToPrimer3/connectorToPrimer3.pl -i SSRMMDOUT/example1.fa-and-example2.fa.compare -o poly.txt -s 2
5.2在windows中运行
类似地,脚本connectorToPrimer3.pl也有一个已经编译好的版本connectorToPrimer3.exe,位于bin目录中。
若用默认的参数配置设计引物,则命令行可以这样输入:
#Design primers for all SSRs (SSRMMD option '-p' = 0)
connectorToPrimer3.exe -i SSRMMDOUT/example1.fa.SSRs -o all.txt
#Design primers for candidate polymorphic SSRs (SSRMMD option '-p' = 1)
connectorToPrimer3.exe -i SSRMMDOUT/example1.fa-and-example2.fa.compare -o poly.txt -s 2
使用linux批量引物设计,使用SSRMMD便捷、迅速与准确地进行:SSR位点检测,多态性SSR筛选,与批量SSR引物设计...相关推荐
- UI设计灵感|手机钱包App怎样设计才更便捷?
手机支付应该给用户带来的是方便和快捷,支付和财务管理等常用功能也该用最少的步骤来实现. 相信大多数人如今出门已经不带现金了,无现金生活为我们带来了巨大的便利,但众多的资产管理和理财类应用仍然在进行不断 ...
- 《安富莱嵌入式周报》第283期:全开源逆向“爆破”硬件工具,Linux内核6.1将正式引入RUST语言,I3C培训教程,80款市场成熟的电感式位置传感器设计
往期周报汇总地址:嵌入式周报 - uCOS & uCGUI & emWin & embOS & TouchGFX & ThreadX - 硬汉嵌入式论坛 - P ...
- 【检测技术课案】直流低电阻测试仪的设计与制作(一)
1 引言 1.1 研究背景与意义 在生产和实践中,常常需要对低电阻进行精确测量,如电机和变压器的线圈电阻,电闸接线端子和电缆插座的接触电阻,大功率电器开关的接触电阻等[1],这些电阻一般在MΩ级甚至μ ...
- 学生网页作业:商城网站设计——仿天猫商城(9页) HTML+CSS大作业_ 网页制作作业_疫情防控网页设计...
HTML5期末大作业:商城网站设计--仿天猫商城(9页) HTML+CSS大作业 文章目录 HTML5期末大作业:商城网站设计--仿天猫商城(9页) HTML+CSS大作业 一.作品展示 二.文件目录 ...
- HTML5期末大作业:动漫网站设计——迪斯尼公主(6个页面) HTML+CSS+JavaScript 动漫网页HTML代码 学生网页课程设计期末作业下载 动漫大学生网页设计制作成
HTML5期末大作业:动漫网站设计--迪斯尼公主(6个页面) HTML+CSS+JavaScript 动漫网页HTML代码 学生网页课程设计期末作业下载 动漫大学生网页设计制作成 临近期末, 你还在为 ...
- 无线电能传输 wpt 磁耦合谐振 过零检测 matlab simulink仿真 pwm MOSFET,过零检测模块 基于二极管整流的无线电能传输设计
无线电能传输 wpt 磁耦合谐振 过零检测 matlab simulink仿真 pwm MOSFET,过零检测模块 基于二极管整流的无线电能传输设计 基于同步整流的无线电能传输设计(含过零比较) 两个 ...
- HTML5期末大作业:动漫网站设计——迪斯尼公主(6个页面) HTML+CSS+JavaScript 动漫网页HTML代码 学生网页课程设计期末作业下载 动漫大学生网页设计制作成...
HTML5期末大作业:动漫网站设计--迪斯尼公主(6个页面) HTML+CSS+JavaScript 动漫网页HTML代码 学生网页课程设计期末作业下载 动漫大学生网页设计制作成 临近期末, 你还在为 ...
- 批量检测商品是否下架,批量查询商品信息接口
批量检测商品是否下架,批量查询商品信息接口代码展示如下: ###[请求参数,支持高并发] 参数: ids 参数:必选,你要检测的批量商品ID列表,用英文逗号相连.如 646724659836,6503 ...
- 软件设计是怎样炼成的(8)——用户感觉好才是真的好(用户体验设计)
摘要: 软件的UI是客户直接体验软件的地方,好的用户体验可能比追求新技术和追求漂亮设计更加重要,优秀的软件应该是既有外在美也有内在美的.但大部分用户是"外貌协会"的,你和用户说你的 ...
- HTML5期末大作业:爱折扣电商购物网站设计——清新爱折扣电商购物网站(5页) HTML+CSS+JavaScript 关于电商的HTML网页设计-----网购
HTML5期末大作业:爱折扣电商购物网站设计--清新爱折扣电商购物网站(5页) HTML+CSS+JavaScript 关于电商的HTML网页设计-----网购 常见网页设计作业题材有 个人. 美食. ...
最新文章
- 厉害!从电影花瓶到 Wi-Fi 之母,这才是乘风破浪的姐姐!
- Openstack组件部署 — 将一个自定义 Service 添加到 Keystone
- 信息学奥赛一本通(1218:取石子游戏)
- Facebook何恺明团队提出SlowFast网络,视频识别无需预训练
- 832计算机专业基础,2019年考试科目832计算机组成原理与数据结构.doc
- Q124:PBRT-V3,“路径追踪”积分器(14.5章节)
- Docker入门者手册
- Tomcat—logs文件夹中不再产生日志文件
- oppoa1计算机记录删了怎么办,捡到oppoA1怎么解锁
- 2×125MW发电厂升压站电气部分设计
- android+桌面+横屏,安卓怎么强制桌面横屏
- elasticjob-配置手册
- ctf-wp-mmmmpy
- 怎么打开设备管理器(有图 超详细)
- WEB服务器搭建与配置
- 从一个NPU失效问题看Linux PM Domain Framework的实现逻辑
- python 火车票订票系统(12306)
- 国内众多云桌面你更中意谁
- 关于EF多线程更新数据的一个报错
- python制作会动的表情包_Python自动生成表情包