| 查看: 1575 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
张俊亚木虫 (正式写手)
|
[求助]
求人写个程序,生物信息学分析方面的要求用perl或者python来写!已有4人参与
|
||
|
要求其实很简单,具体的理论我就拿个例子来讲吧!是生物信息学方面的! fasta格式的文件 >1 ATCAGATATATATCGGCGG................. >2 .... >45000 ATCAGTATCACAGTG...................,....... 一般能够有几万条序列,从这样的fasta格式文件中,找到以ATCAG(假设,要求这个序列可以自由设定)为开头的序列,并将其变成其反向互补序列(就是把整条序列反过来然后按A-T,C-G转换就行)。最好用python和perl来写,求高手指点!谢谢! |
» 收录本帖的淘帖专辑推荐
交叉知识 |
» 猜你喜欢
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
AI论文写作工具:是科研加速器还是学术作弊器?
已经有3人回复
孩子确诊有中度注意力缺陷
已经有6人回复
2026博士申请-功能高分子,水凝胶方向
已经有6人回复
论文投稿,期刊推荐
已经有4人回复
硕士和导师闹得不愉快
已经有13人回复
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
同一篇文章,用不同账号投稿对编辑决定是否送审有没有影响?
已经有3人回复
ACS Applied Polymer Materials投稿
已经有10人回复
RSC ADV状态问题
已经有4人回复
» 本主题相关价值贴推荐,对您同样有帮助:
yuehedou
木虫 (小有名气)
- 应助: 26 (小学生)
- 金币: 2171.4
- 红花: 2
- 帖子: 215
- 在线: 190.5小时
- 虫号: 1996646
- 注册: 2012-09-12
- 性别: GG
- 专业: 生物信息学

3楼2014-06-08 21:05:25
yuehedou
木虫 (小有名气)
- 应助: 26 (小学生)
- 金币: 2171.4
- 红花: 2
- 帖子: 215
- 在线: 190.5小时
- 虫号: 1996646
- 注册: 2012-09-12
- 性别: GG
- 专业: 生物信息学

2楼2014-06-08 21:02:55
【答案】应助回帖
|
用python试了一下,交流一下 =================================================== #!/usr/bin/python fasta=open('yourfile') fastalist=fasta.readlines() fasta.close() linenum=0 def rev(a): b=a[::-1] c='' for ii in b: if ii=='A': c=c+'T' if ii=='T': c=c+'A' if ii=='C': c=c+'G' if ii=='G': c=c+'C' return c while linenum<len(fastalist): print fastalist[linenum][:-1] if fastalist[linenum+1][:5]=='ATCAG': print rev(fastalist[linenum+1]) else: print fastalist[linenum+1][:-1] linenum+=2 |
4楼2014-06-12 15:23:04
bovine
铜虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1050.1
- 红花: 1
- 帖子: 366
- 在线: 137.5小时
- 虫号: 1825174
- 注册: 2012-05-19
- 性别: GG
- 专业: 微生物遗传育种学
5楼2014-06-13 21:22:26













回复此楼
)))木虫会把右括号转化为笑脸??