24小时热门版块排行榜    

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

erlei

木虫 (正式写手)

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

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

» 猜你喜欢

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

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

erlei

木虫 (正式写手)

引用回帖:
4楼: Originally posted by genomelin at 2013-04-23 07:06:38
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 ...

谢谢!
5楼2013-04-23 08:47:53
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 跳槽后在研项目怎么办? +4 简单化xn 2026-08-22 8/400 2026-08-23 07:20 by ilovexiaobin
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见