博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
php 统计fasta 序列长度和GC含量
阅读量:6843 次
发布时间:2019-06-26

本文共 740 字,大约阅读时间需要 2 分钟。

最近php7的消息铺天盖地, 忍不住想尝试下。星期天看了下语法, 写个小脚本练下手;

这个脚本读取fasta 文件, 输出序列的长度和GC含量;

/", $line)) { if ($id) { $meta[$id] = $seq; $id = $line; $seq = ''; } else { $id = $line; } } else { $seq .= $line; } } $meta[$id] = $seq; fclose($file_handle); return $meta;}function write_res($meta) { foreach ($meta as $key => $value) { $len = cal_length($value); $gc = cal_gc($value); echo "$key\t$value\t$len\t$gc\n"; }}function cal_length($seq) { return strlen($seq);}function cal_gc($seq) { $gc = array(); preg_match_all("/G|C|g|c/", $seq, $gc); return count($gc[0]) / strlen($seq);}?>

 

转载地址:http://smbul.baihongyu.com/

你可能感兴趣的文章
iOS概念攻坚之路(二):Runtime
查看>>
关于前端请求发送时间时而长时而短问题(stalled a lot)
查看>>
Python 工匠:编写条件分支代码的技巧
查看>>
记一次前端面试经历
查看>>
带你探索JUnit 5.4
查看>>
<暗时间> 时间, 不在于你拥有多少, 而在于你怎样使用
查看>>
单例 - iOS
查看>>
戛纳电影节百花齐放,中国明星衣着品味紧跟时尚前沿
查看>>
mysql 存储emoji表情
查看>>
10年测试总监经验分享,你与优秀工程师的距离!
查看>>
2019年在哪里找好的高层次人才扶持政策?
查看>>
解决代码报红:Cannot resolve symbol 'xxx'
查看>>
第71节:Java中HTTP和Servlet
查看>>
Linux开源CommunityBridge平台 提供资金、安全以及人员三项关键
查看>>
Python爬虫入门教程 5-100 27270图片爬取
查看>>
Day1:html和css
查看>>
开源如何在云上存活?
查看>>
Android 网络基础之 HTTP
查看>>
ES6实现继承
查看>>
有擎企业系统v1.0.0 积木式搭建网站,页面构建更灵活
查看>>