24小时热门版块排行榜    

查看: 2727  |  回复: 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的回帖

mycaas

铁杆木虫 (正式写手)

【答案】应助回帖

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

xiaokoupu

银虫 (小有名气)

【答案】应助回帖

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

tw7649116

木虫 (小有名气)

【答案】应助回帖

目前我也在做水稻叶绿体基因组组装,原始数据如下是水稻基因组重测序数据(rawdata.fastq.gz),根据相关文献已经进行了以下分析:
利用bwa得到mapping到reference的.bam文件
利用samtools从bam文件提取mapped reads.bam(20M左右)
利用soapdenovo或者velvet进行组装获得了scaffold序列,但是只有100k不到,而水稻叶绿体基因组大约130多k
下面就不知道怎么做了,也不知道前面的步骤对不对!

楼主问题解决了吗?是怎么做的呢?
6楼2014-04-28 18:50:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nn219

新虫 (初入文坛)

我们是先把原始数据做一个质控,然后用SOAPdenovo或CLC软件拼接成contig或scafford,然后找参考序列,挑出叶绿体序列,再手动拼接,gap一般可以用不同的拼接软件试下看能不能补齐。

另外,楼主可以告诉我SOAPdenovo的具体操作么
7楼2015-01-13 10:28:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赖苏亮

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by nn219 at 2015-01-13 10:28:44
我们是先把原始数据做一个质控,然后用SOAPdenovo或CLC软件拼接成contig或scafford,然后找参考序列,挑出叶绿体序列,再手动拼接,gap一般可以用不同的拼接软件试下看能不能补齐。

另外,楼主可以告诉我SOAPdeno ...

你好!
      我现在也在做叶绿体基因组拼接,我的思路跟你的很相似,请问如果先拼接然后再筛选出叶绿体序列,会不会由于软件的误差将核的跟线粒体的序列也拼进去,而导致拼接错误?
      我看许多文献都是先用一个近缘物种的叶绿体基因组作为reference筛选叶绿体reads,然后再拼接的。
8楼2016-06-05 17:13:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

筱杨i

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by 赖苏亮 at 2016-06-05 17:13:48
你好!
      我现在也在做叶绿体基因组拼接,我的思路跟你的很相似,请问如果先拼接然后再筛选出叶绿体序列,会不会由于软件的误差将核的跟线粒体的序列也拼进去,而导致拼接错误?
      我看许多文献都是先用 ...

能请问一下是怎么筛选的叶绿体基因组么

发自小木虫Android客户端
9楼2019-03-05 16:23:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cgdsbsd 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂 +15 于海海海海 2026-03-19 15/750 2026-03-19 22:41 by 学员8dgXkO
[考研] 085600材料与化工 求调剂 +14 enenenhui 2026-03-13 15/750 2026-03-19 21:53 by allen-yin
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 求调剂 +3 Mqqqqqq 2026-03-19 3/150 2026-03-19 14:11 by peike
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 0703化学调剂,求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 一志愿华中科技大学,080502,354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 085410人工智能专硕317求调剂(0854都可以) +3 xbxudjdn 2026-03-18 3/150 2026-03-18 22:14 by zhq0425
[考研] 330求调剂 +3 小材化本科 2026-03-18 3/150 2026-03-18 21:55 by 无懈可击111
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 298-一志愿中国农业大学-求调剂 +7 手机用户 2026-03-17 7/350 2026-03-18 14:34 by vgtyfty
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
信息提示
请填处理意见