使用mirDeep2進行miRNA seq資料分析

2021-09-29 03:43:58 字數 1868 閱讀 4245

git clone  mirdeep2.0.1.2

cd mirdeep2.0.1.2/

使用install.pl指令碼進行安裝

perl install.pl
會有如下的提示資訊

提示資訊

可以按照他的要求,直接使用source ~/.bashrc載入環境變數,然後再次執行perl install.pl就會幫你解決依賴關係,依賴工具如下

bowtie short read aligner

vienna package with rnafold

squid library

randfold

perl package pdf::api2

參考基因組的fasta檔案

mirbase中該物種的成熟mirna

mrbase中該物種的前體mirna

高通量測序結果的fasta檔案

假如你已經有了如下檔案

檔名描述資訊

cel_cluster.fa參考基因組的fasta檔案

mature_ref_this_species.famirbase中該物種的成熟mirna

mature_ref_other_species.famirbase中該物種鄰近物種的成熟mirna

precursors_ref_this_species.famrbase中該物種的前體mirna

reads.fa高通量測序結果的fasta檔案

第一步: 建立索引

bowtie-build cel_cluster.fa cel_cluster
第二步: 將read回帖到參考基因組

-s reads_collapsed.fa -t reads_collapsed_vs_genome.arf -v

各個引數的含義如下:

第三步(可選): 快速進行定量。如果不需要**新的mirna, 可以用直接用mirbase資料庫進行定量

quantifier.pl -p precursors_ref_this_species.fa -m mature_ref_this_species.fa \

-r reads_collapsed.fa -t cel -y 16_19

輸出結果為mirna_expressed.csv, 記錄每個樣本的每個mirna的count數,結果同樣可以用網頁開啟expression_16_19.html檢視

第四步: 鑑定新的mirna,並進行定量

mirdeep2.pl reads_collapsed.fa cel_cluster.fa reads_collapsed_vs_genome.arf \

mature_ref_this_species.fa mature_ref_other_species.fa \

precursors_ref_this_species.fa -t c.elegans 2> report.log

這一步要求的參考基因組的序列不能有'atcgn'以外的字元,沒遇到報錯就萬事大吉,遇到報錯就用tr解決吧

第五步: 瀏覽結果

最後可以開啟results.html檢視結果。

使用C 和Excel進行報表開發(2)

內容摘要 下面的 使用乙個現有的統計圖excel檔案,載入這個檔案,設定某一項的值和統計圖的表示方式,然後儲存到另乙個excel檔案裡。如果是vs2005 office2003,也可以選擇使用visual studio tools for office來作。下面的 使用乙個現有的統計圖excel檔案...

使用Retrofit2進行HTTP請求設定請求超時

採用retrofit2本身可以進行優雅的restful請求,但是無法設定請求超時時間,需要配合okhttp3來設定請求超時.新增依賴 compile com.squareup.okhttp3 logging interceptor 3.0.1 新增配置資訊並設定超時時間 connecttimeout...

使用WinSock2 SPI進行網路控制訪問內容

編者按 與傳統的包過濾防火牆技術不同,本文從應用層閘道器技術入手,深入 了利用winsock2 spi進行網路內容訪問控制的問題。這是網路安全的一項新內容,或者說,它為網路安全技術的愛好者和研發人員提供了乙個新的思路。防火牆可以實施和執行網路訪問策略,但是,傳統的防火牆技術集中於如何防範外部網路對內...