• 博客园logo
  • 会员
  • 周边
  • 新闻
  • 博问
  • 闪存
  • 众包
  • 赞助商
  • Chat2DB
    • 搜索
      所有博客
    • 搜索
      当前博客
  • 写随笔 我的博客 短消息 简洁模式
    用户头像
    我的博客 我的园子 账号设置 会员中心 简洁模式 ... 退出登录
    注册 登录
 






总是空大

 
 

Powered by 博客园
博客园 | 首页 | 新随笔 | 联系 | 订阅 订阅 | 管理

随笔分类 -  生物信息-进化

我在生物信息之微进化研究中写的一些小程序
 
以12SNP为阈值判断结核菌的传播链
摘要:在结核菌中,以12SNP为阈值判断菌株间是否为近期传播。利用生物信息学方法从原始测序结果得到各菌株的SNP文件,通过互相比较得到两个文件:(1)各菌株间相互比较的情况,格式为:菌株AVS菌株B SNP距离,总行数为n*(n-1)/2 (2)利用文件1得到近期传播的菌株列表接下来要把这些菌株串起来,形 阅读全文
posted @ 2017-05-06 10:06 总是空大 阅读(389) 评论(0) 推荐(0)
判断某突变在NJ法构建的进化树中是否呈现平行进化的趋势
摘要:关于结核分枝杆菌的进化中,最简单的判断菌株受到选择压力的方法就是看是否存在明显的平行进化特征。因此,利用进化树和相关测序数据,可以进行相关判断。import sysimport refrom collections import Iterable#将进化树进行简化tree=""for row in 阅读全文
posted @ 2017-05-06 10:04 总是空大 阅读(564) 评论(0) 推荐(0)