確定RNA seq鏈特異性

2021-08-27 16:10:12 字數 1587 閱讀 5204

當需要確定鏈特異性時,可如下操作,需要安裝bedops, rseqc兩個程式。

安裝bedops, rseqc

conda install bedops

conda install rseqc

將比對過程中使用的gtf檔案轉為12列的bed檔案

在bedops中有gtf2bed可以將gft轉換為12列的bed檔案。

輸入命令:

gtf2bed < mus_musculus.grcm38.93.chr.gtf > mouse.bed

其中mus_musculus.grcm38.93.chr.gtf為輸入檔名稱,mouse.bed為輸出檔名。發現會報錯如下:

error: potentially missing gene or transcript id from gtf attributes (malformed gtf at line [1]?)

原因在連線中:

該錯誤表明第一行缺少該transcript_id欄位。gene_id正如你所注意到的,它有乙個字段,但沒有transcript_id欄位。gtf 2.2規範指示此字段是必填字段,但其值可以是空字串。

解決方案我選擇第二個:

輸入命令:

awk '' mus_musculus.grcm38.93.chr.gtf  > mus_musculus.grcm38.93.chr.change.gtf

然後重新執行命令:

gtf2bed < mus_musculus.grcm38.93.chr.change.gtf > mouse.bed

順利生成12列的mouse.bed檔案。

確定鏈特異性

輸入命令:

infer_experiment.py -r mouse.bed -i *.bam

結果:

結果解釋在網頁:

這個結果說this is single-end, strand specific rna-seq data. strandness of reads are concordant with strandness of reference gene.

建庫方式是fr-secondstrand。

參考:

此圖1++,1--,2+-,2-+建庫方式為: fr-secondstrand

此圖1+-,1-+,2++,2--建庫方式為: fr-firststrand

css選擇器特異性計算

翻譯自w3c 中的一段 是關於css選擇器優先順序的 計算乙個選擇器的特異性 選擇器的特異性計算如下 1.計算選擇器中id選擇器的數量 a 2.計算選擇器中class選擇器,屬性選擇器和偽類選擇器的數量 b 3.計算型別選擇器和偽元素選擇器的數量 c 4.忽視通用選擇器 連線三個數字a b c就是特...

敏感性 特異性,查準率 查全率

總結 敏感性 查全率 召回率 查全率又稱召回率 查準率即精度。準確率與之均不同,定義如下 被正確診斷的案例之和與總案例之比即為準確率。在癌症示例中,敏感性和特異性指 查準率和查全率的定義如下 從這裡可以看出,敏感性就是查全率,但特異性並不是查準率。敏感性和特異性是這個矩陣中的行。更具體地說,如果我們...

CSS選擇器優先順序(特異性)

1 行內樣式 important 2 id選擇器 id3 class 屬性 偽類選擇器 title,input type text hover4 型別和偽元素選擇器 div,before 0,0,0,0 行內樣式,id選擇器,class 屬性 偽類選擇器 型別 偽元素 權重相同,定義靠後優先 con...