| 查看: 1864 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
如何用perl读取序列文件后，找出GC最长的一条 已有7人参与
|
|||
|
perl的初学者，觉得很难，需要处理大量的数据，自学中，要写一个程序，我有一个序列文件，计算有多少条序列，平均长度，平均GC值，这些都比较简单。 最难得是需要找到GC含量最高的一条，并返回序列的名称。感觉需要用到hash，用while循环来比较，每比一次用GC值大的那一条序列的名称来赋值，但是从这我开始不会了，求助  最后一个真的不会了，求大神指教  序列大概长这样，很多条    我自己写的前几个问题的程序 |
回复此楼» 猜你喜欢
PRMT5·MTA 抑制剂 MRTX1719 的科研应用
已经有0人回复
-
SW033291 科研使用解析：分子特征、储存与实验获取要点
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有236人回复
-
泛 RAS 抑制新思路：大环肽拟似物如何同时牵住突变型与野生型 RAS？
已经有0人回复
-
TYK2/JAK1 双靶点工具化合物 Brepocitinib：激酶选择性、研究进展与实验应用全解析
已经有0人回复
-
矢车菊素-3-O-葡萄糖苷质控实录：从一张HPLC图谱读懂纯度与稳定性
已经有0人回复
-
恩诺沙星（Enrofloxacin）科研解析：分子特征、LC-MS/MS 检测与实验室使用要点
已经有0人回复
-
Calderasib理化性质与实验操作技术笔记
已经有0人回复
-
RMC-6291怎么用？RAS(ON) 三重复合物工具化合物的实验室实操要点
已经有0人回复
-
苯并呋咱全解析：NBD 荧光衍生化试剂母核的化学特性与实验室应用
已经有0人回复
-
生物质炭与活化剂氢氧化钾混合注意事项
已经有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
