| 9 | 1/1 | 返回列表 |
| 查看: 2747 | 回復(fù): 8 | ||
cgdsbsd新蟲 (初入文壇)
|
[求助]
葉綠體序列拼接 已有1人參與
|
|
我現(xiàn)在有基因組數(shù)據(jù)、轉(zhuǎn)錄組數(shù)據(jù),怎么提取出葉綠體基因,并且拼接成一個(gè)圓環(huán)? 我開始是想先用一個(gè)從NCBI上下載的一個(gè)相近物種的葉綠體基因做reference,然后用soap2把基因組序列mapping到reference上,篩選出可能是葉綠體的短read,用soapdenovo將這些read拼接成contig,然后用BWA mapping在reference上,看有多少gap ,在gap兩端設(shè)計(jì)引物做PCR,將gap補(bǔ)滿,不知道這個(gè)思路對(duì)不對(duì)? 現(xiàn)在我得出的短read都是單序列,沒有成對(duì)的,然后從頭合成拼接的時(shí)候生成2500多條contig,而且最長的才300多bp,沒有scafer生成,下一步不知道該怎么繼續(xù)了,很茫然,請(qǐng)各位大牛幫忙,謝謝。。。 |

鐵桿木蟲 (正式寫手)
銀蟲 (小有名氣)
新蟲 (初入文壇)

新蟲 (初入文壇)

木蟲 (小有名氣)
|
目前我也在做水稻葉綠體基因組組裝,原始數(shù)據(jù)如下是水稻基因組重測(cè)序數(shù)據(jù)(rawdata.fastq.gz),根據(jù)相關(guān)文獻(xiàn)已經(jīng)進(jìn)行了以下分析: 利用bwa得到mapping到reference的.bam文件 利用samtools從bam文件提取mapped reads.bam(20M左右) 利用soapdenovo或者velvet進(jìn)行組裝獲得了scaffold序列,但是只有100k不到,而水稻葉綠體基因組大約130多k 下面就不知道怎么做了,也不知道前面的步驟對(duì)不對(duì)! 樓主問題解決了嗎?是怎么做的呢? |
新蟲 (初入文壇)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +5 | 某某某某位 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +7 | 紙魚ly 2026-03-21 | 8/400 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 11/550 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +17 | chaojifeixia 2026-03-19 | 19/950 |
|
|
[考研] 327求調(diào)劑 +5 | prayer13 2026-03-23 | 5/250 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +6 |  ̄^ ̄゜汗 2026-03-19 | 9/450 |
|
|
[考研] 一志愿中國石油大學(xué)(華東) 本科齊魯工業(yè)大學(xué) +4 | 石能偉 2026-03-17 | 4/200 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 268求調(diào)劑 +9 | 簡單點(diǎn)0 2026-03-17 | 9/450 |
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|