24小时热门版块排行榜    

查看: 1862  |  回复: 13
本帖产生 1 个 BioEPI ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

erlei

木虫 (正式写手)

[求助] 求助有关生物信息学的

本人刚接触生物信息学软件,现在老师让我从线粒体基因组全序列中将所有单个的编码基因析出,一个一个的析出太慢了,时间不够,老师自己对这块不懂,让我自己弄,我没学过,希望有人可以指点帮忙
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

genomelin

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
erlei: 金币+5, ★★★很有帮助 2013-04-24 07:59:48
引用回帖:
8楼: Originally posted by erlei at 2013-04-23 05:06:25
不好意思,前面电脑卡了,没传上去,这次传上去了...

将 tag “gene” 改为 “product”
因为你的文件中有些不是统一的mito GBK格式
推荐只用 NC 类的,其它的有可能有错误
另外NC 类也是有注释错误的,需要你做个简单clustaw,看是不是有的注释到了相反链上

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 123.fa
  • 2013-04-23 17:09:25, 1.69 M
9楼2013-04-23 17:12:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

137167741

至尊木虫 (著名写手)

★ ★
感谢参与,应助指数 +1
wizardfan: 金币+2, 应助指数-1, 谢谢参与,不过回帖更倾向于交流,而不是应助 2013-04-23 07:05:01
第一、什么物种
第二、什么叫“析出”
总有一天我要修改用户名。
2楼2013-04-22 09:37:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erlei

木虫 (正式写手)

引用回帖:
2楼: Originally posted by 137167741 at 2013-04-22 09:37:32
第一、什么物种
第二、什么叫“析出”

双壳纲
析出是就是将十二个编码基因从线粒体全基因组序列中截取下来
3楼2013-04-22 10:24:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

genomelin

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 272求调剂 +9 材紫有化 2026-02-28 9/450 2026-03-02 20:22 by hypershenger
[考研] 0856材料求调剂 +12 hyf hyf hyf 2026-02-28 13/650 2026-03-02 20:19 by hypershenger
[考研] 085600材料工程一志愿中科大总分312求调剂 +9 吃宵夜1 2026-02-28 11/550 2026-03-02 20:14 by hypershenger
[考研] 材料284求调剂,一志愿郑州大学英一数二专硕 +15 想上岸的土拨鼠 2026-02-28 15/750 2026-03-02 20:13 by hypershenger
[考研] 理学,工学,农学调剂,少走弯路,这里欢迎您! +4 likeihood 2026-03-02 7/350 2026-03-02 20:13 by 崔wj
[考研] 材料类考研调剂 +3 gemmgemm 2026-03-01 4/200 2026-03-02 20:02 by hypershenger
[考研] 302材料工程求调剂 +5 Doleres 2026-03-01 6/300 2026-03-02 19:53 by 张晓芳0105
[考研] 材料085601调剂 +5 多多子. 2026-03-02 5/250 2026-03-02 19:15 by zhukairuo
[考研] 289求调剂 +8 yang婷 2026-03-02 9/450 2026-03-02 19:08 by zhukairuo
[考研] 289求调剂 +3 BrightLL 2026-03-02 3/150 2026-03-02 17:28 by 杨杨杨紫
[考研] 一志愿东北大学材料专硕328,求调剂 +3 shs1083 2026-03-02 3/150 2026-03-02 17:27 by houyaoxu
[考研] 0856化工专硕求调剂 +15 董boxing 2026-03-01 15/750 2026-03-02 15:06 by 晃晃不许晃
[考研] 26考研报考西工大材料308分求调剂 +4 weizhong123 2026-03-01 4/200 2026-03-02 12:46 by 无际的草原
[考研] 292求调剂 +7 yhk_819 2026-02-28 7/350 2026-03-02 12:43 by 无际的草原
[考研] 0856材料与化工,270求调剂 +8 YXCT 2026-03-01 9/450 2026-03-02 11:01 by 无际的草原
[考研] 264求调剂 +4 巴拉巴拉根556 2026-02-28 4/200 2026-03-02 10:48 by yuchj
[考研] 322求调剂 +3 熊境喆 2026-03-01 3/150 2026-03-02 08:44 by houyaoxu
[考研] 化工299分求调剂 一志愿985落榜 +5 嘻嘻(*^ω^*) 2026-03-01 5/250 2026-03-01 19:47 by 无际的草原
[硕博家园] 2025届双非化工硕士毕业,申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
信息提示
请填处理意见