| 查看: 1836 | 回复: 13 | ||
| 本帖产生 1 个 BioEPI ，点击这里进行查看 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助有关生物信息学的
|
||
| 本人刚接触生物信息学软件，现在老师让我从线粒体基因组全序列中将所有单个的编码基因析出，一个一个的析出太慢了，时间不够，老师自己对这块不懂，让我自己弄，我没学过，希望有人可以指点帮忙 |
回复此楼» 猜你喜欢
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有6人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有6人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有9人回复
-
有院领导为了换新车，用横向课题经费买了俩车
已经有10人回复
-
遇见不省心的家人很难过
已经有24人回复
-
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有6人回复
-
酰胺脱乙酰基
已经有13人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
erlei11楼2013-04-27 07:46:53
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