| 24小時(shí)熱門版塊排行榜 |
| 9 | 1/1 | 返回列表 |
| 查看: 2509 | 回復(fù): 8 | |||
幽幽之蘭鐵蟲(chóng) (初入文壇)
|
[交流]
【求助/交流】懇求:如何進(jìn)行序列拼接 已有6人參與
|
|
請(qǐng)問(wèn)各位如何將測(cè)得的序列進(jìn)行拼接啊,要是有詳細(xì)的圖文解析那就更感謝啦, O(∩_∩)O~ |
銀蟲(chóng) (小有名氣)
|
此名為ContigExpress的軟件可用于做序列拼接,主要使用方法如下: 1. 解壓縮下載的壓縮文件contig.zip文件,保證文件CExpress.exe,Gexudat.def在同一個(gè)目錄下,打開(kāi)Cexpress.exe應(yīng)用程序,進(jìn)入ContigExpress操作界面,如圖1。 圖1 2. 點(diǎn)擊菜單上的“Project”選擇“Add Fragments”,一般我們發(fā)給您的是AB1文件,如果您有其它格式的文件,也可以選擇,在這里我們選擇AB1文件,以其為例,如圖2。 圖2 3. 選擇您存放AB1文件(即我們Email給您的測(cè)序結(jié)果的彩圖文件)的目錄,選擇文件類型為ALL FILES, 之后打開(kāi)要拼接的AB1,從而添加進(jìn)ContigExpress軟件。在此以A、B兩個(gè)序列為例,如果有多個(gè)序列的也可以同時(shí)添加進(jìn)入。 圖3 4. 選中要拼接的序列,再選菜單 “Assemble”欄下的“Assemble Selected Fragments”命令,或用工具欄上的按鈕,如圖3。若兩個(gè)結(jié)果能夠拼接起來(lái)的,會(huì)得到一個(gè)Assemble1下的contig1的結(jié)果,如圖4。 圖4 5. 雙擊contig1,打開(kāi)拼接后的結(jié)果,選中菜單“VIEW”欄,進(jìn)入VIEW OPTION,將SHOW ALIGNMENT AS 由 TEXT 改為 GRAPH.,點(diǎn)擊OK 后得到結(jié)果如圖5。此時(shí)可能會(huì)因?yàn)閮蓷l序列的測(cè)序結(jié)果誤差,會(huì)有不同的地方,在拼接圖片框中的綠色豎杠就表示了這些不同的地方,如圖所示。接著可點(diǎn)擊綠色豎杠找到有誤差的地方,進(jìn)行修改。 6. 在修改過(guò)程中,遇到有誤差的地方,可以根據(jù)峰形來(lái)判斷是多讀還是漏讀來(lái)進(jìn)行修改,此時(shí)電腦認(rèn)為是漏讀堿基的地方會(huì)以點(diǎn)來(lái)表示,如圖5,此處很明顯是A序列上多讀了一個(gè)G堿基,可將其刪除。(注:因?yàn)檐浖旧淼膯?wèn)題,只有在拼接過(guò)程中是正向的序列才能進(jìn)行修改操作,若在反向上修改堿基,保存時(shí)會(huì)產(chǎn)生錯(cuò)誤而直接關(guān)閉程序。所以若要修改反向序列上的堿基,可先保存后,把原有的Assemble1的結(jié)果拆開(kāi),點(diǎn)序列圖標(biāo)上的“Name”,如圖3,所選中的序列上的一個(gè)“name”橫欄,使序列按Name的升降次序來(lái)排列,把要作為正向的序列放到要作為反向序列上面即可。以此序列為例,將其改變方向后可實(shí)現(xiàn)反向拼接,如圖6。此時(shí)反向拼接中,電腦也會(huì)將整條序列自動(dòng)進(jìn)行反向互補(bǔ)。若有多條序列只要讓反向序列在正向序列上即可。建議最好是修改一個(gè)堿基后就保存一下,然后拆開(kāi)了再拼接,這樣可以最大程度地避免誤差。) 圖5 圖6 7. 全部修改完成后,所有的綠色豎杠就消失了,但電腦對(duì)那些個(gè)別堿基不同,即沒(méi)有多讀或漏讀的現(xiàn)象,不會(huì)用綠色豎杠標(biāo)出,這樣也就是說(shuō),在整個(gè)序列中還有誤差存在,所以這時(shí)我們就需要用到工具欄上的另一個(gè)按鈕“find next ambiguous”,如圖7,來(lái)找序列中因兩條序列堿基不同而出現(xiàn)的N值,需要記住此時(shí)也只能在正鏈上進(jìn)行修改。如圖,此處依據(jù)峰形可確定A鏈上的C應(yīng)該為A,把它校正即可。 圖7 8. 經(jīng)過(guò)校正和修改之后,兩條序列的重疊部分已經(jīng)完全一致,序列已經(jīng)拼通,可以把整個(gè)拼接結(jié)果保存下來(lái),若想拷貝并保存整條序列,可以右擊contig1選取“Select All”,再選“Copy sequence from … bp to … bp”,如圖8。然后打開(kāi)可存放文本的應(yīng)用程序,如記事本,把整個(gè)序列粘貼進(jìn)去,保存即可。 圖8 |

銀蟲(chóng) (小有名氣)

鐵蟲(chóng) (初入文壇)
鐵蟲(chóng) (初入文壇)
金蟲(chóng) (正式寫(xiě)手)
金蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
銀蟲(chóng) (初入文壇)

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +7 | zjl的號(hào) 2026-03-16 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 【0856】化學(xué)工程(085602)313 分,本科學(xué)科評(píng)估A類院;瘜W(xué)工程與工藝,誠(chéng)求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
[考研] 274求調(diào)劑0856材料化工 +13 | z2839474511 2026-03-11 | 14/700 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +4 | 梨花珞晚風(fēng) 2026-03-17 | 4/200 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] 341求調(diào)劑 +3 | 番茄頭--- 2026-03-10 | 3/150 |
|
|
[考研] 332求調(diào)劑 +3 | zjy101327 2026-03-11 | 6/300 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|
|
[考博] 2026年博士申請(qǐng) +3 | QwQwQW10 2026-03-11 | 3/150 |
|