24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2503  |  回复: 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的回帖

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的回帖
相关版块跳转 我要订阅楼主 cgdsbsd 的主题更新
信息提示
请填处理意见