| 5 | 1/1 | 返回列表 |
| 查看: 2515 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
幽幽之蘭鐵蟲 (初入文壇)
|
[交流]
【求助/交流】懇求:如何進(jìn)行序列拼接 已有6人參與
|
||
|
請問各位如何將測得的序列進(jìn)行拼接啊,要是有詳細(xì)的圖文解析那就更感謝啦, O(∩_∩)O~ |
銀蟲 (小有名氣)

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

鐵蟲 (初入文壇)
鐵蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +6 | yy要上岸呀 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +3 | 睿08 2026-03-17 | 5/250 |
|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么! 2026-03-15 | 18/900 |
|
|
[考研] 293求調(diào)劑 +6 | 世界首富 2026-03-11 | 6/300 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動化所實(shí)習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[基金申請]
有必要更換申報口嗎
20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
[考研] 一志愿中科院,化學(xué)方向,295求調(diào)劑 +4 | 一氧二氮 2026-03-11 | 4/200 |
|
|
[考研] 308求調(diào)劑 +5 | 是Lupa啊 2026-03-11 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 考研調(diào)劑 +4 | 芬達(dá)46 2026-03-12 | 4/200 |
|