| 查看: 1887 | 回复: 9 | |||
[交流]
如何用perl读取序列文件后，找出GC最长的一条 已有7人参与
|
|
perl的初学者，觉得很难，需要处理大量的数据，自学中，要写一个程序，我有一个序列文件，计算有多少条序列，平均长度，平均GC值，这些都比较简单。 最难得是需要找到GC含量最高的一条，并返回序列的名称。感觉需要用到hash，用while循环来比较，每比一次用GC值大的那一条序列的名称来赋值，但是从这我开始不会了，求助  最后一个真的不会了，求大神指教  序列大概长这样，很多条    我自己写的前几个问题的程序 |
回复此楼» 猜你喜欢
科研人必读：新一代 KRAS G12C 抑制剂 Olomorasib 的研究进展
已经有0人回复
-
BT-11（Omilancor）研究化合物全解析：从分子结构到实验应用要点
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有185人回复
-
富马酸贝达喹啉：二芳基喹啉类工具化合物的化学特征
已经有0人回复
-
BT投稿with editor状态快两周了，正常吗
已经有25人回复
-
Silevertinib（BDTX-1535）结构拆解：炔基双环片段为何成为合成瓶颈
已经有0人回复
-
RO4929097 解析：双苯并氮杂䓬骨架与 Notch 通路实验设计要点
已经有0人回复
-
分子量 1157 的 PROTAC 怎么配不析出：SGK3-PROTAC1 的溶解、分装与图谱复核
已经有0人回复
-
logP 6.1 的羧酸怎么配母液：GSK 650394 的溶解窗口与稀释顺序
已经有0人回复
-
从"水不溶"到 71 mg/mL：雷特格韦的 pH 依赖溶解与螯合色谱边界
已经有0人回复
-
133 Ų 的极性表面积，如何决定 TNG462 的母液配制与色谱条件
已经有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
