| 查看: 1710 | 回复: 13 | ||
| 本帖产生 1 个 BioEPI ，点击这里进行查看 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助有关生物信息学的
|
||
| 本人刚接触生物信息学软件，现在老师让我从线粒体基因组全序列中将所有单个的编码基因析出，一个一个的析出太慢了，时间不够，老师自己对这块不懂，让我自己弄，我没学过，希望有人可以指点帮忙 |
回复此楼» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
送红花一朵
12楼2013-04-27 08:21:15
2楼2013-04-22 09:37:32
3楼2013-04-22 10:24:48
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
wizardfan: 金币+10, BioEPI+1, 感谢热心的回答 2013-04-24 06:31:58
感谢参与，应助指数 +1
wizardfan: 金币+10, BioEPI+1, 感谢热心的回答 2013-04-24 06:31:58
|
GO to "http://www.ncbi.nlm.nih.gov/genomes/OrganelleResource.cgi?opt=organelle&taxid=6544" click "Download"-> "download the list of accessions"-> fileA GO to "http://www.ncbi.nlm.nih.gov/sites/batchentrez" upload the fileA and got the all GBk webpage click "send to" to save file with "Genebank full " -> file B install perl and bioperl package run "perl GetCds.pl fileB FileYouWanted" GetCds.pl : #!/usr/bin/perl use Bio::SeqFeatureI; use Bio:  erl;use Bio::Seq; my %species; @seq=read_all_sequences("$ARGV[0]",'genbank'); $flag=1; foreach $seqobj (@seq) { @features=$seqobj->get_SeqFeatures(); foreach $feat (@features) { if($feat->primary_tag eq 'CDS'){ print "$feat\n"; @name=$feat->get_tag_values('gene'); $name=join "",@name; open OUT,">>$ARGV[1]"; print OUT '>'.$seqobj->desc().$name."\n"; print OUT $feat->seq->seq()."\n"; close OUT; } } } |
4楼2013-04-23 07:06:38