| 查看: 1593 | 回复: 9 | |||
[交流]
如何用perl读取序列文件后，找出GC最长的一条 已有7人参与
|
|
perl的初学者，觉得很难，需要处理大量的数据，自学中，要写一个程序，我有一个序列文件，计算有多少条序列，平均长度，平均GC值，这些都比较简单。 最难得是需要找到GC含量最高的一条，并返回序列的名称。感觉需要用到hash，用while循环来比较，每比一次用GC值大的那一条序列的名称来赋值，但是从这我开始不会了，求助  最后一个真的不会了，求大神指教  序列大概长这样，很多条    我自己写的前几个问题的程序 |
回复此楼» 猜你喜欢
筑牢营养安全线：以精准检测，护健康基石
已经有0人回复
-
推荐一些20种氨基酸检测的实际应用案例
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有150人回复
-
不合理蛙科研实验中的趣事：实验器材的 “乌龙”
已经有0人回复
-
不合理蛙科研实验中的趣事：和实验材料的 “斗智斗勇”
已经有0人回复
-
蛋白质检测：精准分析，解锁生物分子的密码
已经有0人回复
-
不合理蛙科研实验之小鼠实验：严谨设计，解析生命机制的重要载体
已经有0人回复
-
不合理蛙科研实验之重金属检测：精准筛查，守护健康与环境的防线
已经有0人回复
-
不合理蛙科研实验之“蛙测重金属我背锅三千”
已经有0人回复
-
不合理蛙科研实验之“鼠逃三次我发三篇SCI”
已经有0人回复
2楼2016-09-06 22:18:30
3楼2016-09-07 10:56:50
4楼2016-09-07 17:22:35
ygidtu
新虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1690.7
- 散金: 798
- 红花: 3
- 帖子: 419
- 在线: 117.7小时
- 虫号: 2854989
- 注册: 2013-12-07
- 性别: GG
- 专业: 生物信息学
5楼2016-09-10 12:32:21
6楼2016-09-10 19:18:05
7楼2016-09-10 20:51:33
yong5060269
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 308.9
- 散金: 12
- 红花: 1
- 帖子: 50
- 在线: 19.1小时
- 虫号: 3590596
- 注册: 2014-12-12
- 性别: GG
- 专业: 生物大分子结构与功能
8楼2016-09-14 19:06:12
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
10楼2016-12-09 22:03:09