歡迎關注”生信修煉手冊”!
在mirdeep軟件的分析結果中觅够,會提供miRNA前體的二級結構纳令,這個結果實際上是通過調用RNAfold
來實現(xiàn)的香拉,該軟件是一個經(jīng)典的預測RNA二級結構的軟件穗熬,網(wǎng)址如下
http://rna.tbi.univie.ac.at//cgi-bin/RNAWebSuite/RNAfold.cgi
官網(wǎng)提供了在線服務,只需要上傳fasta格式的序列即可内狗,示意如下
默認參數(shù)會輸出以下兩種二級結構
1. optimal secondary structure
最佳二級結構穆壕,保證對應的自由能最小,最小自由能簡稱MFE, 結果示意如下
2. centroid secondary structure
自由能表征改變這個結構需要注入的能量大小其屏,對應的數(shù)值越小喇勋,該結構越穩(wěn)定。
同時給出了可視化結果偎行,示意如下
這個程序也是可以下載到本地運行的川背,基本用法如下
RNAfold < hsa.hairpin.fa -noPS > precursors.str
-noPS
參數(shù)代表不產(chǎn)生二級結構對應的postscript文件,這種文件可以轉換為PDF格式蛤袒,產(chǎn)生的str
文件內容如下
>hsa-let-7a-1
UGGGAUGAGGUAGUAGGUUGUAUAGUUUUAGGGUCACACCCACCACUGGGAGAUAACUAUACAAUCUACUGUCUUUCCUA
(((((.(((..((((((((((((((((...(((.....))).((....))....))))))))))))))))..)))))))) (-35.60)
采用dot-bracket表示法標記二級結構熄云,上述用法只給出了最佳的二級結構預測結果和對應的自由能。
·end·
—如果喜歡妙真,快分享給你的朋友們吧—
掃描關注微信號缴允,更多精彩內容等著你!