dna组装验证的引物选择+数据分析

primavera的Python项目详细描述


primavera是一个python库,用于计划和分析基于引物的dna组装验证,使用sanger测序或验证pcr。它实现了设计和选择引物的方法,以确保相关的装配片段将被覆盖,并可以产生简单的(但近似的)情节总结了一批桑格测序的结果。

底漆选择示例

下面的代码假设文件available_primers.fa包含实验室中所有可用引物的标签和序列,并且要序列验证的程序集具有注释,指示序列应覆盖的区域和应避免底漆退火的区域。

https://raw.githubusercontent.com/Edinburgh-Genome-Foundry/Primavera/master/docs/_static/images/annotated_genbank.png
fromprimaveraimportPrimerSelector,Primer,load_recordimportos# LOAD ASSEMBLIES RECORDS AND AVAILABLE PRIMERSrecords=[load_record(file_path,linear=False)forfile_pathin['my_record_1.gb','my_record_2.gb'...]]available_primers=Primer.list_from_fasta("example_primers.fa")# SELECT THE BEST PRIMERSselector=PrimerSelector(tm_range=(55,70),size_range=(16,25))selected_primers=selector.select_primers(records,available_primers)# PLOT THE COVERAGE AND WRITE THE PRIMERS IN A SPREADSHEETselector.plot_coverage(records,selected_primers,'coverage.pdf')selector.write_primers_table(selected_primers,'selected_primers.csv')

返回的selected_primers包含一个引物列表(每个构造一个列表)。返回的pdf报告如下:

https://raw.githubusercontent.com/Edinburgh-Genome-Foundry/Primavera/master/docs/_static/images/annotated_primer_selection.png

信息

pip安装:

pip install primavera

web文档:

https://edinburgh-genome-foundry.github.io/Primavera/

github页

https://github.com/Edinburgh-Genome-Foundry/Primavera

现场演示

http://cuba.genomefoundry.org/select_primers

许可证:麻省理工学院,爱丁堡基因组铸造厂版权所有

更多生物软件

https://raw.githubusercontent.com/Edinburgh-Genome-Foundry/Edinburgh-Genome-Foundry.github.io/master/static/imgs/logos/egf-codon-horizontal.png

primavera是用于dna设计、制造和验证的合成生物学软件套件的一部分。

欢迎加入QQ群-->: 979659372 Python中文网_新手群

推荐PyPI第三方库


热门话题
java找不到当前线程的会话   同步Java同步方法调用非同步方法   这段代码的java时间复杂度是如何列出一个集合的所有子集的?   java Jetty SslConnectionFactory nextProtocol   为什么在java集合中使用游标而不是for循环   java如何对长函数进行单元测试?   如何从java中的文本输入创建数组   java工具栏内容未显示,尽管出现了片段内容   java MediaPlayer无法处理本地文件   java在多次读取时,从TCP服务器读取的数据将挂起   java组织。springframework。数据领域无法将PageImpl强制转换为   java如何使用@WithMockUser并从属性文件中传递用户名和密码?   JavaSpringBootGradle插件2。x如何使用自定义启动器布局工厂   sockets Java客户端和服务器之间使用字符串等进行通信的“行业标准”方法是什么   java变量引用不存在的资源:${project\u loc:project\u name}   搜索如何返回LDAP Java中所有用户的特定属性?   JDK1.8.016x的java密钥斗篷代理问题   spring如何使用Java中的Exchange获取邮件   java Kotlin+Dagger 2:Dagger*文件不会生成   如何在传入的double上执行操作,然后在Java中返回