24小时热门版块排行榜    

查看: 3376  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cgdsbsd

新虫 (初入文坛)

[求助] 叶绿体序列拼接 已有1人参与

我现在有基因组数据、转录组数据,怎么提取出叶绿体基因,并且拼接成一个圆环啊?
我开始是想先用一个从NCBI上下载的一个相近物种的叶绿体基因做reference,然后用soap2把基因组序列mapping到reference上,筛选出可能是叶绿体的短read,用soapdenovo将这些read拼接成contig,然后用BWA mapping在reference上,看有多少gap  ,在gap两端设计引物做PCR,将gap补满,不知道这个思路对不对?
现在我得出的短read都是单序列,没有成对的,然后从头合成拼接的时候生成2500多条contig,而且最长的才300多bp,没有scafer生成,下一步不知道该怎么继续了,很茫然,请各位大牛帮忙,谢谢。。。
回复此楼

» 猜你喜欢

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

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

xiaokoupu

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
呵呵,我觉得是不是还可以确定些保守基因,然后利用这些基因通过相似性钓一些short reads,因为mappin的mismatch一般<2,用基因去钓取可以适当放宽些条件,然后把钓取的序列添加进去再拼接;没这么做过,如果lz有兴趣这么做,完了告诉我下效果怎么样啊,哈哈
3楼2012-07-28 21:53:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

mycaas

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
gap可以先用软件补一补,不一定要用PCR
没有pair-end信息做scaffold确实挺困难
2楼2012-07-27 17:06:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cgdsbsd

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by mycaas at 2012-07-27 17:06:54
gap可以先用软件补一补,不一定要用PCR
没有pair-end信息做scaffold确实挺困难

xie  xie   
求知。。。
4楼2012-07-30 10:07:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cgdsbsd

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by xiaokoupu at 2012-07-28 21:53:02
呵呵,我觉得是不是还可以确定些保守基因,然后利用这些基因通过相似性钓一些short reads,因为mappin的mismatch一般<2,用基因去钓取可以适当放宽些条件,然后把钓取的序列添加进去再拼接;没这么做过,如果lz有 ...

怎么确定保守基因?hai   you  soapdenovo    de   config_file   zen   me   she   zhi   can  shu   a   ?   xie   xie  zhi   jiao  ...
求知。。。
5楼2012-07-30 11:19:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见