組裝基本流程_第1頁
組裝基本流程_第2頁
組裝基本流程_第3頁
全文預覽已結束

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領

文檔簡介

重測序序列組裝流程圖參考基因1、bwaindex參考基因庫 目的基因3、bwasampe生成sam文件4、Samtoolsview生成bam文件1f5、samtoolssort生成排過序后的bam文件6、samtoolsmpileup8、vcfutils.plvcf2fq生成內容為fastq格式的vcf文件用戶名:student密碼:1234561、 bwaindex就是給參考序列建索引,參數(shù)用的是默認參數(shù)。bwaindexref.fa2、 bwaaln是把目標序列和參考基因序列進行比對bwaaln-t6ref.fareads.fq>reads.sai3、 bwasampe比對pair-endreads,輸出SAM格式文件。因為我們的基本上是pair-endreads,所以用sampe。如果是比對single-endreads,就要用bwasamse了。用的參數(shù)用的是默認參數(shù)。reads_2.saireads_1.fqbwasamperef.fareads_1.saireads_2.saireads_1.fqreads_2.fq>reads.sam4、samtoolsview把sam文件轉換成bam文件。samtoolsview-bSreads.sam-oreads.bam5、samtoolssort對bam文件進行排序,默認按染色體名稱和位置大小順序。samtoolssortreads.bamreads.sort6、 samtoolsmpileup生成位點信息文件,可以通過不同的參數(shù)來的到不同的結果,mapping用的參數(shù)是-uS,之后找過snp,用的參數(shù)是-DSugf。把bam文件轉換成bcf格式文件。samtoolsmpileup-usreads.sort.bam>reads.bcf7、 bcftoolsview是把bcf文件輸出為vcf格式。bcftoolsview-cgreads.bcf>reads.vcf8、vcfutils.plvcf2fq把結果的一致性序列整理成

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論