| 查看: 1843 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
如何用perl读取序列文件后，找出GC最长的一条 已有7人参与
|
|||
|
perl的初学者，觉得很难，需要处理大量的数据，自学中，要写一个程序，我有一个序列文件，计算有多少条序列，平均长度，平均GC值，这些都比较简单。 最难得是需要找到GC含量最高的一条，并返回序列的名称。感觉需要用到hash，用while循环来比较，每比一次用GC值大的那一条序列的名称来赋值，但是从这我开始不会了，求助  最后一个真的不会了，求大神指教  序列大概长这样，很多条    我自己写的前几个问题的程序 |
回复此楼» 猜你喜欢
Wnt-C59：PORCN靶向化学探针，Wnt信号通路阻断机制与科研应用全解析
已经有0人回复
-
Fexaramine：肠道限制性FXR化学探针，结构特征与信号通路研究全解析
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有137人回复
-
NVS-STG2：靶向STING跨膜域的分子胶类激动剂，结构特征与寡聚化机制研究全解析
已经有0人回复
-
DOTA-GPC3肽（CAS 3066439-94-6）：双功能分子探针前体的结构特征与标记应用解析
已经有0人回复
-
科研用 TLR7/8 双激动剂怎么选？现货与质控要点
已经有0人回复
-
SEL120-34A HCl：CDK8/CDK19高选择性激酶抑制剂
已经有0人回复
-
氢氧化钾活化操作
已经有18人回复
-
散金
已经有68人回复
-
生物质超级电容器
已经有0人回复
hxswdl
至尊木虫 (文坛精英)
Tortoise
- 应助: 15 (小学生)
- 贵宾: 2.473
- 金币: 25960.2
- 散金: 52382
- 红花: 399
- 沙发: 46
- 帖子: 13021
- 在线: 1879.7小时
- 虫号: 1764206
- 注册: 2012-04-18
- 专业: 植物学研究的新技术、新方
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
你需要一条一条的序列来处理，共三个思路： 1. 用">"来判断是否为一条序列； 2. 用$/来一条一条的读取序列； 3. 直接使用bioperl来读取每条序列对应的结果。 可以学学华大一位大牛很早写的一个处理fasta的小脚本 https://github.com/JinfengChen/S ... /align/fastaDeal.pl 祝好！ |
9楼2016-11-17 11:33:43
2楼2016-09-06 22:18:30
3楼2016-09-07 10:56:50
4楼2016-09-07 17:22:35
10