| 查看: 2208 | 回复: 13 | ||
| 本帖产生 1 个 BioEPI ，点击这里进行查看 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助有关生物信息学的
|
||
| 本人刚接触生物信息学软件，现在老师让我从线粒体基因组全序列中将所有单个的编码基因析出，一个一个的析出太慢了，时间不够，老师自己对这块不懂，让我自己弄，我没学过，希望有人可以指点帮忙 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有6人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有6人回复
-
现代”学阀”该如何界定
已经有14人回复
-
求合成方法
已经有8人回复
-
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有4人回复
-
有机合成以后会不会被AI改变？做科研的虫友怎么看
已经有8人回复
-
我的奶奶
已经有3人回复
13楼2013-04-28 13:53:41
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